Improved fracs parameterizations for cross-sections of isotopes near the proton-drip line in projectile fragmentation reactions

Improved fracs parameterizations for cross-sections of isotopes near the proton-drip line in projectile fragmentation reactions 1 first-author 1 1 2 corresp Corresponding author, machunwang@126.com Corresponding author, machunwang@126.com machunwang@126.com 3 Institute of Particle and Nuclear Physics, Henan Normal University, Xinxiang 453007, China Institute of Particle and Nuclear Physics, Henan Normal University, Xinxiang 453007, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China College of Information and Electronic Engineering, Shangqiu Institute of Technology, Shangqiu 476000, China College of Information and Electronic Engineering, Shangqiu Institute of Technology, Shangqiu 476000, China The fracs parameterizations, labelled as fracs-c, have been improved in order to predict the presence of isotopes near the proton drip line produced in projectile fragmentation reactions. By investigating the cross-sections for proton-rich isotopes in a series of reactions with energies ranging from intermediate to relativistic, it is shown that the fracs-c parameterizations can predict isotopes near the proton drip line considerably well. The fracs-c parameterizations are suggested to serve as an effective tool for predicting the presence of proton-rich isotopes with large asymmetry in a projectile fragmentation reaction. Different reactions have been investigated to check these results. Drip line isotope Proton-rich isotope fracs Projectile fragmentation 1 introduction 1 It has been predicted that more than 8,000 nuclides should be bounded, which have lifetimes longer than 1 microsecond [ 1 1 - 4 4 1-4 ]. Most of these predicted nuclides are radioactive, and those with greatest asymmetry form the boundaries of the nuclides chart. Isotopes near the boundaries of the nuclides chart have attracted significant attention from both theoretical and experimental scientists as their properties remain unclear. Within the intermediate mass range, the study of isotopes near the proton drip line is very attractive, because radioactive nuclear beam (RNB) facilities have opened new research opportunities. Some exotic phenomena have been illustrated in proton-rich isotopes, such as the proton skin, proton halo, core deformation [ 5 5 , 6 6 ], and ( β -)delayed one or two proton emissions [ 7 7 - 12 12 7-12 ]. At present, the new generation of RNB facilities will extend the research to the more extremely asymmetric isotopes near the drip line region [ 13 13 - 16 16 13-16 ]. Improved techniques also provide the new opportunities to learn about their properties in greater detail. The production probability, i.e., cross-section, is one of the most important physical quantities to study proton-rich nuclides in reactions while the production low rates make them difficult to measure. In experiments, it is important to be able to predict the yields or cross-sections for proton-rich isotopes, which are not well-explained in the existing methods. Many methods have been proposed to predict isotopic yields based on empirical formulas [ 17 17 - 19 19 17-19 ]. Recently, Mei proposed the fracs parameterizations [ 20 20 ] based on ideas from the epax parameterizations [ 21 21 - 23 23 21-23 ] by incorporating an extra term that depends on the incident energy. These parameterizations predict significantly better results compared with those from the existing models. Though many advantages have been illustrated, the fracs parameterizations are currently not the best option for predicting the presence of proton-rich isotopes, especially for those near the proton drip line. The main reasons are that, the neutron-rich projectile nuclei, as well as the neutron-rich fragments, have been elaborately treated. However, the discussion for the proton-rich fragments (with neutron excess I ≡ N - Z < -1) have been omitted. In this article, the fracs parameterizations are improved by re-adjusting the parameters for proton-rich isotopes, which we call the fracs-c parameters, and we show that the parameters predict more accurate cross-sections for isotopes near the proton drip lines. This article is organized as follows. Methods are described in Sec. 2. The results are presented and discussed in Sec. 3. A summary is presented in Sec. 4. 2 Methods 1 Readers are referred to the original description of the fracs parameterizations [ 20 20 ]. In this article, we would like to introduce the most important portion related to the modifications applied to the original description. Three parts have been included in the fracs parameterizations to predict the cross-section of an isotope in projectile fragmentation reactions, i.e., the mass yield, isobaric (charge) distribution, and odd-even staggering. The cross-section for an isotope is described by Ref. [ 20 20 ], <math display="block" id="M1"><mrow><mi>σ</mi><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>Z</mi><mo stretchy="false">)</mo><mo>=</mo><mi>Y</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mrow><mtext>prob</mtext></mrow></msub><mo>−</mo><mi>Z</mi><mo stretchy="false">)</mo><msub><mi>Δ</mi><mrow><mtext>OES</mtext></mrow></msub><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>Z</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where Y ( A ) is the mass yield and Y ( Z prob - Z ) represents the corresponding isobaric distribution. Δ OES ( A,Z ) is adopted to deal with odd-even staggering. The isobaric distribution is assumed to have a Gaussian distribution, with the remaining nuclear generation being related to the statistical evaporation from excited pre-fragments. The isobaric distribution of fragments with a given mass number A is borrowed from the epax parameterizations [ 21 21 - 23 23 21-23 ], <math display="block" id="M2"><mrow><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mrow><mtext>prob</mtext></mrow></msub><mo>−</mo><mi>Z</mi><mo stretchy="false">)</mo><mo>=</mo><mi>n</mi><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mi>R</mi><mo>|</mo><mo stretchy="false">(</mo><msub><mi>Z</mi><mrow><mtext>prob</mtext></mrow></msub><mo>−</mo><mi>Z</mi><msup><mrow><mo stretchy="false">)</mo><mo>|</mo></mrow><mi>U</mi></msup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where Z prob is the most probable charge, n is a normalization factor. and R is the width parameter. In fracs, ideas from the epax3 parameterizations have been borrowed [ 23 23 ], which allowed improved predictions for proton-rich isotopes. The projectile asymmetry is also corrected by introducing <math id="M3"><mrow><msubsup><mi>Δ</mi><mtext>m</mtext><mtext>p</mtext></msubsup></mrow></math> and <math id="M4"><mrow><msubsup><mi>Δ</mi><mtext>m</mtext><mtext>n</mtext></msubsup></mrow></math> terms, and so on. Extra exponent parameters U p and a transition point Z exp ( A ) are introduced for the proton-rich isotopes. Efforts have been also expended to overcome the disadvantages of the epax3 parameterizations, in which the dependence on incident energy was not considered [ 23 23 ]. In the following results, it is shown that the fracs parameterizations are still not sufficient for predicting the presence of proton-rich isotope of large asymmetry. To overcome this disadvantage of fracs, the Z exp term, which is the only term that affects the very proton rich isotopes, should be adjusted. The Z exp reads as: <math display="block" id="M5"><mrow><msub><mi>Z</mi><mrow><mtext>exp</mtext></mrow></msub><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>Z</mi><mrow><mtext>prob</mtext></mrow></msub><mo>+</mo><mo stretchy="false">[</mo><msub><mi>l</mi><mn>1</mn></msub><mo>+</mo><msub><mi>l</mi><mn>2</mn></msub><msup><mrow><mo stretchy="false">(</mo><mi>A</mi><mo>/</mo><mn>2</mn><mo stretchy="false">)</mo></mrow><mrow><mn>0.3</mn></mrow></msup><mo stretchy="false">]</mo><mi>ln</mi><mo stretchy="false">(</mo><mn>10</mn><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><mi>R</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where R is a parameter to describe the width isobaric distribution, <math display="block" id="M6"><mrow><mi>R</mi><mo>=</mo><msub><mi>R</mi><mn>0</mn></msub><mi>exp</mi><mo stretchy="false">[</mo><msub><mi>r</mi><mn>1</mn></msub><mi>A</mi><mo>+</mo><msub><mi>r</mi><mn>2</mn></msub><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><mo stretchy="false">(</mo><msub><mi>r</mi><mn>3</mn></msub><mo>+</mo><msub><mi>Z</mi><mtext>p</mtext></msub><mo>−</mo><msub><mi>Z</mi><mrow><mi>β</mi><mtext>p</mtext></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msub><mi>r</mi><mn>4</mn></msub><mo>+</mo><msub><mi>r</mi><mn>5</mn></msub><mi>A</mi><mo>/</mo><msub><mi>A</mi><mtext>p</mtext></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mn>,</mn></mrow></math> where r i is parameters, which are the same for the symmetric and asymmetric projectile nuclei. R 0 is different for the proton rich projectile ( <math id="M7"><mrow><msubsup><mi>R</mi><mn>0</mn><mtext>p</mtext></msubsup></mrow></math> ) and neutron rich projectile ( <math id="M8"><mrow><msubsup><mi>R</mi><mn>0</mn><mtext>n</mtext></msubsup></mrow></math> ), respectively. For a proton rich projectile, <math id="M9"><mrow><msubsup><mi>R</mi><mn>0</mn><mtext>p</mtext></msubsup></mrow></math> reads: <math display="block" id="M10"><mrow><msubsup><mi>R</mi><mn>0</mn><mtext>p</mtext></msubsup><mo>=</mo><msub><mi>r</mi><mn>0</mn></msub><mi>exp</mi><mo stretchy="false">[</mo><msub><mi>r</mi><mn>7</mn></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mtext>p</mtext></msub><mo>−</mo><msub><mi>Z</mi><mrow><mi>β</mi><mtext>p</mtext></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mn>,</mn></mrow></math> In Eq. (3), l 1 , l 2 , and the width parameter R are fixed by fitting a large amount of experimental data in the present fracs. l 2 is the parameter to adjust the mass of isotope. The relevant parameter to adjust <math id="M11"><mrow><msubsup><mi>R</mi><mn>0</mn><mi>p</mi></msubsup></mrow></math> is r 7 . Thus l 1 and r 7 should be the key parameters to improve the fracs for the proton-rich isotopes. As the description in Ref. [ 20 20 ], r 7 adjusts the width of the quasi-Gaussian shape of the isotopic distribution, and l 1 tunes the transition point where the quasi-Gaussian shape turns into an exponential decay, respectively (more parameters are used together in fracs) [ 20 20 ]. With a series of experimental results being investigated (as the following results in Sec. 3), r 7 and l 1 are refitted for the proton-rich isotopes near the proton drip line, which are listed in Table 1 Table 1 and shown in Fig. 1 Fig. 1 . With the decreasing I , the isospin dependence of the proton rich isotope is also not well considered. An extra term is introduced to describe the dependence of the cross-section on the asymmetry of very proton-rich isotopes with I < -2: 1 1 Parameter 1 1 E (MeV/u) 1 1 Value 1 1 r 7 1 1 / 1 1 0.19 1 1 l 1 1 1 140 1 1 1.2 1 1   1 1 345 1 1 2.0 1 1   1 1 650 1 1 0.4 1 1   1 1 1000 1 1 0.2 1 1   1 1 1015 1 1 0.2 <math display="block" id="M12"><mrow><msup><mi>σ</mi><mo>′</mo></msup><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>Z</mi><mo stretchy="false">)</mo><mo>=</mo><mi>σ</mi><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>Z</mi><msup><mrow><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow><mi>a</mi></msup><mn>,</mn><mi/><mi/><mi/><mi/><mi/><mi>a</mi><mo>=</mo><mo>|</mo><mi>I</mi><mo>|</mo><mo>−</mo><mn>3</mn><mi/><mtext> </mtext><mi>i</mi><mi>f</mi><mtext> </mtext><mi/><mi>I</mi><mo><</mo><mo>−</mo><mn>2,</mn></mrow></math> The fitting result (see Fig. 1 Fig. 1 ) to the correlation between l 1 and E shows that it obeys a Gaussian function, which reads <math id="M13"><mrow><msub><mi>l</mi><mn>1</mn></msub><mo>=</mo><mn>0.200</mn><mo>+</mo><mfrac><mrow><mn>730.878</mn></mrow><mrow><mn>318.283</mn><msqrt><mrow><mi>π</mi><mo>/</mo><mn>2</mn></mrow></msqrt></mrow></mfrac><mo>×</mo><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mn>2</mn><mo stretchy="false">(</mo><mfrac><mrow><mi>E</mi><mo>−</mo><mn>315.108</mn></mrow><mrow><mn>318.283</mn></mrow></mfrac><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></math> . 3 results and discussion 1 In this section, the predictions for the proton-rich isotopes by frac-c will be compared to the epax2 and fracs methods. The measured cross-sections for Ge ( Z = 32) and Se ( Z = 34) isotopes in the 140 MeV/u 78 Kr + 9 Be reaction [ 24 24 ] are plotted in Fig. 2 Fig. 2 . The predictions from the epax2 [ 22 22 ], fracs [ 20 20 ], and fracs-c parameterizations are compared to the measured data. The fracs parameterizations show improved predictions for these isotopes with I ranging from -4 to 0, although the measured results are overestimated. By adopting the new fitted parameters and the extra term (shown in Table 1 Table 1 ), the fracs-c can predict the measured cross-sections quite well. More measured isotopes in projectile fragmentation reactions will be tested with the fracs-c parameterizations. In Fig. 3 Fig. 3 , measured isotopes with I =-4 and -5 in the 345 MeV/u 78 Kr + 9 Be [ 11 11 ] are plotted together with predictions from the fracs and fracs-c parameterizations. One can see that the fracs parameterizations significantly overestimate the measured results, while the results by the fracs-c parameterizations show better agreement with the measured data compared to the fracs parameterizations. The cross-sections of very proton-rich isotopes with Z = 21–28 and I ≥ -6 in the 650 MeV/u 58 Ni + 9 Be reaction have been reported [ 25 25 ]. The cross-sections of these measured isotopes are shown in Fig. 4 Fig. 4 , where the predicted results by fracs and fracs-c parameterizations are compared. For the I = 0, 1, and 2 nuclides, fracs can reproduce the measured results well. While the fracs tends to underestimate the measured results as the nuclide become more proton-rich. The fracs-c parameterizations reproduce the results for I ≥ 0 and I < 0 isotopes, even when I = -6. The cross-sections of the I =-1 and -2 nuclides have been measured in the 1.0 GeV/u 124 Xe + Pb [ 26 26 ]and 136 Xe + Pb, and 1.015 GeV/u 112 Sn + Be [ 27 27 ]reactions. The measured and predicted cross-sections by fracs and fracs-c parameterizations are compared in Fig. 5 Fig. 5 . One can see that, in the three reactions, the fracs predictions agree with the measured data for the I =-1 nuclide, while the results for the I =-2 nuclide are underestimated. The fracs-c parameterizations improve the prediction quality and reproduce the cross-sections well both for the I =-1 and -2 nuclides. It has been claimed that the fracs parameterizations are more appropriate for reactions with incident energy higher than 140 MeV/u because the parameters are fitted for the series of reactions with higher energy. In this article, we have investigated the reaction systems with intermediate (Ni + Be) to the heavy ones (Xe + Pb) masses with incident energy above 140 MeV/u. It is shown that the fracs-c parameterizations improve the quality of the fracs parameterizations for proton-rich isotopes and reproduce the measured results quite well. To quantitatively estimate the deviation factors for the predictions, an A -factor technique is adopted to show the deviation between the predicted and measured results in this work. The A -factor is defined as [ 28 28 ] <math display="block" id="M14"><mrow><mi>A</mi><mo>≡</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><mfrac><mrow><mo>|</mo><msubsup><mi>σ</mi><mi>i</mi><mrow><mtext>exp</mtext></mrow></msubsup><mo>−</mo><msubsup><mi>σ</mi><mi>i</mi><mrow><mtext>cal</mtext></mrow></msubsup><mo>|</mo></mrow><mrow><mo>|</mo><msubsup><mi>σ</mi><mi>i</mi><mrow><mtext>exp</mtext></mrow></msubsup><mo>+</mo><msubsup><mi>σ</mi><mi>i</mi><mrow><mtext>cal</mtext></mrow></msubsup><mo>|</mo></mrow></mfrac><mn>,</mn></mrow></math> <math display="block" id="M15"><mrow><msubsup><mi>Δ</mi><mi>i</mi><mrow><mtext>exp</mtext></mrow></msubsup><mo>→</mo><msub><mi>Δ</mi><mi>i</mi></msub><mo>≡</mo><msqrt><mrow><msup><mrow><mo stretchy="false">(</mo><msubsup><mi>Δ</mi><mi>i</mi><mrow><mtext>exp</mtext></mrow></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><msubsup><mi>Δ</mi><mi>i</mi><mrow><mtext>cal</mtext></mrow></msubsup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mn>,</mn></mrow></math> where <math id="M16"><mrow><msubsup><mi>Δ</mi><mi>i</mi><mrow><mtext>cal</mtext></mrow></msubsup></mrow></math> should be the statistical error of the predicted cross-section by the fracs-c parameterizations. Statistical error has been ignored in this work since it could not be calculated by fracs-c. The A -factor is independent and invariant with respect to exchange of the experimental and calculated cross-sections. A prediction will be more accurate for smaller A -factor values. The predicted A -factor results using the fracs and fracs-c parameterizations are shown in Fig. 6 Fig. 6 . It is easy to see that the values of the A -factor for the fracs-c parameterizations, in general smaller than 0.3, are quite smaller than those for the fracs parameterizations. The improved fracs-c parameterizations for fracs make it more suitable for cross-section predictions regarding proton-rich isotopes near the drip-line area, which fulfills the urgent requirement for experimental designation of rare nuclei in projectile fragmentation reactions. 4 summary 1 The fracs parameterizations have been improved, which we call fracs-c, and have been used to predict the cross-sections for isotopes produced in projectile fragmentation reactions. The parameters ( r 7 and l 1 in fracs) used to determine the cross-sections of proton-rich isotopes have been re-adjusted, and an extra term is also introduced to better reproduce the cross-sections for proton-rich isotopes. By analyzing the measured cross-sections of proton-rich isotopes produced in a series of reactions, including the 140 and 345 MeV/u 78 Kr + 9 Be, 650 MeV/u 58 Ni + 9 Be, 1 GeV/u 124, 136 Xe + Pb, 1.015 GeV/u 112 Sn + 9 Be reactions, it is shown that fracs-c can better reproduce the measured cross-sections of nuclides with I =-6, which is considerably close to the proton-drip line of intermediate-mass nuclides. References X.W. Xia , Y. Lim , P.W. Zhao et al ., The limits of the nuclear landscape explored by the relativistic continuum Hartree-Bogoliubov theory . Atomic Data and Nuclear Data Tables , 121-122 , 1 ( 2018 ) doi: 10.1016/j.adt.2017.09.001 http://doi.org/10.1016/j.adt.2017.09.001 . The limits of the nuclear landscape explored by the relativistic continuum Hartree-Bogoliubov theory H. Wang , G. Audi , F.G. Kondev et al ., The AME2016 atomic mass evaluation (II) . Tables, graphs and references. Chin. Phys. 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) R. Utama , J. Piekarewicz , Validating neural-network refinements of nuclear mass models . Phys. Rev. C 97 , 014306 ( 2018 ). doi: 10.1103/PhysRevC.97.014306 http://doi.org/10.1103/PhysRevC.97.014306 Validating neural-network refinements of nuclear mass models N. Wang , M. Liu , X. Wu et al ., Surface diffuseness correction in global mass formula . Phys. Lett. B 734 , 215 ( 2014 ). doi: 10.1016/j.physletb.2014.05.049 http://doi.org/10.1016/j.physletb.2014.05.049 For WS4 mass table see http://www.imqmd.com/mass/WS4.txt http://www.imqmd.com/mass/WS4.txt Surface diffuseness correction in global mass formula G. Sawhney , R.K. Gupta , M.K. Sharma , Importance of Deformations in Dynamical Evolution of Proton-halo Nuclei . Acta Physica Polonica B 47 , 959 ( 2016 ). doi: 10.5506/APhysPolB.47.959 http://doi.org/10.5506/APhysPolB.47.959 Importance of Deformations in Dynamical Evolution of Proton-halo Nuclei D.Q. Fang , W. Guo , C.W. Ma et al ., Examining the exotic structure of the proton-rich nucleus Al-23 . Phys. Rev. C 76 , 031601 ( 2007 ). doi: 10.1103/PhysRevC.76.031601 http://doi.org/10.1103/PhysRevC.76.031601 Examining the exotic structure of the proton-rich nucleus Al-23 L.J. Sun , C.J. Lin , X.X. Xu et al ., Experimental Study of Beta-Delayed Proton Emission of Ca-(36,37) . Chin. Phys. Lett . 32 , 012301 ( 2015 ). doi: 10.1088/0256-307X/32/1/012301 http://doi.org/10.1088/0256-307X/32/1/012301 Experimental Study of Beta-Delayed Proton Emission of Ca-(36,37) Y.G. Ma , D.Q. Fang , X.Y. Sun et al ., Different mechanism of two-proton emission from proton-rich nuclei Al-23 and Mg-22 . Phys. Lett. B 743 , 306 ( 2015 ). doi: 10.1016/j.physletb.2015.02.066 http://doi.org/10.1016/j.physletb.2015.02.066 Different mechanism of two-proton emission from proton-rich nuclei Al-23 and Mg-22 X.X. Xu , L.J. Sun , C.J. Lin et al ., Observation of Beta-delayed two-proton emission in the decay of Si-22 . Phys. Lett. B 766 , 312 ( 2017 ). doi: 10.1016/j.physletb.2017.01.028 http://doi.org/10.1016/j.physletb.2017.01.028 Observation of Beta-delayed two-proton emission in the decay of Si-22 L.J. Sun , X.X. Xu , D.Q. Fang et al ., Beta-decay study of the Tz=-2 proton-rich nucleus Mg-20 . Phys. Rev. C 95 , 014314 ( 2017 ). doi: 10.1103/PhysRevC.95.014314 http://doi.org/10.1103/PhysRevC.95.014314 Beta-decay study of the Tz=-2 proton-rich nucleus Mg-20 B. Blank , T. Goigoux , P. Ascher et al ., New neutron-deficient isotopes from Kr-78 fragmentation . Phys. Rev. C 93 , 061301(R) ( 2016 ). doi: 10.1103/PhysRevC.93.061301 http://doi.org/10.1103/PhysRevC.93.061301 New neutron-deficient isotopes from Kr-78 fragmentation T. Goigoux , P. Ascher , B. Blank et al ., Two-Proton Radioactivity of Kr-67 . Phys. Rev. Lett . 117 , 162501 ( 2016 ). doi: 10.1103/PhysRevLett.117.162501 http://doi.org/10.1103/PhysRevLett.117.162501 Two-Proton Radioactivity of Kr-67 T. Kubo , In-flight RI beam separator BigRIPS at RIKEN and elsewhere in Japan . NIM B 204 , 97 ( 2003 ). doi: 10.1016/S0168-583X(02)01896-7 http://doi.org/10.1016/S0168-583X(02)01896-7 In-flight RI beam separator BigRIPS at RIKEN and elsewhere in Japan B.-H. Sun , J.-W. Zhao , X.-H. Zhang et al ., Towards the full realization of the RIBLL2 beam line at the HIRFL-CSR complex . Science Bulletin 63 , 78 ( 2018 ). doi: 10/1016/j.scib.2017.12.005 http://doi.org/10/1016/j.scib.2017.12.005 Towards the full realization of the RIBLL2 beam line at the HIRFL-CSR complex J.W. Xia , W.L. Zhan , B.W. Wei et al ., The heavy ion cooler-storage-ring project (HIRFL-CSR) at Lanzhou . NIM A 488 , 11 ( 2002 ). doi: 10.1016/S0168-9002(02)00475-8 http://doi.org/10.1016/S0168-9002(02)00475-8 The heavy ion cooler-storage-ring project (HIRFL-CSR) at Lanzhou M. Hausmann , A.M. Aaron , A.M. Amthor et al ., Design of the Advanced Rare Isotope Separator ARIS at FRIB . NIM B 317 , 349 ( 2013 ). doi: 10.1016/j.nimb.2013.06.042 http://doi.org/10.1016/j.nimb.2013.06.042 Design of the Advanced Rare Isotope Separator ARIS at FRIB C.W. Ma , J.L. Xu , An empirical formula for isotopic yield in Fe + p spallation reactions . J. Phys. G: Nucl. Part. Phys . 44 , 125101 ( 2017 ). doi: 10.1088/1361-6471/aa90e6 http://doi.org/10.1088/1361-6471/aa90e6 An empirical formula for isotopic yield in Fe + p spallation reactions M. Yu , H.-L. Wei , Y.-D. Song , C.-W. Ma , Chin. Phys. C 41 , 094001 ( 2017 ). Chin. Phys. C M. B. Tsang , W. G. Lynch , W. A. Friedman et al ., Fragmentation cross sections and binding energies of neutron-rich nuclei . Phys. Rev. C 76 , 041302(R) ( 2007 ). doi: 10.1103/PhysRevC.76.041302 http://doi.org/10.1103/PhysRevC.76.041302 Fragmentation cross sections and binding energies of neutron-rich nuclei B. Mei , Improved empirical parameterization for projectile fragmentation cross sections . Phys. Rev. C 95 , 034608 ( 2017 ). doi: 10.1103/PhysRevC.95.034608 http://doi.org/10.1103/PhysRevC.95.034608 Improved empirical parameterization for projectile fragmentation cross sections K. Sümmerer , B. Blank , Modified empirical parametrization of fragmentation cross sections . Phys. Rev. C 61 , 034607 ( 2000 ). doi: 10.1103/PhysRevC.61.034607 http://doi.org/10.1103/PhysRevC.61.034607 Modified empirical parametrization of fragmentation cross sections K. Sümmerer , Improved empirical parametrization of fragmentation cross sections . Phys. Rev. C 86 , 014601 ( 2012 ). doi: 10.1103/PhysRevC.86.014601 http://doi.org/10.1103/PhysRevC.86.014601 Improved empirical parametrization of fragmentation cross sections K. Sümmerer , Erratum: Improved empirical parametrization of fragmentation cross sections [ Phys. Rev. C 86 , 014601 ( 2012 )]. Erratum: Improved empirical parametrization of fragmentation cross sections K. Sümmerer , Erratum: Improved empirical parametrization of fragmentation cross sections Phys. Rev. C 87 , 039903 ( 2013 ). doi: 10.1103/PhysRevC.87.039903 http://doi.org/10.1103/PhysRevC.87.039903 Erratum: Improved empirical parametrization of fragmentation cross sections A. Stolz , T. Baumann , N.H. Frank et al ., First observation of Ge-60 and Se-64 . Phys. Lett. B 627 , 32 ( 2005 ). doi: 10.1016/j.physletb.2005.08.130 http://doi.org/10.1016/j.physletb.2005.08.130 First observation of Ge-60 and Se-64 B. Blank , S. Andriamonje , R.D. Moral et al ., Production cross sections and the particle stability of proton-rich nuclei from Ni-58 fragmentation . Phys. Rev. C 50 , 2398 ( 1994 ). doi: 10.1103/PhysRevC.50.2398 http://doi.org/10.1103/PhysRevC.50.2398 Production cross sections and the particle stability of proton-rich nuclei from Ni-58 fragmentation D. Henzlova , K.-H. Schmidt , M. V. Ricciardi et al ., Experimental investigation of the residues produced in the 136Xe + Pb and 124Xe + Pb fragmentation reactions at 1A GeV . Phys. Rev. C 78 , 044616 ( 2008 ). doi: 10.1103/PhysRevC.78.044616 http://doi.org/10.1103/PhysRevC.78.044616 Experimental investigation of the residues produced in the 136Xe + Pb and 124Xe + Pb fragmentation reactions at 1A GeV A. Stolz , T. Faestermann , J. Friese et al ., Projectile fragmentation of Sn-112 at Elab=1 A GeV . Phys. Rev. C 65 , 064603 ( 2002 ). doi: 10.1103/PhysRevC.65.064603 http://doi.org/10.1103/PhysRevC.65.064603 Projectile fragmentation of Sn-112 at Elab=1 A GeV S.K. Sharma , B. Kamys , F. Goldenbaum et al ., Ranking and validation of spallation models for isotopic production cross sections of heavy residua . Eur. Phys. J. A 53 , 150 ( 2017 ). doi: 10.1140/epja/i2017-12334-2 http://doi.org/10.1140/epja/i2017-12334-2 Ranking and validation of spallation models for isotopic production cross sections of heavy residua We thank Dr. B. Mei for sending us his well written code for fracs, which cuts the time for us to finish this work. 10.1007/s41365-018-0439-4 This work was partially supported by the National Natural Science Foundation of China (No. U1732135), the Key Research Program of Frontier Sciences of CAS (No. QYZDJSSW-SLH002), and the Natural and Science Foundation in Henan Province (No. 162300410179). 2018-04-12 2018-05-15 2018-05-17 2018-05-29 2018-07-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-0439-4 http://dx.doi.org/10.1007/s41365-018-0439-4 . 2018 Nuclear Science and Techniques 07 29