Production of heavy neutron-rich nuclei with radioactive beams in multinucleon transfer reactions

Production of heavy neutron-rich nuclei with radioactive beams in multinucleon transfer reactions 1 2 corresp first-author licheng@mail.bnu.edu.cn licheng@mail.bnu.edu.cn licheng@mail.bnu.edu.cn 1 2 1 2 1 2 1 2 1 2 3 corresp fszhang@bnu.edu.cn fszhang@bnu.edu.cn fszhang@bnu.edu.cn The Key Laboratory of Beam Technology and Material Modification of Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China The Key Laboratory of Beam Technology and Material Modification of Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China Beijing Radiation Center, Beijing 100875, China Beijing Radiation Center, Beijing 100875, China Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator of Lanzhou, Lanzhou 730000, China Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator of Lanzhou, Lanzhou 730000, China The production mechanism of heavy neutron-rich nuclei is investigated by using the multinucleon transfer reactions of 136,148 Xe+ 208 Pb and 238 U+ 208 Pb in the framework of a dinuclear system model. The evaporation residual cross sections of target-like fragments are studied with the reaction system 148 Xe+ 208 Pb at near barrier energies. The results show that the final isotopic production cross sections in the neutron-deficient side are very sensitive to incident energy while it is not sensitive in the neutron-rich side. Comparing the isotopic production cross sections for the reactions of 208 Pb bombarded with stable and radioactive projectiles, we find that neutron-rich radioactive beams can significantly increase the production cross sections of heavy neutron-rich nuclei. Multinucleon transfer reactions DNS model Heavy-ion collisions Neutron-rich nuclei 1 Introduction 1 Production of new heavy neutron-rich nuclei has attracted widespread interests in recent years [ 1 1 - 11 11 1-11 ]. The production of these exotic nuclei located far from the stability line is very important for exploring their structure and decay properties [ 12 12 , 13 13 ], as well as a measurement of spins, neutron skin thickness [ 14 14 ], magnetic dipole, and electric quadrupole moments and charge radii, et al . In addition, it is also extremely important for understanding the r-process in nuclear astrophysics. The r-process path is located (and interrupted by fission) just in the region of unknown heavy nuclei with a large neutron excess. The neutron closed shell N = 126 is probably the last “waiting point” in the r-process of nucleosynthesis [ 15 15 ]. It is helpful for answering the origin of heavy elements from iron to uranium in the Universe. Many new neutron-rich nuclei are produced by experiment with fragmentation method of projectile [ 16 16 , 17 17 ]. For example, there are 60 new neutron-rich isotopes in the atomic number range of 60 ≤ Z ≤ 78 produced by using a 1 GeV/nucleon 238 U beam impinged on a 9 Be target by Kurcewicz et al. at GSI Darmstadt [ 16 16 ]. They found that the fragmentation plays a dominant role in the production of the isotopes of the elements above Z = 72. However, production cross sections of these nuclei with the fragmentation method are very small and are about a few ten pb orders of magnitude. It is difficult for further production of more neutron-rich nuclei with the fragmentation method. Recently, a multinucleon transfer (MNT) reaction experiment of 136 Xe+ 198 Pt was performed by Watanabe et al. at GANIL in Caen [ 9 9 ]. They found that the MNT approach has huge advantages for the production of very neutron-rich nuclei with Z ⩽77 compared to the above fragmentation approach. It demonstrates that the MNT route is a more promising path to produce new neutron-rich isotopes around N =126 shell. Radioactive ion beam facilities can provide very neutron-rich projectiles [ 18 18 , 19 19 ] which could significantly improve the production cross section of heavy neutron-rich nuclei [ 20 20 , 21 21 ]. Some facilities, such as the high intensity heavy ion accelerator facility (HIAF) in China [ 22 22 , 23 23 ], have a plan to produce new heavy neutron-rich nuclei with intense secondary beams of exotic radioactive nuclei by the MNT reactions. Theoretical support for these very time-consuming and expensive experiments is vital for choosing the optimum target-projectile-energy combinations and for the estimation of cross sections. The dinuclear system (DNS) model is a semiclassical model which has shown reasonable success in predicting the evaporation residual cross section in fusion reactions [ 24 24 - 26 26 24-26 ] and multinucleon transfer reactions [ 27 27 , 28 28 ]. In this work, we apply the DNS model to study the multinucleon transfer process with stable projectiles 136 Xe, 238 U, and radioactive projectile 148 Xe impinged on a 208 Pb target. The structure of this paper is as follows. In Sect. 2, we briefly introduce the DNS model. In Sect. 3, the isotope production cross sections of 136 Xe+ 208 Pb, 238 U+ 208 Pb, and 148 Xe+ 208 Pb systems have been studied by using the DNS model. Finally, the conclusion is given in Sect. 4. 2 Theoretical descriptions 1 In the framework of the DNS model, the distribution probability can be obtained by solving a set of master equations numerically in the potential energy surface of the DNS. The probability distribution function P ( Z 1 , N 1 , E 1 , t ) for fragment 1 with proton number Z 1 , neutron number N 1 and excitation energy E 1 is described by the following master equation [ 27 27 , 29 29 ]: <math display="block" id="M1"><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mi>d</mi><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mrow/></mtd></mtr></mtable><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup></mrow></munder></mstyle><msub><mi>W</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mo>;</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">[</mo><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mi>P</mi><mo stretchy="false">(</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>E</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>−</mo><msub><mi>d</mi><mrow><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup></mrow></munder></mstyle><msub><mi>W</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mo>;</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">[</mo><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msubsup><mi>E</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo>−</mo><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup></mrow></msub><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>−</mo><mo stretchy="false">[</mo><msubsup><mi>Λ</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow><mrow><mi>q</mi><mi>f</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>Θ</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>Λ</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow><mrow><mi>f</mi><mi>i</mi><mi>s</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>Θ</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mn>.</mn></math> The sum is taken over all possible proton and neutron numbers that fragments <math id="M2"><mrow><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup></mrow></math> and <math id="M3"><mrow><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup></mrow></math> may take, but only one nucleon transfer is considered in the model with the relation <math id="M4"><mrow><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mo>=</mo><msub><mi>Z</mi><mn>1</mn></msub><mo>±</mo><mn>1</mn></mrow></math> and <math id="M5"><mrow><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup><mo>=</mo><msub><mi>N</mi><mn>1</mn></msub><mo>±</mo><mn>1</mn></mrow></math> . <math id="M6"><mrow><msub><mi>W</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mo>;</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow></math> denotes the mean transition probability from channel ( Z 1 , N 1 , E 1 ) to ( <math id="M7"><mrow><msubsup><mi>z</mi><mn>1</mn><mo>′</mo></msubsup></mrow></math> , N 1 , <math id="M8"><mrow><msubsup><mi>E</mi><mn>1</mn><mo>′</mo></msubsup></mrow></math> ), which can be read as <math display="block" id="M9"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>W</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mo>;</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><msub><mi>τ</mi><mrow><mi>m</mi><mi>e</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mo>;</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>E</mi><mn>1</mn><mo>'</mo></msubsup><mo stretchy="false">)</mo></mrow><mrow><msup><mi>ℏ</mi><mn>2</mn></msup><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><msub><mi>d</mi><mrow><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>×</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><msup><mi>i</mi><mo>'</mo></msup></mrow></munder></mstyle><mo>|</mo><mo>〈</mo><msubsup><mi>Z</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>E</mi><mn>1</mn><mo>'</mo></msubsup><mn>,</mn><msup><mi>i</mi><mo>'</mo></msup><mo>|</mo><mi>V</mi><mo>|</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>i</mi><mo>〉</mo><msup><mo>|</mo><mn>2</mn></msup><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where, τ mem and V denote memory time and single-particle potential, respectively. i denotes all remaining quantum numbers [ 30 30 ]. <math id="M10"><mrow><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow></math> denotes the microscopic dimension corresponding to the macroscopic state ( Z 1 , N 1 , E 1 ) <math display="block" id="M11"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>d</mi><mo stretchy="false">(</mo><msub><mi>m</mi><mn>1</mn></msub><mn>,</mn><msub><mi>m</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo>=</mo><mrow><mo>(</mo><mrow><mtable><mtr><mtd><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>m</mi><mn>1</mn></msub></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mtable><mtr><mtd><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>m</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> The n 1 = g 1 Δ ε and n 2 = g 2 Δ ε denote the numbers of valence states of fragments 1 and 2, respectively. g 1 and g 2 are mean single-particle level densities. <math id="M12"><mrow><mi>Δ</mi><mi>ε</mi><mo>=</mo><msqrt><mrow><mn>4</mn><msup><mi>ε</mi><mo>*</mo></msup><mo>/</mo><mi>g</mi></mrow></msqrt></mrow></math> is an energy range around Fermi energy, ε F , of both fragments. m 1 = n 1 /2 and m 2 = n 2 /2 denote corresponding numbers of valence nucleons [ 30 30 , 31 31 ]. The quasifission rate <math id="M13"><mrow><msubsup><mi>Λ</mi><mrow><msub><mrow/><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow><mrow><mi>q</mi><mi>f</mi></mrow></msubsup></mrow></math> and the fission rate <math id="M14"><mrow><msubsup><mi>Λ</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow><mrow><mi>f</mi><mi>i</mi><mi>s</mi></mrow></msubsup></mrow></math> of the heavy fragment are estimated with the one-dimensional Kramers formula [ 32 32 ]. <math id="M15"><mrow><mi>Θ</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><msqrt><mrow><msup><mi>ε</mi><mo>*</mo></msup><mo>/</mo><mi>a</mi></mrow></msqrt></mrow></math> is the local temperature. a = A i /12 is the level-density parameter. We assume that ε * is divided between the two primary fragments according to their mass ratio. Hence the excitation energy of fragment 1 is expressed as E 1 = A 1 /( A 1 + A 2 ) ε * . A 1 and A 2 are mass numbers of fragment 1 and fragment 2, respectively. The local excitation energy, ε * , is defined as <math display="block" id="M16"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msup><mi>ε</mi><mo>*</mo></msup><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><msub><mi>E</mi><mrow/></msub><mi>diss</mi><mo>−</mo><mo stretchy="false">[</mo><mi>U</mi><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo>−</mo><mi>U</mi><mo stretchy="false">(</mo><msub><mi>A</mi><mi>P</mi></msub><mn>,</mn><msub><mi>A</mi><mi>T</mi></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> Here E diss is the energy dissipation from the relative kinetic energy loss, which can be seen in Ref. [ 33 33 ]. The dissipation of the relative motion and angular momentum of the DNS is described by the Fokker-Planck equation. The U ( A 1 , A 2 ) and U ( A P , A T ) are the driving potentials at the configuration of ( A 1 , A 2 and at the entrance point of DNS, which can be read as <math display="block" id="M17"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mi>U</mi></mtd><mtd columnalign="left"><mrow><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mn>,</mn><mi>J</mi><mn>,</mn><mi>R</mi><mo>;</mo><msub><mi>β</mi><mn>1</mn></msub><mn>,</mn><msub><mi>β</mi><mn>2</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>1</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mi>B</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mo>+</mo><mi>B</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>2</mn></msub><mn>,</mn><msub><mi>N</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo>−</mo><mo stretchy="false">[</mo><mi>B</mi><mo stretchy="false">(</mo><mi>Z</mi><mn>,</mn><mi>N</mi><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>V</mi><mrow><mi>r</mi><mi>o</mi><mi>t</mi></mrow><mrow><mi>C</mi><mi>N</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>J</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><mi>V</mi><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mn>,</mn><mi>J</mi><mn>,</mn><mi>R</mi><mo>;</mo><msub><mi>β</mi><mn>1</mn></msub><mn>,</mn><msub><mi>β</mi><mn>2</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>1</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> Where, B ( Z 1 , N 1 ), B ( Z 2 , N 2 ), and B ( Z , N ) are the binding energies of the fragment 1, fragment 2, and the compound nucleus, respectively, in which the shell and the pairing corrections are included reasonably. <math id="M18"><mrow><msubsup><mi>V</mi><mrow><mi>r</mi><mi>o</mi><mi>t</mi></mrow><mrow><mi>C</mi><mi>N</mi></mrow></msubsup></mrow></math> is the rotation energy of the compound nucleus. The β 1 and β 2 are quadrupole deformation parameters of the two fragments at ground state which is static in the nucleon transfer process. θ 1 and θ 2 are angles between the collision orientations and the symmetry axes for fragment 1 and fragment 2, respectively. In Ref. [ 27 27 ], the calculations show that the tip-tip and side-side collisions give a similar result in the isotope distributions of 238 U+ 238 U collisions at incident energy near 1.1 V C . Hence in this work, only tip-tip cases are considered for calculations. The interaction potential between two fragments includes nuclear, Coulomb, and centrifugal parts <math display="block" id="M19"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>V</mi><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mn>,</mn><mi>J</mi><mn>,</mn><mi>R</mi><mo>;</mo><msub><mi>β</mi><mn>1</mn></msub><mn>,</mn><msub><mi>β</mi><mn>2</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>1</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><msub><mi>V</mi><mi>N</mi></msub><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mn>,</mn><mi>R</mi><mo>;</mo><msub><mi>β</mi><mn>1</mn></msub><mn>,</mn><msub><mi>β</mi><mn>2</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>1</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>V</mi><mi>C</mi></msub><mo stretchy="false">(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>,</mn><msub><mi>A</mi><mn>2</mn></msub><mn>,</mn><mi>R</mi><mo>;</mo><msub><mi>β</mi><mn>1</mn></msub><mn>,</mn><msub><mi>β</mi><mn>2</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>1</mn></msub><mn>,</mn><msub><mi>θ</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><mi>J</mi><mo stretchy="false">(</mo><mi>J</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><msup><mi>ℏ</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>μ</mi><msup><mi>R</mi><mn>2</mn></msup></mrow></mfrac><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where μ is the reduced mass of system. The nuclear potential is calculated by using the double-folding method. The Coulomb potential is obtained by Wong’s formula [ 31 31 ]. The cross sections of the primary fragments ( Z 1 , N 1 ) are calculated as follows <math display="block" id="M20"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>σ</mi><mrow><mtext>pr</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo stretchy="false">)</mo><mo>≈</mo></mrow></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mi>π</mi><msup><mi>ℏ</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>μ</mi><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>J</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>J</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mo stretchy="false">(</mo><mn>2</mn><mi>J</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mi>T</mi><mo stretchy="false">(</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mn>,</mn><mi>J</mi><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>×</mo><mo stretchy="false">[</mo><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>τ</mi><mrow><mtext>int</mtext></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><msub><mi>Y</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mo stretchy="false">]</mo><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where T ( E c.m. , J ) is the capture probability. P ( Z 1 , N 1 , τ int ) is the formation probability of primary fragments from the equilibrium state of DNS. The contribution of fission rate is ignored in this work. The <math id="M21"><mrow><msub><mi>Y</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub></mrow></math> denotes the primary yields of fragments from the quasifission rate [ 32 32 , 34 34 ], which is written as <math display="block" id="M22"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>Y</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msub><mi>τ</mi><mrow><mtext>int</mtext></mrow></msub></mrow></msubsup></mrow></mstyle><msubsup><mi>Λ</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub></mrow><mrow><mtext>qf</mtext></mrow></msubsup><mi>P</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>N</mi><mn>1</mn></msub><mn>,</mn><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><mi>d</mi><mi>t</mi><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> The interaction time, τ int , in the dissipation process of two colliding partners is dependent on the incident energy, E c.m. , in the center-of-mass frame and the angular momentum, J , which is calculated by using the deflection function method and has the value of a few 10 -20 s. For heavy systems with no potential pocket, the distance, R, between the centers of the two fragments is chosen to be the value at the touching point, in which the DNS is assumed to be formed. The maximal angular momentum is taken as the grazing collisions. More detailed descriptions of the DNS model can be found in Refs. [ 30 30 , 31 31 ]. To obtain the final isotope cross section distributions, we apply the GEMINI to deal with the subsequent de-excitation process of primary fragments. The GEMINI is a Monte Carlo code which follows the decay of a compound nuclei by a series of sequential binary decays until the resulting products are unable to undergo any further decay. It allows not only light-particle emission and symmetric fission, but also the decaying nucleus to emit a fragment of any mass [ 35 35 ]. Nuclear level densities are taken as a Fermi-gas form with default parameters. In this work, we simulate 1,000 times the deexcitation processes for each primary fragment. To test the DNS model for the description of multinucleon transfer reactions, the isotopic production cross sections for 136 Xe+ 208 Pb at E c.m. =450 MeV are calculated by using the DNS model. Figure 1 Figure 1 shows the isotopic production cross sections from Z =78 to 81. The dashed and solid lines denote the calculation results for primary and final fragments, respectively. The experimental data (solid circles) are from Ref. [ 2 2 ]. From Fig. 1 Fig. 1 , one can see that final yields shift to the left after the de-excitation process. The final isotopic cross sections are in good agreement with the experimental data. This indicates that the DNS model is applicable for the study of MNT reactions at near barrier energies. 3 Results and discussion 1 The radioactive neutron-rich nuclei have larger N/Z ratios. Such as 148 Xe is a known isotope of xenon with the largest neutron excess ( N/Z =1.74). It can be as a radioactive beam in the future. First, we discuss the isotopic cross sections in different incident energies. Figure 2 Figure 2 shows the primary isotopic cross sections of Pt at incident energies E c.m. =1.02, 1.06, and 1.10 V C for an 148 Xe+ 208 Pb system. V C is the Coulomb barrier. The solid, dash-dotted and dotted lines denote the calculation results from incident energy E c.m. =1.02, 1.06, and 1.10 V C , respectively. One can see that isotopic cross sections increase with the increasing of incident energies both in neutron-rich and neutron-deficient sides. Larger incident energy improves the transfer probability of nucleons on the same potential energy surface. However, larger incident energy will result in higher excitation energy for primary fragments. Figure 3 Figure 3 shows the final isotopic production cross sections from Z =75 to 80 for the reaction of 148 Xe+ 208 Pb at different incident energies. One can see that in the neutron-deficient side, the isotopic production cross sections increase with increasing incident energies. It is very sensitive to incident energies. There are two reasons for this. One is that isotopic production cross sections in the neutron-deficient side at a larger incident energy are larger for primary fragments. The other one is that a larger incident energy leads to primary fragments obtaining more excitation energy. More neutrons evaporating causes a shift of final yields to the neutron-deficient side. While in the neutron-rich side, the final isotopic production cross sections are not sensitive to incident energies, especially for extremely neutron-rich nuclei. In addition, the values of peak for a few proton stripping channels change very slowly. Because of this, the driving potential of the system is reduced in the configuration of mass symmetry. The transfer of nucleons from target to projectile is dominant near the entrance point of DNS. The N/Z ratio equilibration plays an important role between the projectile and target during the collision, which is one of the most important mechanisms to produce neutron-rich nuclei [ 36 36 ]. In Fig. 4 Fig. 4 , one shows the N =126 isotones distributions for the reactions of 136 Xe+ 208 Pb and 148 Xe+ 208 Pb at different incident energies. From Fig. 4(a) Fig. 4(a) , one can see that the cross sections of proton stripping channels are slightly lower than the pickup channels. The nucleus produced in the 208 Pb+ 136 Xe reaction with the largest neutron excess is only 204 Pt. It is because the N/Z value of 208 Pb (1.54) is larger than that of 136 Xe (1.52). It is hard to transfer protons from 208 Pb to 136 Xe. While in the Fig. 4(b) Fig. 4(b) , this phenomenon is opposite. The cross sections of proton stripping channels are obviously larger than the pickup channels. This reaction system can produce very neutron-rich nuclei (such as 199 Ta). Because the N/Z value of 148 Xe (1.74) is significantly larger than that of 208 Pb (1.54), the 208 Pb is easier to transfer protons to 148 Xe. In Fig. 5 Fig. 5 , we show the isotope production cross sections of the 208 Pb target bombarded with stable projectiles 136 Xe, 238 U, and radioactive projectile, 148 Xe. The dashed, solid, and dotted lines denote the calculation results from 136 Xe+ 208 Pb, 238 U+ 208 Pb, and 148 Xe+ 208 Pb systems at incident energy E c.m. =1.10 V C . The open circles denote unknown neutron-rich nuclei. The N/Z values of 136 Xe, 238 U, and 148 Xe are 1.52, 1.59 and 1.74, respectively. One can see that the isotope production cross sections in the neutron-rich side for each panel increase with the increasing N/Z value of projectile. Especially for larger proton stripping channels, this trend is more significant. In Fig. 5(a) and 5(b) Fig. 5(a) and 5(b) , there are 5 (for Z =75) and 3 (for Z =76) unknown neutron-rich nuclei produced by the 148 Xe+ 208 Pb reaction. This is because the exchange of particles between the two colliding partners is primarily driven by the N/Z ratio equilibration. The very neutron-rich nuclei are easily produced by the reaction of 148 Xe+ 208 Pb via neutrons transfer. It can be seen that the production cross sections for the new neutron-rich nuclei 200-204 Re are 1.710, 1.810, 0.173, 0.060, and 0.002 μ b, 204-206 Os are 1.040, 0.041, and 0.007 μ b, 206-208 Ir are 0.081, 0.006, and 0.001 μ b, 209,210 Pt are 0.012 and 0.004 μ b, and 211 Au is 0.003 μ b, respectively. Based on the above, the MNT reaction with the neutron-rich radioactive beams is a very effective method to produce new heavy neutron-rich nuclei. 4 Summary 1 In summary, the multinucleon transfer reactions of 136 Xe+ 208 Pb, 238 U+ 208 Pb, and 148 Xe+ 208 Pb are investigated using the DNS model. The isotopic production cross section of reaction 148 Xe+ 208 Pb at near barrier incident energies is studied. We find that the final isotopic production cross section of reaction 148 Xe+ 208 Pb in the neutron-deficient side is very sensitive to incident energies. For the neutron-rich nuclei, it is not sensitive to incident energies. The N/Z ratio equilibration plays an important role between the projectile and target during the collision. The cross sections of proton stripping channels are obviously larger than the pickup channels for the reaction of 148 Xe+ 208 Pb on N =126 shell closure. The extremely neutron-rich projectile, 148 Xe, improves significantly the production cross sections of neutron-rich nuclei. References C. Li , F. Zhang , J.J. Li et al ., Multinucleon transfer in the 136Xe+208Pb reaction . Phys. Rev. C 93 , 014618 ( 2016 ). doi: 10.1103/PhysRevC.93.014618 http://doi.org/10.1103/PhysRevC.93.014618 Multinucleon transfer in the 136Xe+208Pb reaction J.S. Barrett , W. Loveland , R. Yanez et al ., 136Xe+208Pb reaction: A test of models of multinucleon transfer reactions . Phys. Rev. C 91 , 064615 ( 2015 ). doi: 10.1103/PhysRevC.91.064615 http://doi.org/10.1103/PhysRevC.91.064615 136Xe+208Pb reaction: A test of models of multinucleon transfer reactions V.I. Zagrebaev , W. Greiner , Production of New Heavy Isotopes in Low-Energy Multinucleon Transfer Reactions . Phys. Rev. Lett . 101 , 122701 ( 2008 ). doi: 10.1103/PhysRevLett.101.122701 http://doi.org/10.1103/PhysRevLett.101.122701 Production of New Heavy Isotopes in Low-Energy Multinucleon Transfer Reactions W. Loveland , A.M. Vinodkumar , D. Peterson et al ., Synthesis of heavy nuclei using damped collisions: A test . Phys. Rev. C 83 , 044610 ( 2011 ). doi: 10.1103/PhysRevC.83.044610 http://doi.org/10.1103/PhysRevC.83.044610 Synthesis of heavy nuclei using damped collisions: A test V.I. Zagrebaev , W. Greiner , Production of heavy trans-target nuclei in multinucleon transfer reactions . Phys. Rev. C 87 , 034608 ( 2013 ). doi: 10.1103/PhysRevC.87.034608 http://doi.org/10.1103/PhysRevC.87.034608 Production of heavy trans-target nuclei in multinucleon transfer reactions V.I. Zagrebaev , W. Greiner , Production of heavy and superheavy neutron-rich nuclei in transfer reactions . Phys. Rev. C 83 , 044618 ( 2011 ). doi: 10.1103/PhysRevC.83.04461 http://doi.org/10.1103/PhysRevC.83.04461 Production of heavy and superheavy neutron-rich nuclei in transfer reactions L. Zhu , J. Su , F.S. Zhang , Influence of the neutron numbers of projectile and target on the evaporation residue cross sections in hot fusion reactions . Phys. Rev. C 93 , 064610 ( 2016 ). doi: 10.1103/PhysRevC.93.064610 http://doi.org/10.1103/PhysRevC.93.064610 Influence of the neutron numbers of projectile and target on the evaporation residue cross sections in hot fusion reactions K. Zhao , Z.X. Li , N. Wang et al ., Production mechanism of neutron-rich transuranium nuclei in 238U+238U collisions at near-barrier energies . Phys. Rev. C 92 , 024613 ( 2015 ). doi: 10.1103/PhysRevC.92.024613 http://doi.org/10.1103/PhysRevC.92.024613 Production mechanism of neutron-rich transuranium nuclei in 238U+238U collisions at near-barrier energies Y.X. Watanabe , Y.H. Kim , S.C. Jeong et al ., Pathway for the production of neutron-rich isotopes around the N=126 shell closure . Phys. Rev. Lett . 115 , 172503 ( 2015 ). doi: 10.1103/PhysRevLett.115.172503 http://doi.org/10.1103/PhysRevLett.115.172503 Pathway for the production of neutron-rich isotopes around the N=126 shell closure K. Sekizawa , K. Yabana , Time-dependent Hartree-Fock calculations for multinucleon transfer processes in 40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions . Phys. Rev. C 88 , 014614 ( 2013 ). doi: 10.1103/PhysRevC.88.014614 http://doi.org/10.1103/PhysRevC.88.014614 Time-dependent Hartree-Fock calculations for multinucleon transfer processes in 40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions K. Sekizawa , K. Yabana , Time-dependent Hartree-Fock calculations for multinucleon transfer and quasifission processes in the 64Ni+238U reaction . Phys. Rev. C 93 , 054616 ( 2016 ). doi: 10.1103/PhysRevC.93.054616 http://doi.org/10.1103/PhysRevC.93.054616 Time-dependent Hartree-Fock calculations for multinucleon transfer and quasifission processes in the 64Ni+238U reaction 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 C.W. Ma , C.Y. Qiao , T.T. Ding et al ., Temperature of intermediate mass fragments in simulated 40Ca + 40Ca reactions around the Fermi energies by AMD model . Nucl. Sci. Tech . 27 , 111 ( 2016 ). doi: 10.1007/s41365-016-0112-8 http://doi.org/10.1007/s41365-016-0112-8 Temperature of intermediate mass fragments in simulated 40Ca + 40Ca reactions around the Fermi energies by AMD model X.F. Li , D.Q. Fang , Y.G. Ma , Determination of the neutron skin thickness from interaction cross section and charge-changing cross section for B, C, N, O, F isotopes . Nucl. Sci. Tech . 27 , 71 ( 2016 ). doi: 10.1007/s41365-016-0064-z http://doi.org/10.1007/s41365-016-0064-z Determination of the neutron skin thickness from interaction cross section and charge-changing cross section for B, C, N, O, F isotopes R. Wang , L.W. Chen , Positioning the neutron drip line and the r-process paths in the nuclear landscape . Phys. Rev. C 92 , 031303(R) ( 2015 ). doi: 10.1103/PhysRevC.92.031303 http://doi.org/10.1103/PhysRevC.92.031303 Positioning the neutron drip line and the r-process paths in the nuclear landscape J. Kurcewicz , F. Farinon , H. Geissel et al ., Discovery and cross-section measurement of neutron-rich isotopes in the element range from neodymium to platinum with the FRS . Phys. Lett. B 717 , 371 ( 2012 ). doi: 10.1016/j.physletb.2012.09.021 http://doi.org/10.1016/j.physletb.2012.09.021 Discovery and cross-section measurement of neutron-rich isotopes in the element range from neodymium to platinum with the FRS O.B. Tarasov , D.J. Morrissey , A.M. Amthor et al ., Evidence for a Change in the Nuclear Mass Surface with the Discovery of the Most Neutron-Rich Nuclei with 17≤Z≤25 . Phys. Rev. Lett . 102 , 142501 ( 2009 ). doi: 10.1103/PhysRevLett.102.142501 http://doi.org/10.1103/PhysRevLett.102.142501 Evidence for a Change in the Nuclear Mass Surface with the Discovery of the Most Neutron-Rich Nuclei with 17≤Z≤25 J.M. D’Auria , A review of radioactive beam facilities in the world . Nucl. Instrum. Methods B 99 , 330 ( 1995 ). doi: 10.1016/0168-583X(95)80083-2 http://doi.org/10.1016/0168-583X(95)80083-2 A review of radioactive beam facilities in the world T. Motobayashi , H. Sakurai , Research with fast radioactive isotope beams at RIKEN . Prog. Theor. Exp. Phys . 2012 , 03C001 ( 2012 ). doi: 10.1093/ptep/pts059 http://doi.org/10.1093/ptep/pts059 Research with fast radioactive isotope beams at RIKEN W. Loveland , Synthesis of transactinide nuclei using radioactive beams . Phys. Rev. C 76 , 014612 ( 2007 ). doi: 10.1103/PhysRevC.76.014612 http://doi.org/10.1103/PhysRevC.76.014612 Synthesis of transactinide nuclei using radioactive beams M.H. Mun , G.G. Adamian , N.V. Antonenko et al ., Toward neutron-rich nuclei via transfer reactions with stable and radioactive beams . Phys. Rev. C 91 , 054610 ( 2015 ). doi: 10.1103/PhysRevC.91.054610 http://doi.org/10.1103/PhysRevC.91.054610 Toward neutron-rich nuclei via transfer reactions with stable and radioactive beams J.C. Yang , J.W. Xia , G.Q. Xiao et al ., High Intensity heavy ion Accelerator Facility (HIAF) in China . Nucl. Instrum. Methods B 317 , 263 ( 2013 ). doi: 10.1016/j.nimb.2013.08.046 http://doi.org/10.1016/j.nimb.2013.08.046 High Intensity heavy ion Accelerator Facility (HIAF) in China L.J. Mao , J.C. Yang , J.W. Xia et al ., Electron cooling system in the booster synchrotron of the HIAF project . Nucl. Instrum. Methods A 786 , 91 ( 2015 ). doi: 10.1016/j.nima.2015.03.052 http://doi.org/10.1016/j.nima.2015.03.052 Electron cooling system in the booster synchrotron of the HIAF project N. Wang , E.G. Zhao , W. Scheid et al ., Theoretical study of the synthesis of superheavy nuclei with Z=119 and 120 in heavy-ion reactions with trans-uranium targets . Phys. Rev. C 85 , 041601(R) ( 2012 ). doi: 10.1103/PhysRevC.85.041601 http://doi.org/10.1103/PhysRevC.85.041601 Theoretical study of the synthesis of superheavy nuclei with Z=119 and 120 in heavy-ion reactions with trans-uranium targets Z.Q. Feng , G.M. Jin , F. Fu et al ., Production cross sections of superheavy nuclei based on dinuclear system model . Nucl. Phys. A 771 , 50 ( 2006 ). doi: 10.1016/j.nuclphysa.2006.03.002 http://doi.org/10.1016/j.nuclphysa.2006.03.002 Production cross sections of superheavy nuclei based on dinuclear system model X.J. Bao , Y. Gao , J.Q. Li et al ., Possibilities for synthesis of new isotopes of superheavy nuclei in cold fusion reactions . Phys. Rev. C 93 , 044615 ( 2016 ). doi: 10.1103/PhysRevC.93.044615 http://doi.org/10.1103/PhysRevC.93.044615 Possibilities for synthesis of new isotopes of superheavy nuclei in cold fusion reactions Z.Q. Feng , G.M. Jin , J.Q. Li , Production of heavy isotopes in transfer reactions by collisions of 238U+238U . Phys. Rev. C 80 , 067601 ( 2009 ). doi: 10.1103/PhysRevC.80.067601 http://doi.org/10.1103/PhysRevC.80.067601 Production of heavy isotopes in transfer reactions by collisions of 238U+238U M.H. Mun , G.G. Adamian , N.V. Antonenko et al ., Production cross section of neutron-rich isotopes with radioactive and stable beams . Phys. Rev. C 89 , 034622 ( 2014 ). doi: 10.1103/PhysRevC.89.034622 http://doi.org/10.1103/PhysRevC.89.034622 Production cross section of neutron-rich isotopes with radioactive and stable beams L. Zhu , Z.Q. Feng , F.S. Zhang , Production of heavy neutron-rich nuclei in transfer reactions within the dinuclear system model . J. Phys. G 42 , 085102 ( 2015 ). doi: 10.1088/0954-3899/42/8/085102 http://doi.org/10.1088/0954-3899/42/8/085102 Production of heavy neutron-rich nuclei in transfer reactions within the dinuclear system model S. Ayik , Schürmann , W. Nörenberg , Microscopic transport theory of heavy-ion collisions I. basic relations and transport coefficients for symmetric fragmentation . Z. Phys. A 277 , 299 ( 1976 ). doi: 10.1007/BF01415605 http://doi.org/10.1007/BF01415605 Microscopic transport theory of heavy-ion collisions I. basic relations and transport coefficients for symmetric fragmentation Z.Q. Feng , G.M. Jin , J.Q. Li et al ., Formation of superheavy nuclei in cold fusion reactions . Phys. Rev. C 76 , 044606 ( 2007 ). doi: 10.1103/PhysRevC.76.044606 http://doi.org/10.1103/PhysRevC.76.044606 Formation of superheavy nuclei in cold fusion reactions G.G. Adamian , N.V. Antonenko , W. Scheid , Characteristics of quasifission products within the dinuclear system model . Phys. Rev. C 68 , 034601 ( 2003 ). doi: 10.1103/PhysRevC.68.034601 http://doi.org/10.1103/PhysRevC.68.034601 Characteristics of quasifission products within the dinuclear system model Z.Q. Feng , Production of neutron-rich isotopes around N = 126 in multinucleon transfer reactions . Phys. Rev. C 95 , 024615 ( 2017 ). doi: 10.1103/PhysRevC.95.024615 http://doi.org/10.1103/PhysRevC.95.024615 Production of neutron-rich isotopes around N = 126 in multinucleon transfer reactions G.G. Adamian , N.V. Antonenko , V.V. Sargsyan et al ., Possibility of production of neutron-rich Zn and Ge isotopes in multinucleon transfer reactions at low energies . Phys. Rev. C 81 , 024604 ( 2010 ). doi: 10.1103/PhysRevC.81.024604 http://doi.org/10.1103/PhysRevC.81.024604 Possibility of production of neutron-rich Zn and Ge isotopes in multinucleon transfer reactions at low energies 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 W. Królas , R. Broda , B. Fornal et al ., Dynamical deformation of nuclei in deep-inelastic collisions: A gamma coincidence study of 130Te + 275 MeV 64Ni and 208Pb + 345 MeV 58Ni heavy ion reactions . Nucl. Phys. A 832 , 170 ( 2010 ). doi: 10.1016/j.nuclphysa.2009.10.159 http://doi.org/10.1016/j.nuclphysa.2009.10.159 Dynamical deformation of nuclei in deep-inelastic collisions: A gamma coincidence study of 130Te + 275 MeV 64Ni and 208Pb + 345 MeV 58Ni heavy ion reactions 10.1007/s41365-017-0266-z This work was supported by the National Natural Science Foundation of China under Grants No. 11635003, No. 11025524 and No. 11161130520; the National Basic Research Program of China under Grant No. 2010CB832903; the European Commission’s 7th Framework Programme (Fp7-PEOPLE-2010-IRSES) under Grant Agreement Project No. 269131; the Project funded by China Postdoctoral Science Foundation (Grant No. 2016M600956). 2017-01-12 2017-04-12 2017-04-16 2017-06-29 2017-08-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. 2017 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-017-0266-z http://dx.doi.org/10.1007/s41365-017-0266-z . 2017 Nuclear Science and Techniques 08 28