Multinucleon transfer dynamics in nearly symmetric nuclear reactions

Multinucleon transfer dynamics in nearly symmetric nuclear reactions 1 first-author 2 3 1 1 corresp Corresponding author, fengzhq@scut.edu.cn Corresponding author, fengzhq@scut.edu.cn fengzhq@scut.edu.cn https://orcid.org/0000-0002-3440-5052 School of Physics and Optoelectronic Technology, South China University of Technology, Guangzhou 510641, China School of Physics and Optoelectronic Technology, South China University of Technology, Guangzhou 510641, China Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China University of Chinese Academy of Sciences, Beijing 100190, China University of Chinese Academy of Sciences, Beijing 100190, China Within the framework of the dinuclear system model, the multinucleon transfer dynamics for nearly symmetric nuclear collisions has been investigated. The reaction mechanism in the systems of 198 Pt+ 198 Pt and 204 Hg+ 198 Pt were investigated at beam energies around the Coulomb barrier. It was found that the isotopic yields are enhanced with increased incident energy in the domain of proton-rich nuclides. However, the production on the neutron-rich side weakly depends on the energy. The angular distribution with the beam energy was also analyzed in the multinucleon transfer reactions. Projectile-like fragments were produced towards the forward emission with increasing incident energy. The target-like fragments manifested the opposite trend in the transfer reactions. Dinuclear system model Multinucleon transfer reactions Symmetric nuclear collisions Angular distribution 1 Introduction 1 The synthesis of new isotopes in the superheavy region is a key issue, motivated by extending the shell closure, discovering a new mechanism for nuclear fission, waiting point of nucleosynthesis during the r-process [ 1 1 ], quantum electrodynamics in strong Coulomb field, and more. Much progress has been made experimentally, and 3473 nuclides had been discovered as of 2018 [ 2 2 ]. Over 8000 nuclides might be bounded in a nuclear chart predicted by nuclear density functional theory [ 3 3 ]. There are several ways to create new isotopes in different mass domains, i.e., projectile fragmentation, fission of transactinide nuclei, fusion–evaporation mechanism, multinucleon transfer (MNT) reactions, etc. [ 4 4 ]. Up to now, the heavy-mass area with neutron number N > 100 is still blank and far from the neutron-drip line [ 5 5 ]. In recent years, the MNT mechanism near Coulomb barrier energies has been studied both experimentally and through different theoretical models, which provides a possible way to reach the heavy neutron-rich nuclei region. Theorists have developed some theoretical models to study MNT reactions, such as the dinuclear system (DNS) model [ 6 6 - 9 9 6-9 ], the GRAZING model [ 10 10 ], and a dynamical model based on the Langevin equations [ 11 11 , 12 12 ]. Furthermore, microscopic approaches based on the nucleon degree of freedom, the time-dependent Ha rtree-Fock (TDHF) approach [ 13 13 - 17 17 13-17 ] and the improved quantum molecular dynamics (ImQMD) model [ 18 18 , 19 19 ] have also been developed for describing the transfer reactions in nuclear collisions. Compared with other models, the DNS model can better consider the structural effects through the potential energy surface (PES), and the calculation time is short. In experiments, attempts to produce neutron-rich nuclei near N = 126 have been made through reactions of 136 Xe + 208 Pb or 64 Ni + 208 Pb at the Flerov Laboratory of Nuclear Reactions in Dubna [ 20 20 ], the Argonne National Laboratory [ 21 21 ], and the Helmholtz Centre for Heavy Ion Research (GSI) [ 22 22 ]. The neutron-rich isotopes via the MNT reaction 136 Xe + 198 Pt were also investigated at Grand Accélérateur National d’Ions Lourds (GANIL) in France [ 23 23 ]. Nearly symmetric reactions of massive nuclei, such as 204 Hg+ 198 Pt [ 24 24 ], also have been used to produce new neutron-rich isotopes at Argonne National Laboratory. Recently, the product yields in the MNT reaction for the system 204 Hg+ 208 Pb at both incident energies of 4.79 and 5.6 MeV/nucleon have been to test the models [ 25 25 ]. Modifications of the theoretical models are still needed to understand the experimental data. The inverse quasi-fission process occurs in nearly symmetric nuclear collisions, which provides a way to explore the competition of shape evolution and nucleon transfer dissipation and is helpful for modifying the theoretical models. The advantage of the DNS model over other models is the inclusion of structure effects in the nucleon transfer process, i.e., the shell and odd-even effects, fission heavy fragments, etc. However, the shape evolution is not considered in the DNS model, in which the individual properties of two fragments at the touching configuration are implemented. Nearly symmetric nuclear collisions will enlarge the shape effect on the observables, such as the fragment yields, energy spectra, and angular distribution. There are still some open problems in terms of transfer reactions, i.e., the mechanism of preequilibrium cluster emission, the mass limit of new isotopes with stable heavy target nuclides, the stiffness of the nuclear surface during the nucleon transfer process, etc. In this work, the dynamical mechanism in MNT reactions of nearly symmetric systems was investigated with the DNS model. In Sect. 2 we give a brief description of the DNS model. Calculated results and discussion are presented in Sect. 3. A summary and perspective on the MNT reactions are presented in Sect. 4. 2 Brief description of the DNS model 1 The DNS concept was proposed by Volkov for describing the deep inelastic heavy-ion collisions near Coulomb barrier energies[ 26 26 ], where a few nucleon transfers were treated under the assumption of individual properties of DNS fragments. Application of the approach to superheavy nucleus formation via massive fusion reactions in competition with the quasifission process was used for the first time by Adamian et al. [ 27 27 , 28 28 ]. The Lanzhou Group performed modifications of the relative motion energy and angular momentum of two colliding nuclei coupling to nucleon transfer within the DNS concept [ 29 29 ], in which the nucleon transfer at the touching configuration of two colliding partners was coupled to the relative motion by solving a set of microscopically derived master equations. The contribution of quasi-fission and fission of heavy fragments was implemented in the nucleon transfer process. The dynamical mechanism in the fusion–evaporation and transfer reactions for producing the superheavy nuclei and exotic heavy nuclei was thoroughly investigated. Recently, the production cross sections of fragments in the MNT reactions around Coulomb barrier energies have been calculated using the DNS model [ 30 30 , 31 31 ], particularly neutron-rich isotopes around the shell closure N = 126 and 162. The production cross sections of the MNT fragments were evaluated by <math display="block" id="M1"><mtable columnalign="left"><mtr><mtd><msub><mi>σ</mi><mrow><mtext>tr</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><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><msub><mi>σ</mi><mrow><mtext>cap</mtext></mrow></msub><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><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mi>f</mi><mo stretchy="false">(</mo><mi>B</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>×</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><msub><mi>J</mi><mn>1</mn></msub><mn>,</mn><mi>B</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>×</mo><msub><mi>W</mi><mrow><mtext>sur</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mn>1</mn></msub><mn>,</mn><msub><mi>J</mi><mn>1</mn></msub><mn>,</mn><mi>s</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>B</mi><mn>.</mn></mtd></mtr></mtable></math> Here, E 1 and J 1 denote the excitation energy and angular momentum for the DNS fragment ( Z 1 , N 1 ), respectively, which are the origin from the center of mass energy E c.m. and the angular momentum J of the colliding system at the barrier B . The maximal angular momentum J max corresponds to the grazing nucleus–nucleus collisions. The emission of γ rays and light particles in competition with binary fission is implemented in the evaluation of survival probability W sur as a statistical model and s is the sum of evaporating particles, such as neutrons, protons, and α [ 32 32 ]. The capture cross section was calculated by the Hill-Wheeler formula with the barrier distribution approach [ 33 33 ]. The barrier distribution was taken as the Ga ussian form <math id="M2"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>B</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mi>exp</mi><mo stretchy="false">[</mo><mo>−</mo><msup><mrow><mo stretchy="false">(</mo><mo stretchy="false">(</mo><mi>B</mi><mo>−</mo><msub><mi>B</mi><mtext>m</mtext></msub><mo stretchy="false">)</mo><mo>/</mo><mi>Δ</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo stretchy="false">]</mo></mrow></math> , with the normalization constant satisfying the unity relation f ( B )d B =1. The quantities B m and Δ were evaluated by B m =( B C + B S )/2 and Δ = ( B C + B S )/2, respectively. B C and B S are the Coulomb barrier at waist-to-waist orientation and the minimum barrier from varying the quadrupole deformation parameters of the colliding partners. The nucleon transfer is coupled to the relative motion energy and rotation energy. The time evolution of the probability P ( Z 1 , N 1 , E 1 , t ) for fragment 1 with proton number Z 1 , neutron number N 1 , and local excitation energy E 1 is described by the following master equations: <math display="block" id="M3"><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mtext>d</mtext><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><mtext>d</mtext><mi>t</mi></mrow></mfrac></mtd></mtr><mtr><mtd><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></mtd></mtr><mtr><mtd><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></mtd></mtr><mtr><mtd><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>×</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></mtd></mtr><mtr><mtd><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><mn>.</mn></mtd></mtr><mtr><mtd><mrow/></mtd></mtr></mtable></math> The mean transition probability <math id="M4"><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> ( <math id="M5"><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><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msubsup><mi>N</mi><mn>1</mn><mo>'</mo></msubsup></mrow></msub></mrow></math> ) is calculated from the single-particle potential in valence space, and the microscopic dimension <math id="M6"><mrow><msub><mi>d</mi><mrow><msub><mi>Z</mi><mn>1</mn></msub><mn>,</mn><msub><mi>Z</mi><mn>1</mn></msub></mrow></msub></mrow></math> is evaluated from the valence states and valence nucleons. We only include sequential nucleon transfer with the master equations by the relation of <math id="M7"><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="M8"><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> . The nucleons around the Fermi surface are active for transferring DNS fragments. The local excitation energy E 1 is related to the dissipation energy of the relative motion and the potential energy surface (PES) of the DNS, which increases with relaxation time. The excitation energy at the equilibrium stage is distributed into the fragments according to mass. The angular momenta of the primary fragments are determined by the moment of inertia. The PES is evaluated by <math display="block" id="M9"><mtable columnalign="left"><mtr><mtd><mi>U</mi><mo stretchy="false">(</mo><mo>{</mo><mi>α</mi><mo>}</mo><mo stretchy="false">)</mo><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><mi>B</mi><mo stretchy="false">(</mo><mi>Z</mi><mn>,</mn><mi>N</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mi>V</mi><mo stretchy="false">(</mo><mo>{</mo><mi>α</mi><mo>}</mo><mo stretchy="false">)</mo><mn>,</mn></mtd></mtr></mtable></math> where Z and N are the proton and neutron numbers of the composite system, respectively; they satisfy the relations Z 1 + Z 2 = Z and N 1 + N 2 = N . The symbol α represents the quantities of Z 1 , N 1 , Z 2 , N 2 , J , R , β 1 , β 2 , θ 1 , and θ 2 . The binding energies B ( Z i , N i ) ( i =1,2) and B ( Z , N ) for the fragments and the compound nucleus (composite system without microscopic corrections for Z > 137) are calculated by the droplet model. The ground-state quadrupole deformations β i are considered in the calculation. The angles θ i are the collision orientations for the deformed nuclei. The interaction potential V ({ α }) of the two DNS fragments is composed of nucleus–nucleus potential, Coulomb potential, and centrifugal force at the touching distance [ 34 34 , 35 35 ]. Figure 1 Figure 1 shows the PES in the reaction 204 Hg+ 198 Pt at the tip–tip orientation. The incident position of projectile and target nuclides is indicated by a pentagram in the chart. The occupation probability is dissipated from the entrance point to possible states once overcoming the local excitation energy. The pocket in the PES appears around the shell closure. The positive local excitation energy enables nucleon diffusion in the valence space. The structure effects of the fragment production in the MNT reactions are embodied in the PES. The low valley in the neutron-rich region is available for new isotope formation because of the large local excitation energy. The emission angle of the MNT fragments is helpful for managing detectors in experiments, which is related to the reaction system and beam energy. We use a deflection function method to evaluate the fragment angle, where the fragment mass number, angular momentum, and local excitation influence the magnitude. The deflection angle is composed of the Coulomb and nuclear interaction as [ 36 36 ] <math display="block" id="M10"><mrow><mi>Θ</mi><mo stretchy="false">(</mo><msub><mi>l</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>=</mo><mi>Θ</mi><msub><mrow><mo stretchy="false">(</mo><msub><mi>l</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mtext>C</mtext></msub><mo>+</mo><mi>Θ</mi><msub><mrow><mo stretchy="false">(</mo><msub><mi>l</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mtext>N</mtext></msub><mn>.</mn></mrow></math> The Coulomb scattering angle is given by the Rutherford function. The nuclear deflection angle is evaluated by <math display="block" id="M11"><mrow><mi>Θ</mi><msub><mrow><mo stretchy="false">(</mo><msub><mi>l</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mtext>N</mtext></msub><mo>=</mo><mi>β</mi><mi>Θ</mi><msubsup><mrow><mo stretchy="false">(</mo><msub><mi>l</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mtext>C</mtext><mrow><mtext>gr</mtext></mrow></msubsup><mfrac><mrow><msub><mi>l</mi><mi>i</mi></msub></mrow><mrow><msub><mi>l</mi><mrow><mtext>gr</mtext></mrow></msub></mrow></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>δ</mi><mi>β</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>l</mi><mi>i</mi></msub><mo>/</mo><msub><mi>l</mi><mrow><mtext>gr</mtext></mrow></msub></mrow></msup><mn>,</mn></mrow></math> where <math id="M12"><mrow><msubsup><mi>Θ</mi><mtext>C</mtext><mrow><mtext>gr</mtext></mrow></msubsup></mrow></math> is the Coulomb scattering angle at the grazing angular momentum with l gr = 0.22 R int [ A red ( E c.m. - V ( R int ))] 1/2 . A red and V ( R int ) correspond to the reduced mass of DNS fragments and interaction potential at the distance R int from the entrance channel, respectively. δ and β are parameterized by fitting the deep inelastic scattering in massive collisions as <math display="block" id="M13"><mtable columnalign="left"><mtr><mtd><mi>β</mi><mo>=</mo><mn>75</mn><mi>f</mi><mo stretchy="false">(</mo><mi>η</mi><mo stretchy="false">)</mo><mo>+</mo><mn>15,</mn><mtext> </mtext><mtext> </mtext><mi>η</mi><mo><</mo><mn>375</mn></mtd></mtr><mtr><mtd><mn>36</mn><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mn>2.17</mn><mo>×</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>3</mn></mrow></msup><mi>η</mi><mo stretchy="false">)</mo><mn>,</mn><mtext> </mtext><mtext> </mtext><mi>η</mi><mo>≥</mo><mn>375</mn></mtd></mtr></mtable></math> and <math display="block" id="M14"><mtable columnalign="left"><mtr><mtd><mi>δ</mi><mo>=</mo><mn>0.07</mn><mi>f</mi><mo stretchy="false">(</mo><mi>η</mi><mo stretchy="false">)</mo><mo>+</mo><mn>0.11,</mn><mtext> </mtext><mtext> </mtext><mi>η</mi><mo><</mo><mn>375</mn></mtd></mtr><mtr><mtd><mn>0.117</mn><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mn>1.34</mn><mo>×</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>4</mn></mrow></msup><mi>η</mi><mo stretchy="false">)</mo><mn>,</mn><mtext> </mtext><mtext> </mtext><mi>η</mi><mo>≥</mo><mn>375</mn></mtd></mtr></mtable></math> with <math display="block" id="M15"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>η</mi><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">[</mo><mn>1</mn><mo>+</mo><mi>exp</mi><mfrac><mrow><mi>η</mi><mo>−</mo><mn>235</mn></mrow><mrow><mn>32</mn></mrow></mfrac><msup><mo stretchy="false">]</mo><mrow><mo>−</mo><mn>1</mn></mrow></msup><mn>,</mn></mrow></math> where <math id="M16"><mrow><mi>η</mi><mo>=</mo><mfrac><mrow><msub><mi>Z</mi><mn>1</mn></msub><msub><mi>Z</mi><mn>2</mn></msub><msup><mi>e</mi><mn>2</mn></msup></mrow><mi>υ</mi></mfrac></mrow></math> , and <math id="M17"><mrow><mi>υ</mi><mo>=</mo><msqrt><mrow><mfrac><mn>2</mn><mrow><msub><mi>A</mi><mrow><mtext>red</mtext></mrow></msub></mrow></mfrac><mo stretchy="false">(</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>−</mo><mi>V</mi><mo stretchy="false">(</mo><msub><mi>R</mi><mrow><mtext>int</mtext></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow></msqrt></mrow></math> . 3 The results and discussion 1 The primary fragments in MNT reactions must be de-excited to obtain the second fragments. This de-excitation process is described using a statistical model [ 32 32 ]. The cross sections that form fragments in this nuclear reaction can reflect the dynamics of a collision and the nuclear decay process of superheavy nuclei in excited states. Therefore, we can study these processes through cross sections. Figure 2 Figure 2 shows the cross sections of primary and secondary fragments formed by the MNT reactions of 204 Hg + 198 Pt and 198 Hg + 198 Pt. Due to the symmetry of the PES, the production cross sections of primary fragments (black lines) and second fragments (red lines)are nearly symmetric. The DNS dissipates into a broad mass region in the nuclear chart for primary fragments and approaches the neutron-rich domain. By comparing the distribution of the secondary fragment and the primary fragment formation cross sections with the mass number, it can be seen that the mass distribution of the secondary fragment formation cross section shifted to the left. This is because the primary fragment has a certain amount of excitation energy during the formation, and deactivation is reached. It should be noted that the production of heavy fragments ( m > 220) decreased drastically with the mass number. This is due to the strong Coulomb repulsion inside the composite system, meaning the composite system can exist only for a short time of several tens of 10 -22 s. The neck configuration is significant in symmetric damped collisions, which softens the potential barrier in the PES and is available for nucleon transfer [ 37 37 , 38 38 ]. The calculated results were compared with the data from Argonne National Laboratory for the nearly symmetric system 204 Hg+ 198 Pt, as shown in Fig. 3 Fig. 3 . The production cross sections of bismuth in the MNT reactions were obviously underestimated by approximately two orders of magnitude in the calculations. It should be noted that the total cross section of all MNT fragments in experiments was larger than the geometric cross section (4.5 barn) owing to the lack of normalization in the thick target. Isotopic cross sections of projectile-like fragments (PLF) and target-like fragments (TLF) were measured at a beam energy of E c.m. = 619 MeV close to the Coulomb barrier V C = 565 MeV. The production cross sections of the transfer fragments were underestimated by 20 times the experimental data. Overall, the isotopic structure is consistent with the measurements but with narrow distributions. Inclusion of the shape evolution in the dissipation process might enlarge the nucleon transfer and enhance the production of MNT fragments away from the PLF or TLF regions, which is missed in the DNS model. This work is in progress. The fragments formed in the symmetric or nearly symmetric transfer reactions undergo the inverse quasi-fission process. The nucleon transfer diffuses into enlarging the mass asymmetry under the PES. It provides the possibility for exploring the dissipation mechanism of two colliding heavy nuclei, stiffness of nuclear surfaces, superdeformation systems, new isotope production, etc. The transfer dynamics in the symmetric system were systematically investigated, as shown in Fig. 4 Fig. 4 for the reaction of 198 Pt + 198 Pt near barrier energies. Four different energies were selected for the calculation, i.e., 500, 542, 570, and 600 MeV, in the center of mass frame of colliding nuclei. A narrow diffusion on the isotopic distribution is obvious at the sub-barrier energy ( V C = 542 MeV). On the proton-rich side, the transfer cross section increases with the incident energy. The energy dependence becomes slight in the neutron-rich region because high incident energy enhances the local excitation energy, which is favorable for nucleon transfer. There, the energy is higher than the Coulomb barrier, and more nuclei can cross the barrier to participate in the reaction and increase the probability of the reaction, so the cross sections of primary fragments increase with the energy. Meanwhile, the high local excitation energy leads to reduction of the survival probability of primary fragments. Both the mass drift and decay process of excited fragments contribute to secondary fragment formation. The fragments produced around the shell closure survive because of the larger separation energy. The beam energy should be selected close to the Coulomb barrier of the reaction system to measure the MNT fragments in experiments. The anisotropy of fragment emission manifests the dynamical characteristics in MNT reactions, which is associated with deformation, dissipation of relative energy and angular momentum, configuration of nucleon transfer, etc. The shape of the angular distribution is an important criteria for the nuclear reaction mechanism. For example, the formation of a composite nucleus has a sufficiently long interaction time to allow the various degrees of freedom in the system to reach equilibrium, so the angular distribution of the reaction products appears to be 90 symmetrical in the center of mass. Estimation of the maximal yields from angular spectra of MNT fragments is beneficial for managing the detector system in experiments. Figure 5 Figure 5 shows the angular distributions of primary fragments of PLFs (black lines) and TLFs (red lines) in the damped collisions of 204 Hg + 198 Pt at the c.m. energies of 520, 580, 703 MeV. The PLFs tend to be emitted towards the forward direction with increasing beam energy. Opposite trends were found for the TLFs. A bump structure appears at high incident energy, which is around the Coulomb scattering angle. A broad distribution of MNT fragments is pronounced at the energy close to Coulomb barrier V C = 565 MeV. The maximal yields of PLFs and TLFs correspond to 55 and 30, respectively. The separation of the emission angle for PLFs and TLFs is enlarged with decreasing angular momentum of the colliding system. 4 Summary and conclusions 1 In summary, the dynamics in the MNT reactions around Coulomb barrier energies was investigated within the DNS model. The calculated isotopic distribution structure is consistent with the experimental data in the reaction of 204 Hg + 198 Pt. The maximal yields of fragments formed in the MNT reactions are located on the line of β stability. The isotopic cross sections are sensitive to the incident energy in the proton-rich region. The mass spectra of fragments in the symmetric system roughly exhibit a symmetric distribution, where the shell effect and odd–even effect contribute to the spectrum structure. The MNT fragments are produced in the forward-angle region along the colliding direction. The energy dependences of the angular distributions for PLFs and TLFs are opposite. The MNT fragments are distributed in the narrow mass region, which might be caused by neglect of shape evolution in the nuclear dissipation. The couplings of more degrees of freedom, i.e., neck degree of freedom, mass asymmetry, elongation evolution, etc., need to be implemented. Further modifications are in progress. References M. Terasawa , K. Sumiyoshi , T. Kajino et al ., New Nuclear Reaction Flow during r-Process Nucleosynthesis in Supernovae: Critical Role of Light Neutron-Rich Nuclei . Nucl. Phys. A 688 , 1 - 2 ( 2001 ). https://doi.org/10.1016/S0375-9474(01)00795-3 https://doi.org/10.1016/S0375-9474(01)00795-3 New Nuclear Reaction Flow during r-Process Nucleosynthesis in Supernovae: Critical Role of Light Neutron-Rich Nuclei M. Thoennessen , 2017 update of the discoveries of nuclides . Int. J. Mod. Phys. E 27 , 1830002 ( 2018 ). https://doi.org/10.1142/S0218301318300023 https://doi.org/10.1142/S0218301318300023 2017 update of the discoveries of nuclides J. Erler , N. Birge , M. Kortelainen et al ., The limits of the nuclear landscape . Nature (London) 486 , 509 ( 2012 ). https://doi.org/10.1038/nature11188 https://doi.org/10.1038/nature11188 The limits of the nuclear landscape F. Niu , P. H. Chen , Y. F. Guo et al ., Multinucleon transfer dynamics in heavy-ion collisions near Coulomb-barrier energies . Phys. Rev. C 96 , 064622 ( 2017 ). https://doi.org/10.1103/PhysRevC.96.064622 https://doi.org/10.1103/PhysRevC.96.064622 Multinucleon transfer dynamics in heavy-ion collisions near Coulomb-barrier energies T. Kurtukian-Nieto , J. Benlliure , K.-H. Schmidt et al ., Production cross sections of heavy neutron-rich nuclei approaching the nucleosynthesis r-process path around A = 195 . Phys. Rev. C 89 , 024616 ( 2014 ). https://doi.org/10.1103/PhysRevC.89.024616 https://doi.org/10.1103/PhysRevC.89.024616 Production cross sections of heavy neutron-rich nuclei approaching the nucleosynthesis r-process path around A = 195 Z. Q. Feng , G. M. Jin , and J. Q. Li , Production of new superheavy Z =108 - 114 nuclei with 238U, 244Pu, and 248,250Cm targets . Phys. Rev. C 80 , 057601 ( 2009 ). https://doi.org/10.1103/PhysRevC.80.057601 https://doi.org/10.1103/PhysRevC.80.057601 Production of new superheavy Z =108 - 114 nuclei with 238U, 244Pu, and 248,250Cm targets 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 ). https://doi.org/10.1103/PhysRevC.81.024604 https://doi.org/10.1103/PhysRevC.81.024604 ; Possibility of production of neutron-rich Zn and Ge isotopes in multinucleon transfer reactions at low energies G.G. Adamian , N.V. Antonenko , V.V. Sargsyan et al . , Predicted yields of new neutron-rich isotopes of nuclei with Z=64-80 in the multinucleon transfer reaction 48Ca+238U . Phys. Rev. C 81 , 057602 ( 2010 ). https://doi.org/10.1103/PhysRevC.81.057602 https://doi.org/10.1103/PhysRevC.81.057602 Predicted yields of new neutron-rich isotopes of nuclei with Z=64-80 in the multinucleon transfer reaction 48Ca+238U G. G. Adamian , N. V. Antonenko , and D. Lacroix , Production of neutron-rich Ca, Sn, and Xe isotopes in transfer-type reactions with radioactive beams . Phys. Rev. C 82 , 064611 ( 2010 ). https://doi.org/10.1103/PhysRevC.82.064611 https://doi.org/10.1103/PhysRevC.82.064611 Production of neutron-rich Ca, Sn, and Xe isotopes in transfer-type reactions with radioactive beams X. B. Yu , L. Zhu , Z. H. Wu et al , Predictions for production of superheavy nuclei with Z=105-112 in hot fusion reactions . Nucl. Sci. Tech 29 , 154 ( 2018 ). https://doi.org/10.1007/s41365-018-0501-2 https://doi.org/10.1007/s41365-018-0501-2 Predictions for production of superheavy nuclei with Z=105-112 in hot fusion reactions A. Winther , Grazing reactions in collisions between heavy nuclei . Nucl. Phys. A 572 , 191 ( 1994 ). https://doi.org/10.1016/0375-9474(94)90430-8 https://doi.org/10.1016/0375-9474(94)90430-8 ; Grazing reactions in collisions between heavy nuclei A. Winther , Dissipation, polarization and fluctuation in grazing heavy-ion collisions and the boundary to the chaotic regime . Nucl. Phys. A 594 , 203 ( 1995 ). https://doi.org/10.1016/0375-9474(95)00374-A https://doi.org/10.1016/0375-9474(95)00374-A Dissipation, polarization and fluctuation in grazing heavy-ion collisions and the boundary to the chaotic regime V. I. Zagrebaev and W. Greiner , Cross sections for the production of superheavy nuclei . Nucl. Phys. A 944 , 257 ( 2015 ). https://doi.org/10.1016/j.nuclphysa.2015.02.010 https://doi.org/10.1016/j.nuclphysa.2015.02.010 Cross sections for the production of superheavy nuclei D. Boilley , B. Cauchois , H. Lü et al , How accurately can we predict synthesis cross sections of superheavy elements? Nucl. Sci. Tech . 29 , 172 ( 2018 ). How accurately can we predict synthesis cross sections of superheavy elements? C. Golabek and C. Simenel , Collision Dynamics of Two U-238 Atomic Nuclei . Phys. Rev. Lett . 103 , 042701 ( 2009 ), https://10.1103/PhysRevLett.103.042701 https://10.1103/PhysRevLett.103.042701 ; Collision Dynamics of Two U-238 Atomic Nuclei C. Simenel , Particle-number fluctuations and correlations in transfer reactions obtained using the Balian-Vénéroni variational principle . Phys. Rev. Lett . 106 , 112502 ( 2011 ). https://doi.org/10.1103/PhysRevLett.106.112502 https://doi.org/10.1103/PhysRevLett.106.112502 Particle-number fluctuations and correlations in transfer reactions obtained using the Balian-Vénéroni variational principle K. Sekizawa and K. Yabana , Time-dependent Hartree-Fock calculations for multinucleon transfer processes in40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions . Phys. Rev. C 88 , 014614 ( 2013 ). https://doi.org/10.1103/PhysRevC.88.014614 https://doi.org/10.1103/PhysRevC.88.014614 ; Time-dependent Hartree-Fock calculations for multinucleon transfer processes in40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions K. Sekizawa and K. Yabana , Time-dependent Hartree-Fock calculations for multinucleon transfer processes in40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions . Phys. Rev. C 93 , 029902(E) ( 2016 ). https://doi.org/10.1103/PhysRevC.93.029902 https://doi.org/10.1103/PhysRevC.93.029902 Time-dependent Hartree-Fock calculations for multinucleon transfer processes in40,48Ca+124Sn, 40Ca+208Pb, and 58Ni+208Pb reactions L. Guo , C. Simenel , L. Shi , et al ., The role of tensor force in heavy-ion fusion dynamic . Phys. Lett. B , 782 , 401 - 405 ( 2018 ). https://doi.org/10.1016/j.physletb.2018.05.066 https://doi.org/10.1016/j.physletb.2018.05.066 The role of tensor force in heavy-ion fusion dynamic X. Jiang , N. Wang , Production mechanism of neutron-rich nuclei around N = 126 in the multi-nucleon transfer reaction 132Sn+ 208Pb . Chin. Phys. C 42 , 104105 ( 2018 ). https://doi.org/10.1088/1674-1137/42/10/104105 https://doi.org/10.1088/1674-1137/42/10/104105 Production mechanism of neutron-rich nuclei around N = 126 in the multi-nucleon transfer reaction 132Sn+ 208Pb K. Sekizawa , Microscopic description of production cross sections including deexcitation effects . Phys. Rev. C 96 , 014615 ( 2017 ). https://doi.org/10.1103/PhysRevC.96.014615 https://doi.org/10.1103/PhysRevC.96.014615 Microscopic description of production cross sections including deexcitation effects J. Tian , X. Wu , K. Zhao et al ., Properties of the composite systems formed in the reactions of 238U+238U and 232Th+250Cf . Phys. Rev. C 77 , 064603 ( 2008 ). https://doi.org/10.1103/PhysRevC.77.064603 https://doi.org/10.1103/PhysRevC.77.064603 Properties of the composite systems formed in the reactions of 238U+238U and 232Th+250Cf K. Zhao , Z. Li , X. Wu et al ., Production probability of superheavy fragments at various initial deformations and orientations in the 238U+238U reaction . Phys. Rev. C 88 , 044605 ( 2013 ), https://doi.org/10.1103/PhysRevC.88.044605 https://doi.org/10.1103/PhysRevC.88.044605 ; Production probability of superheavy fragments at various initial deformations and orientations in the 238U+238U reaction K. Zhao , Z. 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 ). https://doi.org/10.1103/PhysRevC.92.024613 https://doi.org/10.1103/PhysRevC.92.024613 Production mechanism of neutron-rich transuranium nuclei in 238U+238U collisions at near-barrier energies E. M. Kozulin , E. Vardaci , G. N. Knyazheva et al ., Mass distributions of the system 136Xe+208Pb at laboratory energies around the Coulomb barrier: A candidate reaction for the production of neutron-rich nuclei at N = 126 . Phys. Rev. C 86 , 044611 ( 2012 ). https://doi.org/10.1103/PhysRevC.86.044611 https://doi.org/10.1103/PhysRevC.86.044611 Mass distributions of the system 136Xe+208Pb at laboratory energies around the Coulomb barrier: A candidate reaction for the production of neutron-rich nuclei at N = 126 J. S. Barrett , W. Loveland , R. Yanez et al ., 136Xe+208 Pb reaction: A test of models of multinucleon transfer reactions . Phys. Rev. C 91 , 064615 ( 2015 ). https://doi.org/10.1103/PhysRevC.91.064615 https://doi.org/10.1103/PhysRevC.91.064615 136Xe+208 Pb reaction: A test of models of multinucleon transfer reactions O. Beliuskina , S. Heinz , V. Zagrebaev et al ., On the synthesis of neutron-rich isotopes along the N = 126 shell in multinucleon transfer reactions . Eur. Phys. J. A 50 , 161 ( 2014 ). https://doi.org/10.1140/epja/i2014-14161-3 https://doi.org/10.1140/epja/i2014-14161-3 On the synthesis of neutron-rich isotopes along the N = 126 shell in multinucleon transfer reactions 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 ). https://doi.org/10.1103/PhysRevLett.115.172503 https://doi.org/10.1103/PhysRevLett.115.172503 Pathway for the Production of Neutron-Rich Isotopes around the N = 126 Shell Closure T. Welsh , W. Loveland , R. Yanez et al , Modeling multi-nucleon transfer in symmetric collisions of massive nuclei . Phys. Lett. B 771 , 119 - 124 ( 2017 ). https://doi.org/10.1016/j.physletb.2017.05.044 https://doi.org/10.1016/j.physletb.2017.05.044 Modeling multi-nucleon transfer in symmetric collisions of massive nuclei V. V. Desai , A. Pica , W. Loveland et al , Multinucleon transfer in the interaction of 977 MeV and 1143 MeV 204Hg with 208Pb . Phys. Rev. C 101 , 034612 ( 2020 ). https://doi.org/10.1103/PhysRevC.101.034612 https://doi.org/10.1103/PhysRevC.101.034612 Multinucleon transfer in the interaction of 977 MeV and 1143 MeV 204Hg with 208Pb V. V. Volkov , Deep inelastic transfer reactions- The new type of reactions between complex nuclei , Phys. Rep . 44 , 93 ( 1978 ). https://doi.org/10.1016/0370-1573(78)90200-4 https://doi.org/10.1016/0370-1573(78)90200-4 Deep inelastic transfer reactions- The new type of reactions between complex nuclei G. G. Adamian , N. V. Antonenko , W. Scheid et al ., Treatment of competition between complete fusion and quasifission in collisions of heavy nuclei . Nucl. Phys. A 627 , 361 ( 1997 ). https://doi.org/10.1016/S0375-9474(97)00605-2 https://doi.org/10.1016/S0375-9474(97)00605-2 Treatment of competition between complete fusion and quasifission in collisions of heavy nuclei G. G. Adamian , N. V. Antonenko , W. Scheid et al ., Fusion cross sections for superheavy nuclei in the dinuclear system concept . Nucl. Phys. A 633 , 409 ( 1998 ). https://doi.org/10.1016/S0375-9474(98)00124-9 https://doi.org/10.1016/S0375-9474(98)00124-9 Fusion cross sections for superheavy nuclei in the dinuclear system concept Z. Q. Feng , G. M. Jin , F. Fu , and J. Q. Li , Production cross sections of superheavy nuclei based on dinuclear system model . Nucl. Phys. A 771 , 50 ( 2006 ); https://doi.org/10.1016/j.nuclphysa.2006.03.002 https://doi.org/10.1016/j.nuclphysa.2006.03.002 Production cross sections of superheavy nuclei based on dinuclear system model Z. Q. Feng , G. M. Jin , J. Q. Li , and W. Scheid , Formation of superheavy nuclei in cold fusion reactions . Phys. Rev. C 76 , 044606 ( 2007 ). https://doi.org/10.1103/PhysRevC.76.044606 https://doi.org/10.1103/PhysRevC.76.044606 Formation of superheavy nuclei in cold fusion reactions F. Niu , P. H. Chen , Y. F. Guo et al ., Effect of isospin diffusion on the production of neutron-rich nuclei in multinucleon transfer reactions . Phys. Rev. C 97 , 034609 ( 2018 ). https://doi.org/10.1103/PhysRevC.97.034609 https://doi.org/10.1103/PhysRevC.97.034609 Effect of isospin diffusion on the production of neutron-rich nuclei in multinucleon transfer reactions P. H. Chen , F. Niu , W, Zuo , and Z. Q. Feng , Approaching the neutron-rich heavy and superheavy nuclei by multinucleon transfer reactions with radioactive isotopes . Phys. Rev. C 101 , 024610 ( 2020 ). https://doi.org/10.1103/PhysRevC.101.024610 https://doi.org/10.1103/PhysRevC.101.024610 Approaching the neutron-rich heavy and superheavy nuclei by multinucleon transfer reactions with radioactive isotopes P. H. Chen , Z. Q. Feng , J. Q. Li , and H. F. Zhang , A statistical approach to describe highly excited heavy and superheavy nuclei . Chin. Phys. C 40 , 091002 ( 2016 ). https://iopscience.iop.org/1674-1137/40/9/091002 https://iopscience.iop.org/1674-1137/40/9/091002 A statistical approach to describe highly excited heavy and superheavy nuclei P. H. Chen , Z. Q. Feng , F. Niu et al ., Effect of isospin diffusion on the production of neutron-rich nuclei in multinucleon transfer reactions . Eur. Phys, J. A 53 , 95 ( 2017 ). https://doi.org/10.1103/PhysRevC.97.034609 https://doi.org/10.1103/PhysRevC.97.034609 Effect of isospin diffusion on the production of neutron-rich nuclei in multinucleon transfer reactions Z. Q. Feng , G. M. Jin , and J. Q. Li , Production of new superheavy Z=108-114 nuclei with 238U, 244Pu, and 248,250Cm targets . Phys. Rev. C 80 , 057601 ( 2009 ). https://doi.org/10.1103/PhysRevC.80.057601 https://doi.org/10.1103/PhysRevC.80.057601 Production of new superheavy Z=108-114 nuclei with 238U, 244Pu, and 248,250Cm targets Z. Q. Feng , G. M. Jin , and J. Q. Li , Influence of entrance channels on the formation of superheavy nuclei in massive fusion reactions . Nucl. Phys. A 836 , 82 ( 2010 ). https://doi.org/10.1016/j.nuclphysa.2010.01.244 https://doi.org/10.1016/j.nuclphysa.2010.01.244 Influence of entrance channels on the formation of superheavy nuclei in massive fusion reactions P. H. Chen , F. Niu , Y. F. Guo , et al ., Nuclear dynamics in multinucleon transfer reactions near Coulomb barrier energies . Nucl. Sci. Tech 29 , 185 ( 2018 ). https://doi.org/10.1007/s41365-018-0521-y https://doi.org/10.1007/s41365-018-0521-y Nuclear dynamics in multinucleon transfer reactions near Coulomb barrier energies H. Feldmeier , Transport phenomena in dissipative heavy-ion collisions: the one-body dissipation approach . Rep. Prog. Phys . 50 , 915 - 994 ( 1987 ). https://iopscience.iop.org/0034-4885/50/8/001 https://iopscience.iop.org/0034-4885/50/8/001 Transport phenomena in dissipative heavy-ion collisions: the one-body dissipation approach Z. Q. Feng , Nuclear dynamics and particle production near threshold energies in heavy-ion collisions . Nucl. Sci. Tech 29 , 40 ( 2018 ). https://doi.org/10.1007/s41365-018-0379-z https://doi.org/10.1007/s41365-018-0379-z Nuclear dynamics and particle production near threshold energies in heavy-ion collisions 10.1007/s41365-020-00770-1 This work was supported by the National Natural Science Foundation of China (Nos. 11722546 and 11675226) and the Talent Program of South China University of Technology. 2020-03-06 2020-04-14 2020-04-15 2020-05-27 2020-06-01 2021-04-10 © 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. 2020 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-020-00770-1 http://dx.doi.org/10.1007/s41365-020-00770-1 . 2020 Nuclear Science and Techniques 6 31