How to extend the chart of nuclides?

In the past 85 years the number of known nuclides increased by more than a factor of ten, resulting in 4000 presently known isotopes of 118 elements. This considerable progress we owe to the discovery of new reaction types along with the development of powerful accelerators and experimental techniques for separation and identification of reaction products. Model predictions indicate that still about 4000 further nuclides are waiting for their discovery. The vastest unexplored territory is located on the neutron-rich side in the upper half of the chart of nuclides and hides the answers to some of the most fundamental questions of nuclear physics like the limits of nuclear stability, element synthesis in the universe or stellar evolution. The access to these nuclei is presently limited by available beam intensities and/or the lack of appropriate methods for their production and identification. The latter concerns particularly new neutron-rich isotopes of transuranium and superheavy elements. To extend this area, the hope is presently based on multinucleon transfer reactions and on the application of fusion reactions with radioactive ion beams. But how promising are these approaches? Based on a survey of present-day knowledge, we will treat the questions where we currently are on our journey towards new territory on the chart of nuclides, how the chances are to gain new territory in the future and which challenges we will have to face.

illustrates how the number of known nuclides increased in the past 85 years. In 1935, Giorgio Fea published the first compilation of then known nuclides, arranged in a two-dimensional map according to their proton and neutron numbers [1]. It comprised 327 isotopes of elements from hydrogen to uranium, most of them stable or close to the stability line. The base for the enormous increase of new isotopes in the following 85 years was the discovery of new reaction types. But their effective application for isotope production was only enabled by the development of powerful accelerators for heavy ions and sensitive separation and detection techniques for the reaction products (a review on the discovery of nuclides, related reactions and techniques can be found in [2,3]). It was Ernest Rutherford who performed in 1919 the first man-made nuclear reaction in his lab in Manchester. By bombarding nitrogen with α particles from a radioactive source, Noteworthy are the two bulges on the neutron-rich side; the respective nuclei were produced as fission fragments. In the following 60 years, the number of known nuclides was more than doubled and the chart of 2015 contains 4000 different isotopes of 118 elements he triggered the fusion reaction 4 He + 14 N → 18 F * → 1 H + 17 O. After the mid-1920s, the first ion accelerators appeared which enabled nucleosynthesis in fusion reactions with protons and light ions. The discovery of nuclear fission in the late 1930s [4,5] entailed the observation of numerous neutron-rich fission fragments with mass numbers in a wide region around A = 100. The mechanism of spallation/fragmentation was discovered in the mid-1940s [6] and later turned out to be very efficient for the production of isotopes all-over the chart of nuclides up to uranium. Around the same time, the upcoming powerful new accelerator facilities enabled fusion reactions with heavy ions. Their application allowed one to enter the region of transuranium nuclei up to the heaviest known elements. Till today, fusion, fission and fragmentation contribute most effectively to the production of new isotopes. The chart in Fig. 2 gives an overview which nuclides can be populated with these reactions. But it reveals also the vast, still empty area between the presently known isotopes and the expected driplines. The largest unexplored territory is located in the upper half of the chart, on the neutron-rich side. The access to these nuclei is presently limited by available beam intensities for fragmentation and fission reactions and by the bending of the stability line toward the neutron axis for fusion reactions. The increase of beam intensities is a main goal in accelerator laboratories worldwide. But to reach new neutron-rich transuranium nuclei, also new ways for their production are required. Two main approaches are presently discussed: the application of multinucleon transfer (MNT) reactions and the application of fusion reactions with radioactive ion beams (RIBs).
The discovery of MNT reactions dates back to the late 1960s when the upcoming accelerators for heavy ions Fig. 2 Fragmentation, fusion and fission are presently the most effective reactions to produce (new) isotopes in the lab. In fragmentation reactions, neutron-rich as well as neutron-deficient nuclei up to uranium can be produced (yellow area). Fission leads to relatively neutron-rich intermediate heavy nuclei (green area). Fusion is presently the only method to synthesize transuranium and superheavy nuclei, but results in rather neutron-deficient reaction products due to the bending of the stability line (red area) allowed reactions to occur between complex projectile and target nuclei. Until then, the world of fusion reactions was simple. In collisions with very light projectiles, the orbital angular momentum J determines if fusion can take place or not. For angular momenta up to a critical value, J ≤ J crit , the projectile is absorbed by the target nucleus and a compound nucleus (CN) is formed. For J > J crit complete fusion is not possible anymore and only direct reactions or elastic scattering occurs. This picture was suddenly destroyed when scientists at the Joint Institute for Nuclear Research (JINR) in Dubna observed binary reaction products with broad mass, charge and angular distributions which emerged from deepinelastic low-energy heavy-ion collisions [7][8][9][10][11][12][13]. A bit later and independently, similar observations were also made in Orsay and Berkeley [14][15][16][17][18]. The observed characteristics suggested that deep-inelastic collisions (DIC) occupy an intermediate position between such different processes as direct reactions and CN reactions. There is no uniform terminology for this reaction type. Beside the denomination DIC, one can also find the expressions MNT, deep-inelastic transfer, quasi-fission, strongly damped collisions, relaxation phenomena and others. We will use the term MNT reactions.
From the very beginning it became obvious that MNT reactions are suitable to produce new isotopes; several neutron-rich nuclides of elements from carbon to thorium were discovered then. Also studies of MNT reactions with the goal to produce new transuranium and superheavy nuclei were performed in the 1970s up to the 1990s, but no new isotopes were observed in this region and finally, fusionevaporation reactions made the race at that time and the application of MNT reactions for isotope production abated for the moment. Ten years later the topic was raised again by the quest for suitable reaction types to access still unknown areas on the chart of nuclides. It started with new theoretical calculations and the particular interest to synthesize neutronrich medium-mass, heavy, and superheavy nuclei which are not reachable in fusion reactions with stable projectiles . Encouraged by cross-section calculations for mediummass and superheavy nuclei, MNT reactions were afterwards also suggested to synthesize neutron-rich nuclei along the N = 126 shell, which are of great interest for understanding the astrophysical r-process path. These nuclei can principally be produced in fragmentation and fission reactions, but the calculated MNT cross-sections are promising large and it is worth to look closer if they can fulfill their promise for large yields. A completely different approach is the application of fusion reactions with RIBs, which is of particular interest for the synthesis of neutron-rich superheavy nuclei. The attention to this method increases with the upcoming powerful RIB facilities but its success depends crucially on available RIB intensities.
Meanwhile, MNT reactions and fusion reactions with RIBs are topical subjects in various laboratories around the world where appropriate separation and detection techniques are being explored -a challenging task regarding the expected tiny production cross-sections which can easily reach the sub-nanobarn scale. We should not forget that the reaction mechanism alone does not decide the success of the method. Equally important is the sensitivity of available experimental techniques for separation and detection.
But why do we make so big efforts to find new nuclides? A very basic goal is the better understanding and description of the nuclear force. Heavy nuclei are complex manybody systems and the long-range residual force of the strong interaction which binds the nucleons together can so far not be described in the framework of quantum chromodynamics. Therefore, many phenomenological models have been developed to treat nuclei in different regions of the chart of nuclides. Nuclei close to the stability line can be reasonably well described but for nuclei far from stability the lack of understanding increases. This makes obvious why nuclei close to the neutron and proton driplines or superheavy nuclei are highly desired to give new input for the models. New phenomena are expected for very exotic nuclei like for example neutron skins, appearance of new magic numbers and disappearance of the known ones, and new spherical shell closures in the region of superheavy nuclei. We should not forget that most of the astrophysical r-process path, which is responsible for the production of heavy elements in stellar nucleosynthesis, proceeds through terra incognita. According to a report by the National Research Council of the US National Academy of Sciences, the origin of heavy elements from iron to uranium remains one of the 11 greatest unanswered questions of modern physics [56]. So, there are many good reasons for making big efforts to find new nuclides.

Features and model concept
The large variances of mass (charge), energy and angular distributions of binary reaction products from DIC lead to the model concept of a molecule-like dinuclear system (DNS) which is formed by projectile and target nucleus during their interaction and enables the strong exchange of mass (charge) and energy between the reaction partners. This picture was postulated already soon after the discovery of DIC [57][58][59][60][61][62][63][64][65] and is still valid today. The schemes in Fig. 3 from early publications [57,65] illustrate this process. The nuclear force leads to a sticking of projectile and target nucleus, forming the DNS. Sticking times are typically on the order (10 −21 -10 −20 ) s. Once the DNS is formed, the exchange of nucleons and energy between the reaction partners starts. In collisions with non-zero angular momentum, the DNS can also rotate about its center of gravity, where for long sticking times rota- 3 Model concept of nuclear molecule formation and orbiting. The upper sketch is from an early publication of Volkov [57] and illustrates the idea of nuclear molecule formation as first step of deep-inelastic collisions at low energy. During the interaction, the molecule-like DNS evolves by continuous redistribution of nucleons, excitation energy and angular momentum between the nuclei before it finally decays. The lower drawing from Ref. [65] shows the concept of nuclear orbiting in DNS with angular momentum, resulting in three possible scattering trajectories. In all cases, the reaction products are emitted to the same angle, but the underlying reactions are quite different: 1 -quasielastic scattering in a grazing collision, 2 -deep-inelastic collision with sticking and rotation of the dinuclear system by 90 • , 3 -deep-inelastic collision with full rotation of a long-living DNS tions of 360 • or more can take place before the system decays. The DNS rotation explains the broad angular distributions of binary reaction products which resemble the angular distributions of fission fragments. The rotation of the DNS leads also to deviations of the reaction products trajectories with respect to the Coulomb trajectories of elastic scattering. DNS formation is usually accompanied by strong dissipation of kinetic energy which is mainly transformed into internal excitation of the system. This leads to the observed small kinetic energies of the reaction products with respect to the beam energy.
During formation and evolution of the DNS, there is continuous redistribution of nucleons, excitation energy and angular momentum between the two nuclei. The equilibrium between proton and neutron numbers occurs very fast in collisions of light nuclei and continuously proceeds in reactions with heavy nuclei [66][67][68][69]. As it evolves, the DNS can pass with some probability through any macroscopic configuration allowed by the conservation laws for particle number, charge and total energy. In the utmost case, the system can fuse and ends up in a compound nucleus. The two possible pathways of DNS evolution are illustrated in Fig. 4. If the system overcomes the fusion barrier, it reaches complete statistical equilibrium with respect to all degrees of Fig. 4 The two possible evolution paths of nuclear molecules: once the molecule-like DNS is formed, it can undergo complete fusion, resulting in a CN. The de-excitation of the CN proceeds either by evaporating nucleons and γ rays, leaving an evaporation residue (ER) or by fission. Alternatively, the DNS can scission before reaching the CN state. Before scission, large amounts of nucleons can be exchanged between the reaction partners (multinucleon transfer MNT/quasi-fisson, QF). When the system breaks up, the still excited projectile-like and target-like primary transfer products (PTP) are emitted. Their de-excitation takes place via evaporation of nucleons and γ rays leading to the secondary transfer products (STP) or, particularly in the case of heavy nuclei, by fission freedom, forming a CN. The CN can de-excite either by emitting nucleons, nucleon clusters and/or gamma rays or by fission ("fusion-fission"). Alternatively, the DNS can decay before it reaches the CN state. The decay occurs either for static reasons (dominance of Coulomb repulsion between the interacting nuclei) or for dynamical reasons (coupling of modes of motion). Before decay, large numbers of nucleons can be exchanged between the two nuclei and lead to reaction products far from the original projectile and target nuclei (MNT or quasi-fission reactions). This allows us to hope that also new exotic isotopes are among these products. After scission, the projectile-like and target-like MNT products ("primary transfer products") still carry the excitation energy of the DNS, which was distributed among them during the sticking time. Like in the case of a CN, the de-excitation can take place by emission of nucleons and gamma rays leading to residual MNT products in their ground state ("secondary transfer products"). Alternatively, the excited primary MNT products can undergo fission for which the probability increases with increasing mass, spin, and excitation energy of the nuclei. According to this picture, complete-fusion and MNT reactions are closely related processes, which originate both from the formation of a dinuclear system.

Observation of long-living rotating nuclear molecules
A key process to validate the model concept of DIC is the experimental observation of long-living rotating DNS. In one of the earliest experiments, Ar + Th, sticking and orbiting of the nuclei was already indicated by the angular and energy distributions of binary reaction products [11]. Related to these observations, the possible existence of nuclear molecules was postulated then in Ref. [65]: "The idea of nuclear orbiting in collisions with angular momenta close to the critical value implies another interesting effect which eventually could be found experimentally. If some of the colliding systems survive more than half of a revolution, … such a phenomenon would strongly indicate the existence of nuclear molecules rotating in the dynamical force equilibrium".
In the following decades, nuclear orbiting was observed in many further experiments and various collision systems, but signatures from DNS which performed "more than half of a revolution" were only found recently in the velocity spectra of target-like MNT products from collisions of Ni+Pb(Z = 110) at the GSI velocity filter SHIP [70] (Fig. 5). The spectra reveal a double peak structure consisting of a pronounced maximum at velocities larger than 1.5 times the CN velocity v CN and a small peak at velocities v < 0.5v CN . For each observed nuclide the two peaks are located symmetrically with respect to the CN velocity. The peak positions and intensities shift systematically with the proton number of the MNT products. The most consistent explanation for the observed phenomenon is the formation of molecule-like rotating DNS which live sufficiently long to perform rotations by at least 180 • . The high-velocity peak then originates from MNT products which were created in central collisions, without rotation. After scission of the DNS the projectile-like nucleus is emitted in backward direction and the target-like nucleus with velocity v in beam direction, resulting in a peak at v = v + v CN in the lab system. If the DNS has angular momentum, it will rotate about its center of gravity. After half of a revolution, the DNS is again aligned parallel to the beam axis where, however, projectile and target nucleus have changed their positions like indicated by the small drawings on top of Fig. 5. If reseparation of the DNS occurs in this orientation, the target-like nucleus is emitted in forward direction with a reduced velocity v = −v + v CN . Reaction products from other rotation angles are not accepted by the velocity filter which covers only forward angles of (0 ± 2) • .
But the velocity spectra provide still much more information about DNS evolution. Since the peak positions are directly related to the kinetic energy of the measured targetlike MNT products, they reveal the total kinetic energy (TKE) in the exit channel if one assumes a binary reaction process.   64 Ni + 207 Pb at beam energies of 4.80, 5.00, 5.20, 5.40, 5.53 and 5.92 MeV/nucleon. The beam energies in the center-of-mass frame are given in the inset. The data points for a fixed Z but for different beam energies are plotted with an offset for better discrimination. The energies expected for asymmetric fission fragments of the respective compound nucleus 271  The TKE values obtained in this way are shown in Fig. 6 as a function of the proton number Z of the target-like MNT product for six different beam energies from 4.80 to 5.92 MeV/nucleon. In all reactions leading to nuclei with Z ≥ 84, the TKE values are independent of the beam energy and they are located at or below the expected energy of fission fragments ("Viola energy") [72] from the respective compound nucleus 271 Ds. This indicates that the DNS behaves already quite similar to a compound nuclear system. Namely, that full dissipation of kinetic energy took place and it "forgot" the entrance channel properties. In this case, the TKE val-ues, which decrease continuously with increasing number of transferred protons, are determined by the Coulomb barrier of the emitted MNT products at the scission point and reveal a strong deformation of the DNS before reseparation. From the internuclear distance which can be determined from TKE, quadrupole deformations of β 2 ≈ 0.4 are deduced [70]. Further, the intensity ratios of the two velocity peaks reveal the lifetimes (i.e. sticking times) of the DNS if one assumes that the DNS decay follows the decay law (details of the method are described in [70]). With the same "velocity filter method", weak signals from rotating DNS were also observed in the much heavier system Xe + Pb(Z = 136) [73]. This reveals that DNS formation takes still place in systems with proton numbers far beyond the region of known superheavy nuclei, where the fusion probability is already negligibly small. And even one of the heaviest systems, U + U(Z = 184), which was investigated in the late 1970s at GSI [74] and 30 years later at the GANIL VAMOS spectrometer [75], reveals signatures of noticeable time delays during the interaction, correlated with a large mass flow and strong energy dissipation. This behavior of the two heavier systems is particularly interesting if one considers them for synthesis of new exotic MNT products. We will come back to this later.
The direct comparison of these three systems [73] leads to astonishing results. Figure 7 shows the evolution of energy dissipation, interaction times and shapes during nuclear contact as a function of the net number of transferred nucleons d A. The similarity of the qualitative behavior is obvious and reveals two phases of the reaction. In each system, the first phase is characterized by strong dissipation of kinetic energy (total kinetic energy loss, TKEL), which is revealed by the steep slope of the TKEL /E c.m. curves in Fig. 7a. In this phase only a small number of nucleons, not more than 5% of the total A, is exchanged between the two nuclei. Also interaction times and deformations increase strongly in that part of the curves (Fig. 7b, c) indicating that the process is time consuming. In the second phase, the deep-inelastic regime is entered and the situation reverses. Energy dissipation and interaction times become more and more independent of the number of transferred nucleons. Also the shapes change more moderately. It seems that large numbers of nucleons now flow with small friction and small time consumption. The system is obviously starting to evolve freely in the mass asymmetry coordinate. Amazing is that the transition from phase 1 to phase 2 seems to occur after the transfer of about 5% of the total nucleons, independent of the collision system.

DNS evolution pathways
Other important evidence for the model concept of DIC is the experimental verification of the different pathways of DNS evolution as sketched in Fig. 4. This can be done by studying quasi-fission (QF) and fusion-fission (FF) reac- For better comparison of the systems, TKEL is divided by the respective center-of-mass energy E c.m. . The x-axis represents the net number of transferred nucleons d A divided by the total nucleon number A CN of the respective compound system; d A/A CN is given in percent tions, and there are indeed numerous experimental data  which reveal these pathways. The tiny cross-sections of superheavy fusion products [104,105] indicate already that the majority of DNS never ends up in a fusion-evaporation residue. Because the DNS formation cross-sections are large and reach several 100 mb up to barns [76], it is obvious that Fig. 8 Example spectra from Ref. [80] showing the mass (A) versus TKE of binary reaction products from collisions of S + U(Z = 108), Ca + U(Z = 112) and Ni + U(Z = 120) measured with the CORSET setup of JINR. The respective excitation energies of the composite systems are indicated on top of the figures. The red frames mark fission-like events which comprise contributions from QF and FF losses on the way to the fusion-evaporation residue must be huge. Presently, there are groups at JINR, at the Japan Atomic Energy Agency (JAEA) in Tokai, and at the Australian National University (ANU) in Canberra, who make extensive studies of QF and FF reactions. Their experimental setups are based on similar principles. The detection systems allow one to measure the time of flight (TOF) and energy E of the reactions products, from which their mass A and total kinetic energy TKE in the exit channel is obtained. In addition, the setups allow for the measurement of the angular distributions of the reaction products.
At JINR, QF and FF reactions in (super)heavy collision systems are studied for many years with the two-arm TOF-E detection system CORSET [77]. The about 0.5 m long detector arms can be positioned at different angles to measure also the angular distributions of the reaction products. The obtained mass resolution is ΔA/A ≈ 1.5% which corresponds to (3)(4) mass units for very heavy nuclei (also the JAEA and ANU setups allow for similar resolutions). Figure 8 shows example spectra of A−TKE distributions of binary reaction products measured with the CORSET setup in collisions of S + U(Z = 108), Ca + U(Z = 112) and Ni+U(Z = 120) [80]. A pronounced feature of both, QF and FF is the strong dissipation of kinetic energy which accompanies these reactions leading to low TKE values with respect to the incident energy. Their most pronounced difference is in the A distributions. The mass distributions of QF fragments are mainly asymmetric with maxima in the vicinity of the original projectile and target nuclei while FF leads to symmetric mass distributions which peak at ( A 1 + A 2 )/2. However, there is also an overlap of the QF and FF distributions in the sense that QF can contain symmetric components (from long-living DNS) and FF asymmetric ones, leading to uncertainties if one uses the measured distributions to quantify the individual contributions of QF and FF in a specific reaction.
By trend, a decrease of the FF component with increasing A (and Z , respectively) is observed, like in Fig. 8. This Mass-angle distributions of binary reaction products from collisions of U + Ni at a center-of-mass energy of 288 MeV measured at the Australian National University ANU [83]. The masses are given as mass ratio M R between the two reaction products at the scission point of the DNS. The rather narrow angular distribution and M R close to the values of the entrance channel nuclei suggest short DNS lifetimes and QF as origin of the reaction products is explainable with the increasing Coulomb repulsion which more and more prevents CN formation and drives the DNS to scission before full statistical equilibrium is reached. Further it is noticed that the FF branch increases with increasing beam energy. This can be explained with the increasing nuclear overlap which boosts CN formation. Collection of data for a variety of different collision systems and at different beam energies allows for the quantification of the individual contributions of QF and FF pathways (see e.g. [79,82]).
If one represents the A distributions as a function of scattering angle, like it is usually done by the ANU group, one obtains in addition the information on DNS lifetimes. Figure 9 is an example for such mass-angle distributions (MAD) measured with the ANU setup for collisions of Ni + U(Z = 120) [83]. Principally, the width of the angular distributions increases with lifetime of the DNS. The MAD of the quite heavy Ni + U system in Fig. 9 reveals that the angular distributions are relatively narrow and the masses of the reaction products mostly in the vicinity of the projectile and target masses. This indicates rather short DNS lifetimes and only small contributions from FF.

Theoretical approaches to deep-inelastic collisions
The dynamics of the DNS is described by two main degrees of freedom, namely, by (a) the relative motion of the fusing nuclei, and (b) the mass and charge transfer between the nuclei. Additional degrees of freedom are the deformations and orientations, the neck degree of freedom, internal excitations and others. The mass and charge transfer is usually described in coordinates of mass and charge asymmetry where A 1 , A 2 and Z 1 , Z 2 are the mass and charge numbers of the two clusters. The coordinates η A,Z vary in the interval (−1, 1), where η = ±1 means that either A 1 = 0 or A 2 = 0, which is interpreted as complete fusion. In the case η = 0 (A 1 = A 2 , Z 1 = Z 2 ), the DNS consists of two symmetric clusters. The potential energy determines the dynamics of the reaction. It can be calculated with the Strutinsky method where the potentials contain the liquid drop (LD) and shell effects, which are usually determined by the two-center shell model of Maruhn and Greiner [106]. One can distinguish the present theoretical approaches in phenomenological (microscopic and collective) and selfconsistent. Both describe MNT and fusion reactions fairly well, but the phenomenological models are presently more widely used. Their common ground is the formation of a DNS and its trapping in the minimum of the nucleus-nucleus potential as first step of the reaction. Resulting from this, the theoretical study of the DNS dynamics allows one to consider within a unique approach MNT and fusion reactions. Since the nuclei have internal structure, which is quite complicated at the considered excitation energies, some simplifications are required for choosing the most probable path of evolution in the configuration space, instead of considering all possible paths, and for introducing macroscopic characteristics (particle flux, friction, shape evolution, etc.). Dynamical and statistical characteristics are manifested to different degrees at different stages of the collision. This makes it possible to model separately each stage of the collision and also simplifies the theoretical treatment.
The basis for an analysis of the above described effects is the microscopic model. In the microscopic theory, the evolution of the DNS can be determined either by the timedependent Schrödinger equation for the wave function of the system or by the Liouville equation for the density matrix. One cannot solve these equations without making simplifying assumptions. Moreover, the exact solution of these equations would give more information than can be experimentally verified. The successes of relatively simple phenomenological and semi-microscopic models [107][108][109][110][111][112][113][114][115] indicate that not all degrees of freedom have a significant effect on the dynamics of the system. In these models, only certain collective (macroscopic) degrees of freedom are considered, which are chosen a priori, on basis of the requirements for the interpretation of experimental data. The number of collective degrees of freedom taken explicitly into account can be reduced by using experimentally established differences in their relaxation times.
The task of a microscopic theory is to find the effective Hamiltonian of the collective motion and its coupling with the internal degrees of freedom. This coupling is a source of irreversible processes in the DNS such as transfer of mass, charge, etc. The existing models differ in the internal excitations that they introduce explicitly or implicitly. They may be coherent and incoherent particle-hole excitations and transfers of nucleons or clusters of nucleons [116]. The first attempt was made in the model of dissipative diabatic dynamics [117,118]. This model is restricted solely by the form of the coherent motion of the nucleons, their diabatic motion in the changing two-center mean field. The process of excitation of the nuclei includes only particle-hole states whose diabatic levels cross during the motion of the nuclei. Because of the diabatic nature of the single-particle motion, the coupling between the collective and internal degrees of freedom is coherent and nonlinear, in contrast to the linear-response theory. In the dissipative limit, the friction in the dissipative diabatic dynamics model is equivalent to one-body friction ("wall" formula [119]).
There are two approaches, which start from opposite assumptions, to describing the establishment of equilibrium in the DNS. In the first approach, the exchange of nucleons is described as a stochastic process without inertial effects ("overdamped" regime) [66]. In the second approach, the charge asymmetry is regarded as a collective mode [64,120,121]. In the fragmentation theory [120,121], the evolution of the system along the mass asymmetry coordinate is determined by the time-dependent Schrödinger equation. Integrating this equation in the classical trajectory of relative motion of the colliding nuclei, one can determine the width of the mass (charge) distribution. However, the quantum fluctuations, which can be important, are ignored in this adiabatic treatment.
The transport (stochastic) models are widely used to describe DIC [122][123][124][125][126][127][128][129]. The exchange of mass and charge and the transfer of energy and angular momentum from the relative motion to internal degrees of freedom can be successfully interpreted as a diffusion process. By its nature, the interaction of nuclei is a non-equilibrium statistical process: the DNS is not only excited but also relaxes to complete statistical equilibrium. In comparison to the Schrödinger equation, the transport equations describe the irreversible evolution of a probability distribution and not an amplitude. All transport theories proceed from a division of the degrees of freedom of the system into a slow macroscopic degree and a rapidly relaxing internal degree. By using explicitly the statistical hypotheses for the interaction operator of the collective and the internal degrees of freedom, one can deduce the kinetic equations from the Liouville equation of motion and obtain the microscopic transport coefficients. In the macroscopic diffusion models, the transport coefficients are assumed to be proportional to the ratio of the phase spaces of the states of excited system. This means that the direction of nucleon transfer is determined by the balance between the total energies of different configurations of the DNS. However, the disregard of nuclear shell structure limits the description of experimental data. For example, the diffusion models which use only a potential-energy surface cannot explain the strong growth in the variance of the mass distributions along with a weak change of the mean masses of the reaction partners [130].
There are various transport models such as the randommatrix approximation [122,127], the model of one-body dissipation [128,129] and the linear-response model [126]. The main kinetic equation in these approaches is the master equation or Fokker-Planck equation for the distribution function of collective coordinates and conjugate momenta in the phase space of the collective degrees of freedom. In the microscopical models, the transport coefficients of the Fokker-Planck equation have been obtained. In these models, the internal system is described as a sum of independent internal subsystems of each of the nuclei. In Refs. [124,125], there are incoherent particle-hole excitations and exchange of nucleons between the nuclei due to non-diagonal matrix elements of the single-particle potential. The model of Ref. [128] considers only exchange of nucleons and is a microscopic analog of the classical picture of the exchange of particles through a window during the collision of nuclei [119]. The statistical hypothesis enters the model together with the assumption of rapid randomization of nucleon motions in each nucleus. The relative simplicity of the models and their success in describing the loss of kinetic energy ("window" formula) and the widths of the mass (charge) distributions of the reaction products are intriguing.
In the mentioned models, shell effects are not explicitly treated and the details of single-particle spectra of reaction partners are not considered. The matrix elements for nucleon transition are determined by averaging over all shell configurations. An influence of the structure of interacting nuclei is manifested only in a strong dependence of the transport coefficients on the density of the single-particle levels. In Ref. [125], the interacting nuclei are treated in the Fermi-gas approximation, and that means that the structure of nuclei is taken into account only in an averaged manner through the ground-state energy and the level density parameters. The approximation is not based on the perturbation theory, in contrast to the model of Ref. [128], in which the coupling between the collective and internal motions is described in the first order of perturbation theory ("weak-coupling" limit). The transport model of Ref. [128] includes (classically) dynamical effects, whereas the model [124,125] does not consider explicitly the relative motion of the nuclei.
In the approach based on the linear-response theory [126], the main assumption is that in each point of the classical trajectory the internal system is close to thermodynamical equilibrium. In this case one can calculate the deviation of the density matrix of the internal system from equilibrium with the perturbation theory. The linear-response theory is formu-lated in the quasi-adiabatic approximation, i.e. the model is valid only for small collective velocities. The excitation of the internal system (incoherent particle-hole excitations in a mean-field potential) is generated after each infinitesimally short time interval by a change in the mean field of the DNS. The relative motion of the nuclei, the mass (charge) asymmetry, and the shapes of the nuclei are taken into account in this model. The macroscopic and microscopic degrees of freedom are coupled by a friction tensor. The dissipation in the linear-response theory is a quantum-mechanical version of the classical one-body friction ("wall" formula).
The assumption of the local equilibrium in the DNS is a good approximation for the final stage of the reaction but is not correct for the description of the initial stage. At the initial stage there are strong coherent excitations, which decay into complex incoherent states during the time in which the local equilibrium is established. To take into account these effects, the modified Fokker-Planck equation was obtained in [131].
In heavy-ion reactions at energies of a few MeV/nucleon above the Coulomb barrier, the DNS can be characterized as a system of non-interacting particles in a common timedependent mean field, since the average velocity (Fermi velocity) of the nucleons is much larger than the relative velocity of the ions. For the description of the initial stage of the reaction, the most general microscopic approach is through the solution of the multi-particle Schrödinger equation in the time-dependent Hartree-Fock (TDHF) approximation [132][133][134]. The TDHF contains a description of all degrees of freedom that are considered in the transport models. However, there are certain problems to separate the macroscopic observables, owing to the nonlinearity of the dynamical equations. With the TDHF one can obtain a good description of fission and DIC, the deflection functions, the dissipation of energy and relative angular momentum, the centroids of the mass (charge) distributions, etc. However, the calculated variances of the mass (charge) distributions are several times smaller than the experimental ones. The TDHF is also unable to describe well the fusion cross-section. Another serious shortcoming of the TDHF in the theory of nuclear reactions was its inability to take into account fluctuations of the mean field. In spite of the recent progress, the TDHF method is still time consuming for consideration of collisions of deformed nuclei. There are some alternative approaches based on the TDHF, such as the adiabatic time-dependent Hartree-Fock approximation (ATDHF) and semi-microscopic hydrodynamical models [135].
The multistage processes of heavy-ion collisions can be consistently described in the framework of the coupledchannels method [136]. For simplicity, one usually models the coupling of the channels by means of friction forces [129] and considers the collective coordinate coupled with the internal ones in the quantum diffusion approach [137][138][139].
There exist models [140][141][142][143][144][145][146] in which one does not take into account all the effects (for example, microscopic selfconsistency of the density and nuclear potential) included in the TDHF method. In the framework of these models, one has attempted to describe in a unified manner the heavy-ion collisions in a wide range of energies, from the quasi-elastic to the deep-inelastic regime. The model [147][148][149][150][151][152][153][154][155][156][157][158][159] allows one to take explicitly into account the influence of the shell structure of nuclei on the collision dynamics. Constituent elements of the model are the realistic scheme of single-particle levels, nucleon separation energies, single-particle matrix elements of inelastic transitions of nucleons in the nuclei due to the effect of the changing field and the matrix elements of nucleon transfer from one nucleus to another. The singleparticle approximation is improved by the phenomenological inclusion of the residual interaction. The model also takes into account in a simpler manner the fluctuations associated with the distortion of the mean field in the framework of the random-phase approximation. The study of the evolution of single-particle degrees of freedom results in the transport coefficients that provides the microscopic basis for the phenomenological description of DIC with the master equation or Fokker-Planck equation [160]. Using the explicit equations for the single-particle operators, the transition probabilities between the macroscopic states of the system differing in charge (mass) asymmetry were obtained in [147][148][149][150][151][152][153][154][155][156][157][158][159].

Microscopical description of MNT
These models are based on the assumption that at low collision energies the DNS stage plays an important role. For the quantum-mechanical description of the internal degrees of freedom we use the single-particle approximation with realistic schemes of the single-particle levels for each nucleus. Qualitatively, the process is represented as follows: Two potential wells (Woods-Saxon potentials), whose centroids move on classical trajectories, perturb each other giving rise to the nucleon exchange and transitions. This mechanism is mainly responsible for the kinetic energy dissipation. Twonucleon collisions are indirectly taken into account through the Fermi (thermal) occupation numbers. The effects of residual interaction like excitation of high-and low-lying collective states are not explicitly taken into account, although the surface vibration modes contribute to dissipation. The adiabaticity of the relative motion with respect to these vibrations reduces this effect. It is convenient to represent the total Hamiltonian H of the DNS in the form The Hamiltonian of the relative motion is the sum of the kinetic energy operator and the interaction potential V (R) of the nuclei in the ground state: Here, R is the distance between the centers of mass of two nuclei, P is the conjugate momentum and μ is the reduced mass of the system. The last two terms in (2) describe the internal motion of the non-interacting nuclei and the coupling of the relative motion to the internal degrees of freedom, respectively. One can obtain from (2) the classical limit of the equations of motion for the macroscopic collective variables R and P: where the average t| · · · |t is over the internal variables at time t. The classical motion of the nuclei additionally depends on the time-dependent non-conservative coupling potential t|V int |t , for the calculation of which it is necessary to consider the equation of motion for the single-particle density matrix. The single-particle basis is constructed from the asymptotic wave vectors of the single-particle states of the noninteracting nuclei-the projectile "P", |P , and the target nucleus "T ", |T -in the form The orthogonality of this basis is satisfied up to the second order in P|T . It is convenient to take the single-particle Hamiltonian as where m is the nucleon mass, and A = A P + A T is the total mass of the system. The mean single-particle potentials of the projectile, U P , and of the target nucleus, U T , include both, nuclear and Coulomb fields. In the second quantization representation, this Hamiltonian is written in the form where, up to the second order in P|T , Here, ε P(T ) are the energies of the unperturbed single-particle states, which are characterized by the set of quantum numbers P ≡ (n P , j P , l P , The diagonal matrix elements P|U T |P and T |U P |T characterize the shifts of the energies of the single-particle levels due to the interaction of the nuclei. The non-diagonal matrix elements χ (T ) P P and χ (P) T T correspond to particle-hole transitions between the levels in one of the nuclei under the influence of the mean field of the partner nucleus. In turn, the g PT determine the transitions of nucleons between the nuclei because of the action of the DNS mean field. The contributions to the matrix elements from non-inertial effects are not taken into account, since they are small [143].
Information about the evolution of the system can be obtained by solving the equation of motion for the singleparticle density matrix. However, the steady solution of the equation for the occupation numbers does not provide the Fermi thermal occupation numbers in the asymptotic limit (t → ∞) at fixed R, because the residual interaction between the nucleons has not been taken into account. In the reaction, there are simultaneously the excitation of nuclei and their relaxation to statistical equilibrium. Allowance for them in the Hamiltonian (2) leads to the appearance of the well-known two-particle collision term [161,162] in the equations for the occupation numbers. Explicit consideration of the residual interaction requires extensive calculations. As shown in the time-dependent density-matrix theory [163], the two-body dissipation effects play an important role in fusion. For DIC, the fluctuations calculated from the Balian-Vénéroni variational principle [164] are much larger than standard TDHF results and closer to experimental fluctuations. One can take into account the twoparticle collision integral in linearized form (τ approximation) [165,166]: where τ is the parameter of relaxation time, and n eq (R(t)) is the local quasi-equilibrium distribution, i.e., the equilibrium distribution for a fixed value of R. The formal solution for n i (t) = t|a + i a i |t is as follows [167]: where andω ik = (ε i −ε k ]/h. It is convenient to solve Eq. (10) iteratively, dividing the time interval (t − t 0 ) into subintervals: t 0 , t 0 + Δt, t 0 + 2Δt, etc. Then at times Δt < τ, Eq. (10) is approximately written as wherē and The dynamical (n i (t)) and quasi-equilibrium (i.e. thermal, n eq i (R(t))) occupation numbers are calculated at each iterative step Δt. The parameters of the Fermi distribution for proton and neutron subsystems of each nucleus are determined by the conservation laws for energy and nucleon numbers of each type. The relative contributions of the single-and two-particle components are determined by the only parameter τ of our model. Knowing the values of the occupation numbers n i (t), we can readily determine at any time t the mean numbers of protons, Z P(T ) and neutrons, N P(T ) , and their variances σ 2 Z or σ 2 N , and also find the fragment excitation energies E * P(T ) (t): where ε F P (F T ) (R(t)) is the Fermi energy of the light (P) and heavy (T ) fragment. The upper indices, Z and N , in the sum indicate the summation over the proton and neutron singleparticle levels, respectively. At each iterative step, we separate the contribution to the fragment excitation energy from the particle-hole excitations and from the exchange of nucleons.
As originally assumed, the kinetic energy of the relative motion, which is transformed into energy of internal excitation, is distributed between the fragments proportionally to their masses (thermodynamic equilibrium). However, it was shown later in a number of experiments that this assumption is not always correct. In some reactions, the excitation energy is distributed almost equally between the products of binary reactions [168][169][170]. In other reactions, the distribution of the excitation energy is between an equal distribution and one proportional to the masses of nuclei. In the reactions, where one nucleus is magic or close to magic, most of the excitation energy is concentrated in the complementary product at relatively large losses of kinetic energy [171,172].
The calculation of the friction forces requires the explicit treatment of the coupling of the relative motion to the internal motion [140][141][142][143][144][145]173]. The models differ in the internal excitations considered. They may be collective surface vibrations, giant resonances, incoherent particle-hole excitations, or transitions of nucleons between the nuclei. It is clear that the structure of the excited states of the nuclei and the coupling strength of various excitation modes to the relative motion influences the energy distribution between the fragments.

The nucleus-nucleus potential
The description of the DIC dynamics depends on the potential taken between the nuclei. In the phenomenological models, we discriminate between adiabatic and diabatic potentials. Adiabatic potentials represent the minimum of energy for a set of given collective coordinates and have an outside barrier, but in general no further inside barrier versus smaller internuclear distances. They allow for the formation Fig. 10 The diabatic (solid curve), diabatic time-dependent (dashed curve) and adiabatic (dotted curve) nucleus-nucleus potentials for 110 Pd + 110 Pd as a function of the coordinate λ which measures the length of the system in units of the diameter of the spherical nucleus 220 U [175]. The heights of the fusion (B λ f us ) and quasi-fission (B λ q f ) barriers are marked of the compound system as a melting of two nuclei. Adiabatic potentials as a function of η are usually smallest for symmetric fragmentations (η = 0) and therefore obtain the largest fusion cross-sections around η = 0. Diabatic potentials include the effects of relative motion. At every avoided crossing of the single-particle energies in the adiabatic twocenter shell model the nucleons follow the so-called diabatic path, a process depending on the relative velocity of the nuclei and denoted as the Landau-Zener effect for which the Pauli principle between the nuclei is responsible. As an example, we show in Fig. 10 adiabatic and diabatic potentials for 110 Pd + 110 Pd collisions as a function of the length of the nucleus-nucleus system [174,175]. As seen in Fig. 10, the time-dependent diabatic potential almost coincides with the phenomenological diabatic potential near touching configuration (1.45 < λ < 1.7).
The phenomenological diabatic nucleus-nucleus interaction potential between two nuclei is written as which includes the Coulomb interaction, the nuclear interaction, and the centrifugal term. The centrifugal potential V rot in the formed DNS is usually expressed under the assumption that this system has the rigid body moment of inertia. In the calculation of the Coulomb and nuclear interactions, the deformations of the DNS nuclei are important. The Coulomb potential for two quadrupole deformed nuclei is analytically calculated following [176]. The nuclear part of the nucleus-nucleus potential can be obtained by using the double-folding formalism [177] with the density-dependent nucleon-nucleon interaction. The nucleus-nucleus potentials calculated with this V N have a potential pocket resulting from the attractive nuclear and repulsive Coulomb interactions. With increasing Coulomb repulsion in the DNS the depth of the pocket decreases. For the pole-pole orientation of interacting nuclei, the potential pocket is situated at the dis- 5 fm (β P and β T are the quadrupole deformation parameters of the nuclei with radii R P and R T ) and keeps the DNS nuclei in contact. The value of B qf (Z , N ) coincides with the depth of this pocket and decreases with increasing Z and J . For J ≤ 30, the dependence of B qf (Z , N ) is weak and can be disregarded. The decaying DNS with given Z and N has to overcome the potential barrier in R at R b = R m + 1.5 fm on the potential-energy surface. The radius of a nucleus is parameterized as The density-constrained along with the TDHF was used in [178][179][180][181] to obtain the nucleus-nucleus potentials where E DC is the density-constrained energy at the instantaneous separation R(t), while E P and E T are the binding energies of the two nuclei obtained with the same effective interaction. The concept of an adiabatic reference state for a given TDHF state is used in Eq. (20). The calculation of the excitation energy is achieved by dividing the TDHF motion into a collective and intrinsic part [182]. The major assumption is that the collective part is determined by the density ρ(r, t) and current j(r, t). Consequently, the excitation energy is written as where E TDHF is the total energy of the dynamical system and E coll = E kin (ρ(t), j(t)) + E DC (ρ(t)) represents the collective energy of the system. Here, E kin (ρ(t), j(t)) = m 2 drj 2 (t)/ρ(t) is the kinetic part which is asymptotically equivalent to the kinetic energy 1 2 μṘ 2 of the relative motion of two nuclei. The calculated capture cross-sections with this V (R) are found to be in a reasonable agreement with data and other model calculations.

Transport coefficients
Using (4), we obtain the integro-differential equation The second term on the right-hand side contains effective forces, of both, conservative and dissipative nature, due to the coupling of the relative motion with the internal motion. This term depends explicitly not only on the relative separation R(t) but also on the current time t, in contrast to the linearresponse theory [126] in which the asymptotic values are taken. In the adiabatic limit, when the characteristic time of the collective motion is much larger than the relaxation time, we can expand the function V ik (R(t )) in Eq. (22) near the point t = t: The first term of this expansion gives the dynamical correction to the conservative force. Here, One can similarly write down the correction to the conservative potential V (R(t)) due to the particle-hole excitations and the exchange of nucleons. The second term in Eq. (23) leads to an irreversible loss of kinetic energy of the relative motion, which is transformed into internal excitation of the nuclei. The coefficients of radial and tangential friction, respectively, are written in the form where and V L ik (R(t)) is the partial form factor for the given angularmomentum transfer L [157,158].
The third term in Eq. (23) results in the correction to the reduced mass μ(R(t)). The transport coefficients can be similarly calculated for other collective variables such as the mass (charge) asymmetry and the deformation of the nuclei of the DNS. Using the macroscopic reduction procedure of the meanfield theory proposed in Ref. [183], the friction coefficient associated with one-body energy dissipation was extracted in the entrance channel of heavy-ion fusion reactions [184]. The magnitude and form factor of the reduced friction coefficient have a universal property for various reactions investigated and are in a good agreement with those calculated with Eq. (28). The nucleus-nucleus potentials obtained with the density-constrained TDHF method exhibit energy dependence [183] and account for dynamical effects. The magnitude and form factor of the extracted friction coefficient depend on the beam energy as well. At low energies, nuclear dissipation is dominated by the one-body dissipation mechanism (nucleon exchange plus excitation of surface modes) based on the mean-field dynamics.
In the stochastic mean-field approach [185], the meanfield fluctuations are incorporated into the dynamical evolution by including the initial state fluctuations in a stochastic approximation [186][187][188][189]. The diffusion coefficient for nucleon exchange is extracted in a semiclassical approximation and has a form similar to that in the phenomenological nucleon exchange model. Using this similarity, one also infers an expression for the nucleon drift coefficient in the semiclassical approximation. While the diffusion coefficient is insensitive to the neck motion, the nucleon drift coefficient is affected by this motion.

Excitation energy distributions, centroids of charge and mass distributions
To demonstrate the possibilities of the model presented, we show the calculated distribution of excitation energy between the products of deep inelastic heavy-ion collisions in the 56 Fe(505MeV) + 165 Ho reaction [190][191][192]. The results for other reactions can be found in Refs. [167,193]. The theoretical and experimental results for the ratio R P = E * P /(E * P + E * T ) of excitation energy of the projectile-like nucleus to the total excitation energy of the DNS are shown in Fig. 11. We present the relative contributions of nucleon exchange and particle-hole excitations in the process of dissipation of the kinetic energy of relative motion. The calculated ratios R (ex) ) for particle- Fig. 11 a Ratio R P of the excitation energy (E * P ) of the projectile-like fragment to the total excitation energy E loss = E * P + E * T of the system as a function of E loss for the 56 Fe(505MeV) + 165 Ho reaction. The triangles represent experimental results [191]. The results of the model calculations [193] are shown by the solid curve. The dash-dotted curve marks the limit of thermal Dependencies of the calculated ratios R from nucleon exchange and from particle-hole excitations, respectively. The theoretical ratio of the excitation energy E ( ph) T due to nucleon exchange varies from 0.2 to 0.7 in all reactions considered [167,193]. As shown, in the reactions with 58 Ni, 74 Ge, and 56 Fe the nucleon exchange plays the main role, 3, in the dissipation of the kinetic energy. In the reactions with heavier nuclei, the importance of particle-hole mechanism increases, E The mean values of the charge (mass) distributions of the reaction products are also well described with the single-particle matrix elements (8) analytically determined as in Refs. [157,158,160,167]. The calculations showed that variation of the only free parameter τ from 5 × 10 −23 s to 25 × 10 −22 s does not affect significantly the results.
Similar conclusions on the dynamics of DIC have been drawn with other theoretical approaches. Applying the improved quantum molecular dynamics (ImQMD) model [194], the energy dependencies of the nucleus-nucleus potential, the friction parameter, and the random force char-acterizing a one-dimensional Langevin-type description of the heavy-ion fusion process were investigated [195,196]. As shown microscopically, the energy dissipation in the fusion process is governed by two mechanisms: One is caused by the nucleon exchanges between two fusing nuclei, and the other is due to a rearrangement of nucleons in the intrinsic system. The former mechanism monotonically increases the dissipative energy and shows a weak dependence on the incident energy, while the latter depends on both the relative distance between two fusing nuclei and the incident energy. As shown, the latter mechanism is responsible for the energy dependence of the fusion potential. The ImQMD model can be also used to predict the isotopic production cross-sections in multinucleon transfer reactions. While the GRAZING model [197][198][199][200] is suitable to estimate the production cross-sections as a result of the transfer of several nucleons, the approach based on the ImQMD model [40] provides wider isotopic distributions and is suitable for describing the transfer products rather far from the entrance channel.
The competition between collective motion and singleparticle degrees of freedom in dissipative processes was demonstrated in Ref. [201]. The spin-orbit force causes a significant enhancement of the dissipation in DIC. The theoretical fusion cross-section shows reasonably good agreement with the experimental data, considering that no free parameters are adjusted to reaction dynamics in the TDHF approach.

Master equations for nucleon transfer
To describe the mass (charge) distributions in DIC and QF reactions, and the fusion probability, one can apply master equations for the transfer of nucleons between the parts of the DNS [20,167,202]. As coordinates we choose the charge, neutron and mass numbers of both nuclei of the dinuclear system, namely Z P = Z , N P = N and A P = Z + N for the projectile-like (light) nucleus and Z T = Z tot − Z P , N T = N tot − N P and A T = A tot − A P for the target-like (heavy) nucleus. For the derivation of the master equations for the probability P Z ,N (t) to find the DNS in the state (Z , N ) at time t, we start with the Hamiltonian (2) of the DNS at R = R m and finally obtain [167] with the transition rates where the temperature-dependent Fermi occupation numbers for the single-particle states are used. The last line in (30) takes into account the DNS decay in R (Λ qf Z ,N ) and possible fission (Λ fis Z ,N ) of the heavy nucleus in the DNS. The initial condition for the master equations is The stability of the DNS with (Z , N ) is determined not only by the local minima of Δ (±,0) N ). So, the DNS evolution in mass asymmetry depends on the structure of the two nuclei. With increasing temperature of the system, the influence of the shell effects on the nucleon transfer decreases. This is reflected in a certain decrease of the oscillations of Δ when the temperature increases [167].
The transition rates Δ (±,0) Z ,N depend on the temperature only through the Fermi occupation numbers. In the case of strongly asymmetric DNS configurations, the separation between the single-particle levels in the light fragment is larger than the temperature T in the considered region of excitation energies, i.e., the deviation of occupation numbers from 1 or 0 is small, and the shell effects are manifested clearer. For nearly symmetric configurations, the change in the fragment Fermi energies Δε Z F P ,F T (Δε N F P ,F T ) on the transition Z → Z + 1 (N → N + 1) is less than T , the main contribution to the sums in (31) is made by the terms with Similar conclusions follow from consideration of the differences Δ Z ,N +1 . So, the weakening of the influence of shell structure on the nucleon transfer with increasing T occurs slower than the exponential decrease of the shell correction to the binding energy of the nuclei. One can approximately express the ratio, for example Δ Z ,N , through the ratio of corresponding level densities ρ Z ,N /ρ Z +1,N . Though this approximation firstly suggested in Ref. [123] and later used in Refs. [31,33,34,203] washes out the structure peculiarities in the dynamics, it considerably simplifies the calculations. Indeed, only the multidimensional potential-energy surface is required to calculate the transition rates. The master equations (30) can be spread to other collective variables besides Z and N . The evolution in the phase space can be also considered with the Langevin equations. Calculating many classical trajectories with the defined transition rates and random forces, the distribution in collective variables is obtained as in the solution of master equations.
The decay rates of the DNS in R are treated with the onedimensional Kramers rate They depend on the height B R qf of the outer potential barrier at the internuclear distance R b which is nearly independent of the angular momentum at J < 70 if the DNS has a large moment of inertia. The height of the barrier is The mass and charge yields are then obtained: The reaction time t 0 is determined by the balance equation of the probabilities: with the fusion probability P CN where Z BG and N BG are determined by the barrier for fusion in the mass (charge) asymmetry coordinates. The DNS with Z > Z BG evolves to the compound nucleus in a time less than 10 −20 s which is much shorter than the decay time of the compound nucleus. For small angular momenta related to the QF reactions, t 0 ≈ (3 − 4) × 10 −20 s, while in DIC being at large angular momenta the value of t 0 is about ten times smaller. So, the model presented suggest us one way to describe the yield of products of DIC and QF reactions. Indeed, the mechanisms of these processes are similar but differ in time scale. The calculated angular distributions of the cross-sections in Ref. [204], together with the observed experiment-theory agreement indicate that the production of exotic nuclei in DIC at near-barrier energies peaks at 0 degree. Their production is associated, in particular, with the removal of protons and transfer of neutrons. The differential reaction cross-sections seem to be maximal at 0 degrees for the most exotic fragments. The interaction time is considerably reduced in the reactions with two heavy nuclei where the value of B R qf is almost zero. Stable beams as 136 Xe or 192 Os and 208 Pb target were suggested in Ref. [31] for producing neutron-rich nuclei with neutron number close to 126 (important for astrophysics area of the last waiting point in the r -process nucleosynthesis) at energies about 20-30 MeV higher than the corresponding Coulomb barriers. In these reactions, the transfer of several protons from target to projectile has relatively high probability because of small Q value. However, the short interaction time due to the strong Coulomb repulsion could reduce the yield of transfer products. In the calculations, longer interaction time is achieved with larger friction coefficient corresponding to Γ > 3 MeV. However, the friction strength can not considerably exceed the estimate resulting from the TDHF calculations [183,183,184] and analysis of the experimental data on the kinetic energy dissipation. In Ref. [33], the neutron and proton transfer rates scaling factor is the crucial parameter. Indeed, its increase by a factor 2 leads to more than 10 times larger yield of neutron-rich products. The value of the nucleon transfer rate was taken to describe the known experimental data. As proposed, in the 238 U+ 248 Cm reaction at E c.m. = 770 MeV one can produce new neutron-rich isotopes of elements up to Sg with the cross-sections larger than 1 pb.
The mass and charge yields of products are given by The partial and total cross-section for quasi-fission can be calculated as: The fission probability of the heavier nucleus is defined as The capture cross-section, given as (the kinematic angular momentum), which is smaller than the critical angular momentum J crit . Trajectories with J ≥ J crit contribute to deep-inelastic and quasi-elastic collisions.
To explain the experimental TKE of the QF products, one has to regard the large polarizations of the DNS nuclei. For nearly symmetric dinuclear systems with (A P + A T )/2 − 20 ≤ A P ≤ (A P + A T )/2 + 20 we found deformations which are about 3-4 times larger than the deformations of the nuclei in their ground states. Let us assume the distribution of the fragments in charge, mass and deformation to be where we set the distributions of the deformations β P and β T as Gaussian distributions Here, σ 2 β = (hω vib /(2C vib )) coth(hω vib /(2kT )) with the frequency ω vib (Z , N ) and the stiffness parameter C vib (Z , N ) of the quadrupole vibrations are determined from experimental spectra. The average TKE is obtained as a function of the mass number A P of the light fragment As examples, we show here the charge distributions calculated in the reactions 238 U+ 40 Ca (340 MeV) and 238 U+ 48 Ca (425 MeV) (Fig. 12). The calculations were performed for interaction time t 0 = 2 × 10 −21 s and T = 1.5 MeV [205]. These parameters correspond to the collisions with the angular momenta close to J crit being characteristic of DIC. The  [167]. The solid lines show the experimental values of the cross-sections [64] theoretical results agree well with the experimental ones near the maxima of distributions. The underestimation of the tails of distributions is due to the contribution of collisions with J leading to longer interaction times.

Q gg systematics
As established experimentally [57], the cross-sections for production of the isotopes of nuclei with Z ≤ 8 in DIC satisfy the Q gg systematics. The apparent violation of the Q gg systematics at Z > 8 is due to the influence of the evaporation of nucleons and particles from the excited fragments. To justify the use of the Q gg systematics, statistical arguments, which are not included explicitly in the model, are usually employed. It is therefore not obvious in advance that P Z ,N , resulting from the master equations, follows the experimentally observed Q gg systematics. Calculations of P Z N (t 0 = 3 × 10 −21 s, T = 1.0 and 2.0 MeV) for the 197 Au + 40 Ar (292 MeV) reaction (Fig. 13) demonstrate that the Q gg systematics with allowance for corrections for nonpairing of the nucleons in their transfer from the donor to the acceptor nucleus is reproduced by the model. The dashed lines in Fig. 13 are fits to the experimental values.
In accordance with the Q gg systematics, the cross-sections for production of isotopes can be expressed as where ΔE C is the change of the Coulomb energy of the system due to proton transfer, and δ is the correction for the non- reaction [167]. The numbers next to the symbols are the mass numbers of isotopes. The dashed lines correspond to the approximation of the experimental data of Ref. [57] pairing of the nucleons transferred. The parameter T 0 differs from the thermodynamic temperature, i.e., the slope of the lines in the Q gg systematics characterizes T only through the functional dependence of T 0 and T . The slope of the lines in Fig. 13 decreases with increasing T , i.e., the dependence (46) is qualitatively reproduced. However, it changes by less than a factor 2, confirming the nonequivalence of T 0 and T . As found, the calculations, which use only the potentialenergy surface, do not always correctly describe the DNS evolution. For example, these models are unable to explain the large variance of the mass distribution and a small change of the mean value A in the 100 Mo + 40 Ar (270 MeV) reaction [89,130]. The reason is that the DNS evolution is affected not only by the relative positions of the Fermi surfaces of protons and neutrons in the interacting nuclei but also by the details of the shell structure (the number of free and occupied levels near the Fermi surface, the degree of their degeneracy, etc.). The first effect is taken into account by the potential surface. However, the others can be included in the consideration only microscopically. The model, which takes explicitly into account the shell effects for protons and neutrons, well describes the experimental results (Fig. 14).

Simplified statistical method
The production cross-section σ Z ,N of a primary light nucleus in transfer reactions is the product of the capture cross-section σ cap in the entrance reaction channel and the formation-decay probability Y Z ,N of the DNS configuration with charge and mass asymmetries given by Z and N : The primary light neutron-rich nuclei are mainly de-excited by neutron or gamma emissions. We consider only the reac-  [130]. The line shows the calculated results [167] tions leading to the excitation energies of light neutronrich nuclei smaller than their neutron separation energies S n (Z , N ). In this case the primary and secondary yields coincide.
If the projectiles and targets are deformed, the value of E min c.m. , at which the collisions of nuclei at all orientations become possible, is larger than the Coulomb barrier calculated for the spherical nuclei. In the collisions with smaller E c.m. the formation of the DNS is expected to be suppressed. Therefore, we treat E c.m. ≥ E min c.m. for which the capture cross-section is estimated with Eq. (40). Indeed, in the reactions of interest for production of exotic isotopes the values of E c.m. are always larger than E min c.m. to have enough energy for the formation of the DNS with very neutron-rich nucleus. The stability of the light neutron-rich nucleus is expected to be smaller in the excited rotational states than in the ground state. In order to produce the exotic nucleus with small angular momentum, only the partial waves with J ≤ 30 should be considered. Here, we assume that total angular momentum J is distributed in the DNS proportionally to the corresponding moments of inertia and set J cap = 30 in Eq. (40).
Based on the idea of Q gg systematics [57], instead of Eq. (34) one can use the simple statistical method to calculate Y Z N . This method for finding Y Z ,N uses the DNS potential energy calculated as in [167] where B L and B H are the mass excesses of the light and heavy fragments, respectively, which are taken from Ref. [206] for known nuclei and from Ref. [207] for unknown nuclei.
One can conclude from the calculations with Eq. (30) that the quasi-stationary regime is established quite fast in the DNS considered, specially along the trajectory in charge (mass) asymmetry corresponding to N /Z equilibrium in the DNS. For this trajectory, N = N 0 (Z ), i.e. the neutron number follows Z . Therefore, the formation probability for the configuration with Z and N 0 (Z ) is estimated as where C is a normalized constant. The temperature of the initial DNS and with the level-density parameter a = A tot /12 MeV −1 , where A tot is the total mass number of the system.
The formation of the DNS containing the light neutronrich nucleus with given Z is considered as a two-step process. The formation of the DNS with Z and N 0 is firstly treated. Then one should calculate the probability G Z ,N = Λ R Z ,N ,N 0 t 0 of the formation and decay of the DNS with exotic nucleus. Since the DNS with Z and N 0 is in the conditional minimum of the potential-energy surface, we use the Kramers-type expressions for the quasi-stationary which this DNS should overcome to observe the decay of the DNS with Z and N : where the preexponential factor depends on the friction and stiffness of the potential at the minimum and on the barrier. The temperature The main factor which restricts the time t 0 of the reaction and prohibits the formation of DNS containing the exotic nuclei is the evolution of the initial DNS to more symmetric configurations and decay of the DNS during this process. Therefore, the time of decay in R from the initial configuration or from more symmetric configurations mainly determines t 0 . We use again the Kramers-type expression for the quasi-stationary rate Λ R Z i ,N i of decay through the barrier and the rate Λ η sym Z i ,N i of symmetrization of the initial DNS through the barrier B η sym in the direction to more symmetric configurations: This expression leads to Y Z i ,N i ≈ 0.05 that is consistent with our previous calculations [20]. In the calculation of Y Z ,N the uncertainty related to the preexponential factor is estimated within the factor of 1.5. The suggested simplified approach is suitable if the initial DNS point in the reaction is located close to the N /Z equilibrium. If the injection point is considerably displaced from the N /Z equilibrium, the dynamical effects mainly contribute to the production of nuclei near the injection point and our statistical approach underestimates their yields.
The exotic nucleus as well as any nucleus far from the entrance channel of the reaction are the result of MNT between the projectile-like and target-like parts of the DNS. As a result, we can assume the thermal equilibrium in the DNS containing the exotic nucleus or in the DNS which is quite far from the initial DNS in the space (N , Z ). Indeed, for the formation of these DNS one needs quite a long time t 0 ≈ 10 −20 s at J ≤ 30. This allows us to assume the same temperature in the DNS nuclei and to define the excitation energy of the light nucleus with the mass A L as E * N )]A L /A tot . The deviation from the thermal equilibrium is expected only for the DNS decays near the injection point where the temperature of the heavy nucleus is smaller than the temperature of the light nucleus. Thus, assuming the thermal equilibrium in the DNS, we can underestimate the excitation of light primary nucleus and predict the upper limit of E * (Z i , N i ). Note that the partition of excitation energy in the DNS weakly influences Y Z ,N . Since in our calculations of the DNS potential energy the deformations of the nuclei [208] are close to their values in the ground states, the excitation energies of the DNS nuclei remain almost without changes after the DNS decays.
In order to test our method of calculation of σ Z ,N , we treated the production of Ti in the MNT reactions 58 Ni+ 208 Pb (E c.m. = 256.8 MeV) and 64 Ni + 238 U (E c.m. = 307.4 MeV) [209][210][211]. The excitation energies available in these reactions supply 2 neutron evaporation from the primary Ti isotopes having the maximal yields. In the 58 Ni + 208 Pb reaction 50 Ti and 52 Ti are produced with the cross-sections 1 and 0.2 mb [209,210], respectively, which are consistent with our calculated cross-sections 0.6 and 0.35 mb, respec- 15 The calculated (open circles) cross-sections of S isotopes are compared with the experimental ones (solid circles) for the 40 Ca+ 208 Pb reaction at E c.m. = 208.8 MeV [212] tively. In the 64 Ni + 238 U reaction the experimental [211] and theoretical cross-sections for 52 Ti are 0.5 and 1.6 mb, respectively. The description of the yields of S isotopes in the 40 Ca(249MeV) + 208 Pb reaction [212] is also a good verification of the model. Since the initial DNS in this reaction is not in the N /Z -equilibrium, we expect considerable contribution of the non-equilibrium processes to the yields of the products near the entrance channel; the yields of Ca, K, Ar, and Cl seem to be described better with the dynamical models [167,[198][199][200]. The yields of S isotopes, which are four protons from the entrance channel, are expected to be ruled by the potential-energy surface, and the simplified statistical method therefore seems to be suitable. Aside from the yield of 34 S, the calculated yields of S isotopes are in good agreement with the experimental data in Fig. 15. Therefore, the suggested method is suitable for prediction of the crosssections for the products of MNT reactions.
If the excitation energy of the primary products exceed the thresholds for neutron or charged particle emission, or we are interested in the production of neutron-deficient isotopes, the de-excitation of the primary products is taken into account. According to the model, there are nucleon drift and nucleon diffusion between the DNS nuclei which lead to the formation of excited CN and DNS configurations (DNS with different Z and A) with probabilities depending on the potential-energy surface and temperature of the system. The decay of excited CN and DNS configurations is described in a unique way, giving us an opportunity to calculate the production crosssections for residual nuclei in both, light particle and heavy cluster emission channels. Thus, light particle evaporation and complex fragment emissions are treated in the same way and the cluster emission is described under the assumption that the clusters are produced by the collective motion of the nuclear system in the charge asymmetry coordinate, with further thermal escape over the corresponding Coulomb barrier.
The cross-section of the residual nucleus with certain mass number A and charge number Z is given as where σ cap is the partial capture cross-section which defines the transition of the colliding nuclei over the Coulomb barrier and the formation of the initial DNS when the kinetic energy E c.m. and angular momentum J of the relative motion are transformed into the excitation energy and angular momentum of the DNS. The probability for the production of certain residual nucleus (Z , A) from the excited entrance channel DNS in the distinct decay channel is described by , one has to find the formation-emission probability W Z 1 ,A 1 (E c.m. , J ) of a certain light particle or cluster (Z 1 , A 1 ) from the excited system. Here, we consider the decay of the excited nuclear system as a sequential light particle (Z 1 < 2) evaporation, which includes neutrons, protons, deuterons, and tritons, and complex clusters (Z 1 ≥ 2). The maximum value of the angular momentum J max is limited either by the kinematic angular momentum J kin max or by the critical angular momentum J cr depending on which one is smaller: J max = min{J kin max , J cr }. The emission probability W Z 1 ,A 1 (E c.m. , J ) of a light particle or cluster (Z 1 , A 1 ) is calculated as the product of the CN or DNS formation probability and the CN or DNS decay probability: where the indices Z 1 and A 1 go over all possible channels from the neutron evaporation to the symmetric splitting. In the equilibrium limit of the master equation in the charge and mass asymmetries (see [213,214] for details) the probability P Z 1 ,A 1 is given by Here, n, p, d, and t-evaporation channels are taken into consideration with U (R m , Z 1 , A 1 , J ) = 0. The quasi-fission barrier B qf R , calculated as the difference between the bottom of the inner pocket and the top of external barrier, prevents the decay of the DNS in R with the factor P R Z 1 ,A 1 given as In Eqs. (55) and (56),  (Z 1 , A 1 , J ) is equal to the particle binding energy plus the value of the corresponding Coulomb barrier at Z 1 =0. The Fermi-gas model is employed to compute the temperature.
The neutron-deficient isotopes of nuclei are usually produced in very weak decay channels and their absolute production cross-sections are very sensitive to the excitation energy which is available for light particle or cluster emission. After emission of light particles or clusters, the excitation energy of the residual nucleus depends on emission barrier B qf R (Z 1 , A 1 , J ) and kinetic energy of the emitted particle or cluster. In Refs. [213,214], the kinetic energies of emitted particles have been assumed to be 2T Z ,A (J ) in the average. Here, we modified our calculations by taking into account the kinetic energy distributions of emitted light particles (n, p, d, and t) and clusters ( 3 He, 4 He, 12,14 C, …) as where is the kinetic energy of the light particle or cluster in the center-of-mass system. The actual value of in each decay event is chosen by the Monte Carlo method. This modification affects the cross-sections in very weak decay channels and has almost no influence on the major decay channels.
In the calculations, Eqs. (53) and (54) are used to treat the sequential statistical decay (the evaporation of light particles and/or the binary decay) of the hot CN until all fragments become cold (the excitation energy of fragments is smaller than its neutron emission threshold). The number n of the generated events in the Monte Carlo method was chosen according to the smallest decay probability, which is ∼ 1/n. Number n > 10 4 of iterations is large enough for obtaining the calculated results with quite a high accuracy.
To test the method suggested, we treated the production cross-sections of exotic nuclei in very weak decay channels. In Fig. 16, the calculated cross-sections are in a good agreement with the experimental data for the 58 Ni + 54 Fe reaction [215]. Note that our model reproduces the experimental data without adjusting the separation energies of light particles and clusters. So, the model can be employed for the predictions of excitation functions of exotic nuclei in the completefusion and QF reactions.
The results of the calculations [216] show that the maximum production cross-sections of 100 Sn are about 30 nb in the 58 Ni + 46 Ti reaction with stable beam and about 1 µb in the reactions with radioactive beams (Fig. 17). By taking into consideration the beam intensities and optimal production cross-sections, we found that the best candidate for the production of 100 Sn among reactions with radioactive ion beams is the 56 Ni + 46 Ti reaction via 2n emission channel. The production cross-sections of the doubly magic nucleus 100 Sn via light particle evaporation and cluster emission channels are comparable. The production cross-sections for the nuclei 101−103 Sn via xn, pxn, αxn and cluster decay channels are predicted as well in Ref. [216].

Early discoveries of new isotopes in MNT reactions
Shortly after their discovery, MNT reactions revealed already their potential for synthesis of new isotopes. About 75 new nuclides were discovered in the years 1969-1995 in different laboratories. All of them are on the neutron-rich side  [222][223][224][225][226][227][228][229]. Typical for these experiments was the application of relatively heavy and symmetric collision systems like Xe + Ta or W + W. The reaction products were separated according to Bρ with the GSI Online Mass Separator [230] and then stopped and accumulated in a catcher foil inside an ion source. From there, they were extracted as singly charged ions, re-accelerated to 60 keV and then implanted again. Using these samples with implanted reaction products, the isotopes were identified offline in β-, γ -, and X-ray decay spectroscopy. With this setup, also the then heaviest new MNT products 205 Au [231] and 232,234 Ac [232] were discovered, produced in collisions of 208 Pb and 238 U beams with Fig. 19 Dependence of A / Z (lower part) and σ 2 A /σ 2 Z (upper part) on the kinetic energy loss in the 40,48 Ca + 238 U reactions. The dots are the experimental data of Refs. [67,68]. The calculated results [205] are shown by lines W and Ta targets. The Bρ-ΔE-E and TOF-ΔE-E methods lead to mass (charge) resolutions of about 1%. Therefore it is only applicable for lighter nuclei, while the heavy MNT products were identified by decay tagging. This situation is still the same today.
To get an idea about the N /Z -equilibrium in DIC, we show in Fig. 19 exemplarily the results for 40 Ca + 238 U (340 MeV) and 48 Ca + 238 U (425 MeV) collisions. The experimental data [67,68] are compared with DNS model calculations [205]. The calculations were performed for the angular momenta close to the critical value J crit being characteristic of DIC. Good agreement between theoretical and experimental results is achieved when one considers the dependencies of the ratios σ 2 A /σ 2 Z and A / Z on energy loss E loss (Fig. 19). At small energy losses, σ 2 A /σ 2 Z is close to A / Z , and, as is well known [64], this means that proton exchange is dominant in the initial stage of the reaction. With increasing interaction time (increase of E loss ), the ratio σ 2 A /σ 2 Z becomes close to ( A / Z ) 2 , which corresponds to a strong correlation between the proton and neutron transfers. Analyzing the dependencies in Fig. 19, we may conclude that the establishment of N /Z equilibrium is a continuous process.

Early experiments
In parallel to the above described experiments, MNT reactions were explored already in their early days with the goal to synthesize (new) superheavy nuclei. Typical for these experiments was the application of heavy actinide targets from uranium (Z = 92) to einsteinium (Z = 99) in combination with various projectiles from O to U at energies up to ∼ 20% above the Coulomb barrier. Pioneering results were obtained by a collaboration of nuclear chemists in experiments at LBNL and at GSI, where they observed numer-ous MNT products with atomic numbers Z = (94 − 103) [86,[233][234][235][236][237][238]. Their results demonstrated that MNT reactions are suitable to populate also nuclei beyond uranium, although no new isotopes were observed then. The isotope identification was performed with radiochemical methods which one can roughly distinguish in on-line and off-line experiments. In on-line experiments, the nuclei were produced in a He gas filled target chamber and transported by gas-jet technique to chemical separation devices to separate the different chemical elements. After separation, samples were prepared for identification of the isotopes by decay tagging. The minimum time for one such cycle was about 60 s, which determines the minimum required half-lives of accessible nuclei. For long-living reaction products the easier off-line technique was used, where the reaction products were implanted in a catcher foil behind the target, followed by chemical element separation and decay tagging for isotope identification. The minimum time for a full cycle was about 30 min. The radiochemical methods were the most sensitive ones of these first generation experiments, leading to lower limit cross-sections of 20 nb. Figure 20 gives an overview on the most heavy transfer products observed in these experiments with 16,18 O, 20,22 Ne, 40,44,48 Ca, and 238 U projectiles in combination with 238 U, 248 Cm, and 254 Es targets [86,[233][234][235][236][237][238][239]. In all reactions, quite neutron-rich nuclei were populated.

New model predictions
After the mid 1990s, the interest in MNT reactions faded, but already a decade later their comeback started. It was triggered by the quest for appropriate methods to synthesize new neutron-rich superheavy nuclei, which cannot be reached in complete-fusion reactions with stable projectile and target nuclei. The DNS model [22] and the Langevin model [33,51,52] predicted cross-sections which gave cause for optimism and stimulated new experimental activities. However, these models result in two principally different points of view concerning the optimal projectile-target combinations. The DNS model, which uses diabatic internuclear potentials, suggests to use intermediate heavy neutron-rich projectiles, preferentially the neutron-rich 48 Ca, combined with the heaviest available actinide targets like Pu, Cm or Cf. It is obvious that such systems can principally not lead to neutron-rich isotopes of the heaviest known elements, but they provide sufficient neutrons to synthesize nuclei below Z = 110. Figure 21 shows the cross-sections, which are predicted by the DNS model for neutron-rich isotopes of Md to Hs (Z = 101 − 108) in MNT reactions of 48 Ca projectiles on 244,246,248 Cm targets [22]. Note that in these MNT reactions one can also produce new neutron-rich projectile-like fragments such as 84,86 Zn and 90,92 Ge with neutron numbers beyond N = 50 [29] or study the transfer-induced fission of isotopes from Md to Hs [240].
The Langevin model, in contrast, suggests to use the heaviest available collision systems like 238 U+ 238 U, 248 Cm, 254 Es [33,51,52]. Different than in the DNS model, the adiabatic potentials of the Langevin model result also in such giant  [22]. The beam energies were close to the respective Coulomb barrier. The denoted reaction products result after evaporation of one neutron from the excited primary transfer products systems in a sufficient contribution of collisions with long interaction times up to 5 × 10 −21 s [33], which is a necessary requirement to exchange large amounts of nucleons. The advantage of such heavy systems is that they provide large numbers of neutrons and can principally lead to neutron-rich nuclei of the heaviest known elements. But the cross-sections drop by about one order of magnitude with the transfer of every additional proton like revealed by the calculations in Fig. 22 [52][53][54]. According to Fig. 22, in the 238  Also the TDHF theory was adapted to calculate the MNT cross-sections for production of heavy nuclei [30,32,39,43,44,49,164,241]. The inverse QF mechanism, which is similar to that in Ref. [22], was proposed to produce neutron-rich transfermium nuclei [30]. The estimated interaction time is close to that obtained with the DNS model. The production of neutron-rich transfermium nuclei in the reactions 40,48 Ca, 58 Ni + 238 U, 248 Cm was also studied within the model [46] based on the DNS concept.
Extensive studies with 248 Cm targets were performed in the early experiments and also in a recent experiment by us [239]. They reveal the on first glance striking observation that nearly the same region of isotopes is populated with such different projectiles like O, Ne, Ca or U and that the cross-sections and positions of the isotopic distributions differ only little (Fig. 23). But this might be because the identi-fied nuclides are still located in the closer target environment and do therefore not yet reveal the diverging of the different collision systems. But despite still existing uncertainties, the different theoretical approaches and available experimental data leave a common picture. The calculated cross-sections for new neutron-rich superheavy MNT products are on the nanobarn scale or below. If we assume that cross-sections in the range 1 pb to 1 nb are accessible, MNT reactions can be an option to reach new isotopes with Z ≤ 108. For example, the so far not directly produced endpoint nuclei of the decay chains of superheavy isotopes from hot fusion reactions can be reached. The choice of the collision system seems not to be a very critical factor if the envisaged isotopes are not located too far from the entrance channel nuclei.

New experimental approaches and results
The theoretical and experimental state of the art seems encouraging and challenging at the same time. Encouraging, because the expected MNT cross-sections principally allow one to synthesize a good number of new transuranium isotopes. And challenging, because the cross-sections are at the present limits of experimental feasibility which requires the development of (new) separation and detection techniques as well as higher beam intensities. The present experimental activities are spread around the world, where Argonne, JINR, GANIL, GSI, Legnaro, MSU, RIKEN and others are participating. The development of separation and detection techniques is ongoing and also first experiments were performed.
New advances into the transuranium region were so far made at the velocity filter SHIP at GSI. The application of a velocity filter for separation of heavy MNT products is a new approach. Due to the acceptance angle of (0 ± 2) • , only nuclei which are emitted to forward angles can enter the velocity filter. Consequently, some of the detected nuclei originate from central collisions where projectile and target nucleus reach maximum possible overlap and sticking time at the given beam energy. This is the most favorable configuration for the exchange of large numbers of nucleons leading to the creation of isotopes far from the original projectile and target nuclei. Reaction products which are transmitted through SHIP are implanted in a position sensitive silicon strip detector in the focal plane, where they are identified by their α decays. In addition, also gamma decays are recorded with a germanium clover detector which is positioned behind the silicon detector. This enables also the identification of β emitters. This method turned out to be very sensitive leading to lower limit cross-sections of about 1 nb [242,243]. Further, the fast In-flight separation allowed one to detect also short-living nuclei with half-lives down to 20 µs, given by the conversion time plus dead time of the data acquisition system. Meanwhile, the study of MNT reactions in forward   48 Ca + 248 Cm at the velocity filter SHIP [239] direction using In-flight separation, is also tackled by other groups. So far, investigations were also performed at the gasfilled separator TASCA at GSI where heavy MNT products with proton numbers Z = (83 − 90) were identified in collisions of 50 Ti+ 249 Cf [244]. And measurements at zero degree of light neutron-rich projectile-like nuclei from 18 O + 238 U reactions were made at GANIL with the LISE achromatic spectrometer [204].
At SHIP, we investigated reactions of 48 Ca + 248 Cm at an energy slightly above the Coulomb barrier and identified more than 100 different target-like MNT products (Fig. 24). The heaviest ones are two isotopes of Fm and, most likely, the isotope 260 No, the latter with a total cross-section of 1.7 nb, corresponding to three events observed in two days of beamtime. Usually, a contiguous region of isotopes around the target and projectile nucleus is populated in MNT reactions. Therefore, one can assume that also many of the nuclei which are represented by blank squares in Fig. 24 were produced  [242]. Beside 248 Cm with a contribution of 96.85% (black square), the target contained also the isotopes 244−247 Cm with contributions of 0.0007%, 0.031%, 3.10% and 0.015% (black triangles). Five new neutron-deficient isotopes were also identified in this experiment and are marked by dots. One of the observed spontaneous fission activities was attributed tentatively to 260 No, which is indicated by the question mark in the respective square but could not be identified by decay tagging due to their too short or long half-lives and/or unfavorable decay channels like fission or beta decay.
Among the observed MNT products were also five new isotopes with Z ≥ 92. With this, the experiments at SHIP lead for the first time to the discovery of new transuranium isotopes in MNT reactions. Remarkable is that all new isotopes are located on the very neutron-deficient side of the chart. This results from the applied method for isotope identification via α decay tagging, which is so far the only possibility to identify such heavy and slow nuclei. However, the method works only if the nuclei have sufficiently short half-lives and lead to sufficiently long α decay chains. This requirement is mainly fulfilled by neutron-deficient transuranium and superheavy nuclei but largely fails on the neutronrich side of the chart.
Neutron-deficient nuclei around uranium are so far synthesized in fusion-evaporation reactions [245][246][247][248][249][250][251][252]. Therefore, it is interesting to have a closer look and see if MNT reactions might be after all the better alternative to populate the still quite empty neutron-deficient transuranium region. Figure 25 compares measured fusion residue cross-sections of various uranium isotopes [253][254][255][256][257][258][259] with the MNT cross-sections measured at SHIP [243]. Figure 25 can be regarded as a representative example because the behavior is very similar for other elements in the uranium region. If cross-sections from Fig. 25 The product of cross-section and experimental efficiency, σ ε, for uranium isotopes measured in fusion-evaporation reactions [253][254][255][256][257][258][259] (crosses) and in MNT reactions in 48 Ca + 248 Cm [243] (open circles). For 219 U a small offset on the A value is put to avoid an overlap of the data points from complete fusion and MNT more than one experiment were available for the same isotope, we put the largest known value. Instead of the pure production cross-sections, we put the product of cross-section and experimental efficiency, σ ε, because this is the actually relevant parameter which reflects the event count rate (assuming similar beam intensities and target thicknesses). In our example, ε is mainly determined by the angular acceptance of the applied separators. Figure 25 indicates that there is no overall valid answer if MNT or complete-fusion reactions are the better option.
Rather one must distinguish between "neutron-deficient" and "very neutron-deficient" isotopes. For nuclei well above A = 220 the value of σ ε is significantly larger in fusionevaporation reactions. This is, on one hand, caused by the larger angular acceptances of the separators for fusionevaporation residues. Another reason, on the second hand, is that these isotopes can be produced with very asymmetric projectile-target combinations like 22 Ne + 208 Pb → 230 U * where the low entrance channel Coulomb barriers and favorable Q values result in relatively large fusion residue crosssections. In this mass region, fusion-evaporation reactions appear superior to MNT. The situation changes for very neutron-deficient isotopes where the values of σ ε for MNT and fusion tend to approximate each other, despite the small angular efficiency for MNT products. To reach these nuclei in fusion reactions, more symmetric projectile-target combinations like 40 Ar+ 182 W → 222 U * must be applied. The larger entrance channel Coulomb barriers and less favorable Q values lead to relatively small fusion residue cross-sections. Here, a clear advantage of MNT reactions comes now into play. It arises from the broad excitation functions of MNT products which leads to a wide-band population of many different nuclides with sizeable yields in the same experiment, while fusion reactions are only selective on very few specific isotopes. In this region MNT seems indeed in competition with fusion reactions and might become an attractive option for future experiments.

Model predictions
This region is of great interest to understand the astrophysical r-process path which proceeds on the neutron-rich side of the nuclide chart and starts to traverse very neutron-rich, still unknown territory above the Sn region. So far, r-process nuclei are produced in the lab in fragmentation and fission reactions. But in the course of new model calculations, MNT reactions are suggested as an alternative because they reveal cross-sections which are comparable, or even larger than those of fragmentation reactions. The DNS as well as the Langevin model suggest to use targets close to this region like 208 Pb or 198 Pt, but differ in the choice of the projectile. The Langevin model obtains the largest cross-sections with heavy beams like 136 Xe, which have similar A/Z as the target nuclei, while the DNS model pleads for neutronrich but lighter projectiles such as 48 Ca or 64 Ni. As discussed earlier, this reflects the stronger dominance of the Coulomb force over the nuclear force in the DNS approach making very heavy systems unfavorable. Figures 26 and 27 show respective calculations from the DNS [28,35] and the Langevin model [26,41]. In the 136 Xe + 208 Pb reaction the Langevin model predicts larger cross-sections for production  [35] for the neutron-rich platinum isotopes 202 Pt (black symbols), 204 Pt (red symbols), 206 Pt (blue symbols) produced in collisions of the indicated projectiles with 198 Pt targets. The reaction products result from 0n evaporation channels of neutron-rich Pt isotopes than the DNS model predicts for the reactions indicated in Fig. 26. However, the theoretical cross-sections in Ref. [41] are in better agreement with the data on the neutron-deficient side of the mass distribution for each isotope. Therefore, the cross-sections for neutron-rich isotopes might be overestimated in Fig. 27 [260].
Interesting is another approach which is suggested by the DNS model [28]. MNT reactions in 48 Ca + 238 U are suggested to populate neutron-rich nuclei along N = 126. The idea behind it is that the nucleon flow in MNT reactions is oriented along the stability line. This means that with U targets, one can benefit from the preferred south-west flow of nucleons from target to projectile to populate isotopes along N = 126. Indeed, the production of primary neutronrich products occurs with large cross-sections. However, they have enough excitation energy to evaporate at least three neutrons. The neutron evaporation shifts the isotopic distributions to the left and reduces the yields of neutron-rich nuclei of interest. With Pt targets the population of these nuclei depends on the rather suppressed flow in south-east direction, but Pt isotopes are produced "cold" which provides a gain for the yields of neutron-rich isotopes. Figure 28 shows the DNS model calculations for Pt nuclei produced in the MNT reaction 48 Ca + 238 U. The primary isotopic distributions are now shifted toward the neutron-rich side. After neutron evaporation, for example, 204 Pt is produced with about 50 times smaller cross-section (0.4 nb) in the Ca + U reaction than in the Ca + Pt reaction (20 nb). It is up to future experiments to verify this result of the model.
Cross-section calculations for nuclei along N = 126 were also performed with other models [45,47,50] [28]. The yields of primary and secondary isotopes are shown by closed and open squares, respectively. Neutron evaporation channels for neutron-rich primary isotopes are indicated. The heaviest known isotopes are marked by arrows tems like W+U or Xe+Pt by using the model which is based on the DNS approach and the isospin-dependent quantum molecular dynamics model [50]. The results of the various models are quite different and can vary by several orders of magnitude for the same isotope [50,261]. The values in Table 1 give an impression. In Ref. [41], the nucleon transfer before reaching thermal equilibrium was taken into account. This results in larger cross-sections than in Ref. [26]. Also the nucleon transfer rates and relevant degrees of freedom were specified in a more refined way. Though the existing models describe well the maxima of isotopic distributions, there is a large uncertainty in the slope on the right-hand side of the isotopic distributions. Therefore, the experiments on the production of neutron-rich isotopes with N > 126 are desirable in order to reduce the theoretical uncertainties and define the most efficient reactions for production of new neutron-rich nuclei with N > 126.

Experimental state of the art
The population of neutron-rich nuclei along N = 126 in MNT reactions was already studied in different experimental approaches. One can roughly distinguish them in experiments which identified the heavy target-like products directly [260,[262][263][264] and in experiments which detected the projectile-like products and deduced from them the isotopic distributions of the complementary target-like reaction partners [265]. The results of the experiments which performed direct isotope identification reveal an astoundingly uniform picture, despite their different experimental approaches and collision systems. Krolas et al. started their first experiments about 15 years ago applying collisions of 58,64 Ni + 208 Pb at Legnaro [262,263]. In these "thick-target experiments" the reaction products were stopped in the target and their yield distributions were established from in- beam and off-line γ spectroscopy. Actually, these experiments aimed for spectroscopic studies of the MNT products, which are usually created in excited energy and spin states, by applying the γ -γ coincidence technique. Due to the stopping of the produced nuclei in the thick target foils, the emitted γ rays are not affected by Doppler broadening, which enables high-resolution spectroscopy. More recent studies were performed at Argonne National Laboratory (ANL) and by ourselves at GSI Helmholtz Centre. The ANL scientists investigated the 136 Xe + 208 Pb reaction in thick-target experiments and identified a broad region of projectile-like and target-like MNT products from Z = 48 to Z = 88 in off-line γ ray spectroscopy [260]. At GSI, we studied the reaction 64 Ni + 207 Pb [264] with the velocity filter SHIP to separate the target-like MNT products which were emitted to zero degrees. Isotope identification was performed in the focal plane of the velocity filter by γ spectroscopy. We want to stress that all experimental approaches used gamma spectroscopy to identify the target-like nuclei because it is presently the only method to identify heavy, low-energetic beta emitters. It is also noteworthy, that no new isotopes were observed in any of these experiments and that the lowest limit cross-sections were between (1-10) µb which is definitely not sensitive enough to advance deeper toward new neutron-rich nuclei. Figure 29 gives an overview on the isotopes with Z = 70−82 which were so far produced in MNT reactions and directly identified. One can notice that the observed nuclei are still many neutrons far from the most neutron-rich known isotopes of the respective elements. The so far heaviest isotope, 202 Pt, which was observed in MNT reactions was measured in Ref. [260].  136 Xe beams are clearly shifted toward the neutron-rich side. Since isotopes in this area are usually produced in fragmentation reactions, we show in Fig. 30 also measured fragmentation cross-sections [266][267][268][269][270]. Interesting is that the MNT cross-sections are indeed similar or even larger than fragmentation cross-sections, like predicted by the models. And the discrepancy between MNT and fragmentation appears to increase toward the neutron-rich side to the advantage of MNT reactions. Also, a trend is indicated that MNT cross-sections increase with respect to fragmentation cross-sections toward smaller proton numbers. This trend was particularly observed by Watanabe et al. who investigated projectile-like MNT products from 136 Xe+ 198 Pt collisions at the VAMOS spectrometer at GANIL [265]. VAMOS was positioned at an angle of 30 • , close to the grazing angle and supplemented by the EXOGAM array [271,272] around the target to measure prompt γ rays. The A and Z resolutions were ΔA/A = 1/200 and ΔZ /Z = 1/60, allowing for identification of Xe-like fragments. The isotopic distributions of target-like MNT products were then deduced from the measured distributions of projectile-like nuclei.
So far, the results of all direct and indirect studies of MNT products below Pb are in good agreement and indeed indicate that MNT cross-sections overtake fragmentation crosssections toward the neutron-rich side and toward smaller proton numbers of the reaction products. But at this point it is important to consider that the experimental feasibility is not only determined by cross-sections but finally by the resulting yields which depend on experimental conditions like beam intensity, applicable target thickness and efficiency of the  [264]; the denoted cross-section values are corrected for the angular efficiency of SHIP, which was 0.5% in this case. Open circles denote data from a thick target experiment with 64 Ni + 208 Pb performed at Legnaro by Krolas et al. [262]. Open triangles represent the data on Pt isotopes measured by Barrett et al. at ANL in 136 Xe + 208 Pb collisions, equally in a thick-target experiment [260]. The cross-sections measured in fragmentation reactions [266,267] are shown by asterisks. In fragmentation reactions, isotopes with still several more neutrons were discovered so far (see chart in Fig. 29 and Refs. [268][269][270]), but they are not all displayed here. The N = 126 shell for Os, Ir, Pt, Au, Hg, and Tl is located at mass numbers A = 202, 203, 204, 205, 206, and 207 in the respective distributions experimental setup. And these are quite different in MNT and fragmentation reactions. In fragmentation reactions at relativistic energies, targets of 10 g /cm 2 can be used, while in MNT reactions at the Coulomb barrier target thicknesses are about 1mg /cm 2 . In Fig. 31, we show the expected yields for MNT and fragmentation products, taking the cross-sections from Fig. 30 and using the parameters for beam current and target thicknesses given in the figure caption. Note that the yields in Fig. 31 are given at the target and do therefore not yet include the efficiencies of the experimental setups which can vary strongly between the different experiments. The overall trend indicates that the expected yields seem in favor of fragmentation reactions, mainly given by the possibility to use 10,000 times thicker targets. Apart from this, projectile fragmentation reactions combine further advantages, which make them much more efficient than MNT reactions: (i) projectile fragments are emitted in a very narrow forward cone, enabling very efficient and fast In-flight separation, while MNT products are emitted in a wide angular range which makes their collection/separation ineffective; (ii) A and Z identification of the relativistic fragments can be performed with the E-ΔE-TOF-Bρ method which makes identification independent of decay properties; also, the method is very sensitive and is in principle applicable for a single event. In contrast, the lack of appropriate and sensitive identification techniques for heavy MNT products is still a serious bottleneck which makes the method inefficient and cumbersome.

Shell effects in MNT reactions
Experiments on DIC, many of them performed in the 1980s, indicate that shell effects play an important role during the DNS evolution [62,[86][87][88][89][90][91][92] and influence the yield distributions of deep-inelastic reaction products. For example, in the symmetric reaction system 238 U + 238 U [93] one observes an unexpectedly intense diffusion of protons compared with other symmetric systems like 208 Pb + 208 Pb [94]. The variances of mass distributions in 154 Sm + 154 Sm and 144 Sm + 144 Sm reactions strongly differ, which is explained in [95] by the large neutron binding energy in the latter reaction. The cross-sections for compound nucleus formation in the reactions 100 Mo+ 100 Mo and 110 Pd+ 110 Pd differ by four orders of magnitude, although the masses of the interacting nuclei differ by only 10 units [96,97].
The asymmetric systems 16 48 Ca + 238 U reaction, a large mass variance is observed along with a small change in the mean mass number of the light fragment, but it is not the case in the 40 Ca + 238 U reaction [67,68].
Reactions with 238 U beams and targets of 110 Pd and 124 Sn were studied at GSI in the mid 1980s [273,274]. In both cases the bombarding energy was about 7% above the corresponding Coulomb barrier. Strong nuclear structure effects were not only found in the initial phase of the reaction with small energy transfer but also in the final phase with large energy dissipation. In 238 U+ 110 Pd reactions at 5.7 MeV/nucleon, the nucleon flow toward symmetry was strongly oriented along the minimum of the driving potential in the early stage of the reaction until the U-like MNT product reached the Z P ≈ 82, N P ≈ 126 shells and the target-like nucleus, nearly simultaneously, reached the Z T ≈ 56, N T ≈ 82 shells. After this, mainly protons flow from the heavy partner to the lighter one by keeping the neutron numbers of the two nuclei close to magicity. In the initial phase of the reaction 238 U + 124 Sn at 6.02 MeV/nucleon, equal numbers of protons and neutrons were transferred. However, starting from Z P ≈ 82, N P ≈ 135 and Z T ≈ 60, N T ≈ 83, mainly proton transfer was observed as in the 238 U + 110 Pd case. A periodic variation of the multiproton transfer probability as a function of the nuclear charges of the residual light nuclei was also observed in both reactions [274]. Thus, the neutron transfers in these studies reflect driving potentials towards the N T = 82 and N P = 126 closed shells [273,274]. In other words, the driving potential carries information about the paths of evolution of the DNS. Maxima in the mass and charge yields arise from minima in the driving potential and are caused by shell effects in the DNS nuclei [20,29,160,167].
Also in the new generation of theoretical and experimental studies shell effects are discussed as an important factor for the yield distributions of MNT products. In the Langevin model, shell effects in very heavy systems like 238 U + 248 Cm lead to a pronounced minimum at Z P = 82, N P = 126 in the driving potential and enhanced cross-sections for Pb-like MNT products. As a consequence, also the cross-sections of the complemenatary superheavy target-like MNT product (Z T = 106 in the case of U + Cm collisions) will be enhanced [33]. New experimental data on shell effects in heavy systems are still scarce. Our studies of 238 U + 238 U collisions did not reveal sufficiently strong signatures for shell effects to draw convincing conclusions [75]. Only, the measured isotopic distributions indicate that MNT products within the environment of four protons around Pb are created with lower excitation energies than more distant nuclei. Similar effects, but more pronounced, we observed in collisions of 58,64 Ni + 207 Pb at the velocity filter SHIP [71]. The influence of shell effects, particularly of the closed neutron shell at N T = 126, on the excitation energy of the transfer products was seen at the lowest energies, directly at the Coulomb barrier and faded with increasing beam energy. Also, we observed in different collision systems ( 58,64 Ni + 207 Pb, 48 Ca + 238 U, 48 Ca + 248 Cm) at SHIP that the cross-sections of MNT products in the region around Pb decreased only slightly, and the slope became steeper when the products move away from the Z ≈ 82 shell [275]. Similar observations are reported in Refs. [78,79,[276][277][278] for the systems 48 Ca + 232 Th, 238 U, 244 Pu, 248 Cm as well as for 36 S, 48 Te, 64 Ni + 238 U, 136 Xe + 208 Pb, 88 Sr + 176 Yb and 156,160 Gd + 186 W.
All these results suggest that a careful choice of the collision system seems very important to exploit shell effects for enhancing the cross-sections of exotic nuclei in distinct regions of the nuclide chart.

Transfer-type reactions at intermediate energies
Besides MNT reactions at low energies, heavy-ion collisions at intermediate energies is a well established method for the production of rare isotopes [279][280][281][282][283][284][285][286][287][288][289][290][291]. As shown in Refs. [64,292], in the 209 Bi + 136 Xe and 197 Au + 208 Pb reactions at the lower boundary of the Fermi energy domain the total reaction cross-section is almost entirely accounted for by binary collisions, irrespective of a possible further disassembly of the two highly excited primary fragments. Dissipative binary dynamics at these bombarding energies has been also observed in Refs. [293][294][295]. The observed influence of the target isospin on the final isospin of the projectilelike fragments and broadening of the projectile-like fragments charge distributions with increasing energy dissipation is similar to that encountered in MNT reactions at lower bombarding energies up to a few MeV/nucleon above the Coulomb barrier. The transfer process is found in Refs. [296][297][298][299][300][301][302][303][304][305] to be a very strong component for the projectile-like products in peripheral collisions. In the following collision systems nucleon pickup products have been observed among the products of projectile fragmentation reactions at bombarding energies above the Fermi energy: 48  As well it was observed in reactions of 40 Ca, 48 Ca, 58 Ni, 64 Ni beams at 140 MeV/nucleon on 9 Be and 181 Ta targets [308]. As shown in Ref.
[289], in the collisions of nuclei 48 Ca+ 9 Be, nat W at incident energy 142 MeV/nucleon the yields of the most neutron-rich isotopes of light nuclei tend toward the Q gg -systematics [57,62]. This allows us to assume the binary character of the interaction, contributing to the production of neutron-rich nuclei in the reactions 48 Ca + 9 Be, nat W. The larger yields of neutron-rich nuclei with nat W targets compared to 9 Be targets indicates the strong contribution of reaction types other than fragmentation. In fragmentation reactions the sequential evaporation of light particles from the strongly excited nucleus leads to a somewhat uniform distribution of the final products which underlies the semi-empirical EPAX systematics based on the data from many high-energy experiments [309]. The disagreement of the yields of neutron-rich nuclei with the EPAX formula [309] supports the assumption that MNT binary reactions play an important role in the production of exotic nuclei even at quite high bombarding energy. The mechanism of this reaction seems to be the same like the mechanism of DIC. The collisions should occur with large angular momenta (large impact parameters) to supply small excitation energy in the neutron-rich products.
Nucleon transfer reactions at incident energies around 100 MeV/nucleon must occur in nearly peripheral collisions to avoid high excitation in the DNS and fragmentation. The excitation energy of the DNS with an exotic nucleus (Z , N ) is E * (Z , N ). Assuming thermal equilibrium, the excitation energy of the light nucleus with mass It is clear that the formation probability of the DNS with exotic nucleus (Z , N ) increases with E * (Z i , N i , J ). However, this increase is only possible up to the moment when E * P (Z , N ) becomes equal to the neutron separation energy S n (Z , N ). Further increase of E * (Z i , N i , J ) would lead to a strong loss of neutron-rich nuclei due to neutron emission. Taking E * P (Z , N ) ≈ S n (Z , N ), we find the optimal value of E * (Z i , N i , J Z ,N ) and the corresponding angular  N ), and the contribution of collisions with J > J N ,Z to σ Z ,N decreases with increasing J because the formation probability Y Z ,N ,J decreases. Therefore, where μ is the reduced mass of the projectile-target combination and ΔJ is the angular momentum interval above J N ,Z which mainly contributes to the cross-section. In our calculation we set ΔJ = 20 which corresponds to a change of the impact parameter of less than 0.2 fm at the incident energies considered. Only a narrow region of partial waves contributes to the transfer cross-section, strongly selecting the initial condition.
To estimate Y Z ,N ,J , the simple statistical method is used. The diffusion in mass and charge asymmetries is important in the DNS evolution. The simultaneous investigation of the diffusion in these collective coordinates allows us to calculate the formation probability Y Z ,N ,J of the DNS configu-ration. We assume that the average interaction time of two nuclei is much larger than the transient times in mass and charge asymmetries. We approximate the expression for the formation probability rate with the Kramers-type formula (33) [310] and obtain and the average interaction time t 0 = 5 MeV −1h ≈ 3 × 10 −21 s of the two nuclei is assumed to be equal to the characteristic time of DIC. Note that the decay of the DNS from the initial configuration is the dominant decay channel here. As seen, the value of ΔB Z ,N ,J Z ,N contains the corresponding Q-value to expect the Q gg -systematics for the isotopic distribution.
The production cross-sections of neutron-rich isotopes calculated within our approach in the reactions 48 Ca + 181 Ta, nat W at intermediate incident energies are listed in Table 2. The comparison of these results with available experimental data has sense only for those neutron-rich isotopes which are mainly produced as primary products of binary reactions. To the yields of lighter isotopes there are contributions of fragmentation processes as well as of de-excitation by neutron emission of heavier primary isotopes. In order to describe these yields the extension of the present model is necessary.
The calculated results are in good agreement with most of the experimental data. This strongly supports the proposed model. Our calculations clarify that the MNT process is the main process which contributes to the total reaction yields of most of the exotic nuclei in the intermediate energy region. The observed new isotopes with neutron numbers larger than the projectile neutron number in the above named reactions are most probably transfer products from collisions at large angular momentum (peripheral collisions) leading to small excitation energy in the primary neutron-rich products. At J < J Z ,N , the primary neutronrich nuclei are excited and transformed into the secondary nuclei with less number of neutrons due to de-excitation by nucleon emission. The yield of these secondary nuclei follows the Q gg -systematics as well because of the binary character of the reaction. Indeed, in Ref.
[289] the Q ggsystematics fit well the yields of various isotopes. The value of the angular momentum in the entrance channel governs the competition between fragmentation and massive transfer processes.
Instead of the Q gg -systematics, a systematics based on the binding energy per nucleon of the neutron-rich isotope is suggested in Ref. [311]. The binding energy of the neutron-rich isotope correlates with the Q gg -value because the mass excess of the conjugated heavy fragment weakly changes with mass number. The cross-sections calculated for very neutron-rich isotopes with the EPAX [309] 48 Ca projectiles [25] are larger than the experimental ones. The EPAX model describes well the yields of the isotopes of Mg and Si with N − Z < 10. It is apparent that in the binary reaction the projectile must be as close as possible to the region of the nuclide to be produced because in this case a smaller number of nucleons has to be transferred.
Since the predicted production cross-sections for the new exotic isotopes 47 P, 51,53,55,57 Cl, 52,54 Ar, 56,58,60 Ca, 59,61,63 Sc, and 62,64,66 Ti are larger than 0.1 pb, they can be synthesized and detected with present experimental possibilities. The predicted cross-sections seem to be optimistic, specially for the isotopes of Ca, Sc and Ti, in the sense that the predictions are done by assuming that the excitation energy of the DNS is divided proportionally to the mass numbers of the fragments. In transfer reactions the excitation energy would be preferentially generated in the primary pickup products with N + Z > N i + Z i . One should also mention that the production cross-section weakly depends on the bombarding energy. For example, in the reactions 48 Ca + 181 Ta the cross-section at beam energy 64 MeV/nucleon is about of 5% larger than at 140 MeV/nucleon. This is because of the very weak dependence of the ratio (2J Z ,N + 1)/E c.m. on E c.m. in Eq. (58).
One can see that the cross-sections in 48 Ca + nat W reactions are larger than the corresponding cross-sections in reactions of 48 Ca + 181 Ta. Irradiating heavier targets by 48 Ca beams for producing neutron-rich isotopes, we gain in the Q gg value as well as in the value of ΔB Z ,N ,J Z ,N . Therefore, heavier targets are preferable for the production of neutron-rich nuclei. For example, replacing 181 Ta or nat W by 232 Th or 238 U or 248 Cm, one can increase the yield of neutron-rich isotopes. This effect should be taken into consideration in the planned experiments. Summarizing, we see that MNT reactions can occur also at intermediate energies and overlap with the fragmentation reactions. This leaves the question which energy is optimal for the production of new isotopes in MNT reactions.

Incomplete fusion of light weakly bound nuclei
To complete our considerations about nucleosynthesis in transfer reactions, we will discuss here the situation if one of the reaction partners is a very light nucleus. Understanding fusion dynamics of light weakly bound nuclei is very important for describing astrophysical reactions important for nucleosynthesis [312]. The breakup mechanism of weakly bound nuclei is crucial to reveal the dynamics of fusion as well as the consequences of breakup [313,314]. There are three possibilities of reaction processes upon breakup of a weakly bound nucleus. The first is where no fragments are captured and is termed non-capture breakup (ncbu). The second is where not all fragments are captured, which is termed incomplete fusion. The final possibility is where the light nucleus is captured completely by the target nucleus and is termed complete fusion [315]. An important aspect is the interplay between breakup and other reaction processes like the transfer process. However, due to the transfer process's similarities with the icf process it is hard to separate them from an experimental point of view as the fusion products from both processes are the same. Transfer can also cause the breakup of weakly bound nuclei during low-energy collisions [316][317][318][319][320][321][322][323][324].
Various theoretical approaches for understanding the incomplete-fusion process have been developed in the past, spanning a range of concepts and considerations, including breakup fusion, angular momentum window for incomplete fusion, promptly emitted particles, Fermi-jet, exciton, and moving source, thereby explaining the measured energy spectra and angular distribution of the emitted fragments and population of angular momentum in the compound nucleus [325]. New types of models have been used to address low-energy fusion dynamics of weakly bound nuclei, ranging from classical to quantum-mechanical methods. Reference [315] provides a critical survey of different theoretical approaches. New studies on the inclusive non-elastic breakup cross-section may provide a quantum-mechanical route to the calculation of the icf cross-section of weakly bound nuclei [326,327]. Another interesting quantum-mechanical framework is the time-dependent wave-packet (tdwp) method [329,330]. This method calculates the incomplete-and complete-fusion cross-sections unambiguously [330], which is a challenge using the continuum discretized coupledchannels method [331][332][333]. The tdwp approach is currently undergoing further development to be implemented using a three-dimensional reaction model. Some of the challenges of the quantum-mechanical models can be overcome via the use of the three-dimensional classical dynamical model [334][335][336]. This model is implemented using the platypus code [336], which uses classical trajectories in conjunction with stochastic breakup [334,335]. This is done through the input, which includes a breakup function [335,337], determined from sub-barrier breakup measurement [320,322], that undergoes Monte-Carlo sampling [334,335]. This breakup function encodes the effect of the Coulomb and nuclear interactions that cause the breakup, making this approach a quantitative dynamical model for relating the sub-barrier ncbu to the above-barrier incomplete and complete fusion of weakly bound nuclei, rather than a breakup model [334][335][336]. In contrast to most existing models for incomplete fusion, platypus treats the dynamics of incomplete fusion and provides a number of differential cross-sections that are critical for understanding exclusive experimental data [338]. This approach has very recently been extended to the incomplete fusion of complex projectiles [339]. It is also important to note that this fusion model only works at energies above the Coulomb barrier between the projectile and target. This is due to the absence of quantum tunneling that is the primary way of fusion at suband near-barrier energies. There has been a recent attempt to amend this classical model by adding a correction at sub-and near-barrier energies, to take into account quantum tunneling. This was done by incorporating a tunneling factor based on the wkb approximation [340]. This improved the results outputted from the model, relative to experimental sub-barrier fusion measurements [340]. Additional modifications have recently been suggested for interpreting sub-barrier breakup measurements [323,324]. The role of (i) prompt and delayed, direct breakup, and (ii) transfer-triggered breakup modes on incomplete fusion of 6 Li + 209 Bi collisions at energies above the Coulomb barrier has been addressed with an extended version of platypus [341]. For instance, prompt breakup happens in the instant the excitation of the 6 Li projectile is chosen to take place. At this point 6 Li is converted into its cluster fragments (alpha-deuteron) and then the fragments and target propagate according to the defined interactions between them [335]. Delayed breakup is induced by reaching the 1 + , 2 + or 3 + resonant states in 6 Li, which then triggers the dissociation of 6 Li with the delay coming from the half-life of the resonant 6 Li state. The 3 + resonant state has a much longer half-life than the 2 + and 1 + resonant states, so the 6 Li breakup takes place at the outgoing branch of its trajectory, far away from the target nucleus, not affecting fusion. So the effect of the 3 + resonance on fusion can be neglected. Figure 32 shows the incomplete-fusion excitation function in collisions of 6 Li + 209 Bi. A number of reaction processes contributes to the formation of incomplete-fusion products, the dominant process being the neutron-stripping channel (see solid and dashed lines) that involves the breakup (delayed and/or prompt) of the projectile-like nucleus 5 Li [341]. In contrast, most quantum-mechanical fusion calculations assume that the direct breakup of 6 Li into α and deuteron is the dominant incomplete-fusion channel, which is not the main observed breakup channel [322][323][324].  [342] are compared with platypus calculations at above-barrier (arrow) energies. Direct breakup as well as transfer-triggered breakup channels are included. a For delayed direct breakup of 6 Li, delayed breakup of 8 Be after d-pickup, and delayed breakup of 5 Li after n-stripping. b The same but for prompt breakup processes. Clearly, the neutron-stripping process dominates the formation of icf products (i.e., actinium and polonium isotopes) The 2n-stripping channel is important for the production of 211 Bi in collisions of 7 Li + 209 Bi at Coulomb energies [342], as shown in Fig. 33. This process also involves the breakup (delayed and/or prompt) of the projectile-like nucleus 5 Li. Although the 211 Bi production is explained by this extended version of platypus [341], it fails in explaining the incomplete-fusion excitation function for 7 Li+ 209 Bi, as shown in Fig. 34. This is due to the large contribution of the triton-transfer channel that cannot be treated within the  [342] are compared with PLATYPUS calculations at above-barrier energies. Theoretical calculations [341] include either prompt or delayed breakup of both 7 Li and projectile-like nuclei ( 8 Be, 6 Li and 5 Li) after single-particle transfer processes. The triton-transfer channel from the 7 Li ground state cannot be treated within the PLATYPUS model, but it is very important for the formation of specific incomplete-fusion products [328] 7 Li + 198 Pt are compared with platypus + pace2 calculations (lines) [338] platypus model. Very recent measurements [328] have confirmed the importance of this channel for the yield of specific incomplete-fusion products. Figure 35 displays the production cross-section of some evaporation residues of (a) gold and (b) mercury isotopes in the incomplete fusion of 7 Li + 198 Pt [338], after the direct breakup of 7 Li into triton and α fragments. Particleγ coincidence measurements in conjunction with platy-pus+pace2 calculations have provided insights into the formation mechanisms of these incomplete-fusion-evaporation residues [338], suggesting that the breakup fusion mechanism drives their formation. In case of 199 Hg, the γ -ray transitions only above the (13/2 + ) isomeric state were considered, therefore the measured cross-sections only provide a lower limit for this channel [338].
The present results for weakly bound stable 6,7 Li beams are useful for the theoretical developments, and also have implications in predicting production cross-sections of exotic nuclei with radioactive ion beams.

Nucleosynthesis with radioactive ion beams
The application of neutron-rich RIBs would allow one to use fusion-evaporation reactions, and thus to profit from the well defined forward oriented kinematics of the fusionevaporation residues. Besides, fusion reactions with neutronrich RIBs could be a further alternative to reach nuclei along the N = 126 shell. Another possibility would be to use neutron-rich RIBs in MNT reactions [35,343]. Accelerator facilities which provide intense RIBs at Coulomb barrier energies are still in the process of arising. Nevertheless one can make some feasibility considerations based on model calculations. To prepare the discussion, we give in the next section an overview on the state of the art in RIB technologies and the yields which we could expect in such experiments.

Radioactive ion beams: production and yields
Principally, there are two different techniques of exotic ion beam production: the isotope separation on-line (ISOL) [344] and the In-flight method [345]. There are numerous review articles about RIB technologies to which we refer for details while we limit ourselves here to a short overview. The basic principles of the ISOL and In-flight technique are sketched in Fig. 36. The ISOL technique uses mostly proton beams in the 1 GeV energy range to induce fragmentation or fission reactions in target nuclei of heavier elements up to uranium. Characteristic for the ISOL technique is the application of thick targets in which the reaction products are stopped. The targets are usually part of a combined target-ion-source system for production, extraction and ionization of the produced nuclei. From the target-ion-source, the ions are extracted with voltages of several 10 keV. A subsequent magnetic separator selects RIBs of the desired species. The ISOL technique is the oldest technology for RIB production [344]. A typical representative of an ISOL facility is the facility at CERN. The Inflight technique is in some sense complementary to the ISOL method. It uses the principle of projectile fragmentation or fission for which projectiles up to uranium are accelerated to energies of several 100 MeV/nucleon up to 1 GeV/nucleon Due to the different techniques, the elemental and isotopic yields from ISOL and In-flight production are different (Figs. 37 [346], 38 [347]). The principal difference is that the ISOL technique cannot provide isotopes of all elements because the chemical properties of the elements determine how efficiently they can be ionized and extracted from the target-ion-source. Figure 37 reveals the gaps in the isotope (element) yields which emerge from this effect. The highest yields are obtained for heavy alkali elements and nobel gases. Also, the availability of very light RIBs Z <∼ 25) is quite restricted and their intensity is typically much higher in In-flight production. These restrictions do not concern the Inflight technique which can provide isotopes of all elements.

How to get RIBs to Coulomb barrier energies?
Fusion-evaporation reactions require RIBs at Coulomb barrier energy, which is around 5 MeV/nucleon for asymmetric collision systems leading to superheavy nuclei. But neither the ISOL nor the In-flight technique can directly produce RIBs with these energies. Therefore, post-acceleration or deceleration, respectively, is necessary. In the following we discuss established and thinkable techniques and estimate the   [347] beam losses which are connected with them. The individual methods are depicted in Fig. 39.

Post-acceleration of ISOL beams
Radioactive ion beams produced with the ISOL method must be post-accelerated to bring them to Coulomb barrier energies (see scheme in Fig. 39). A classical example for such a facility is ISOLDE at CERN which delivers since recent years RIBs with energies between (5-10) MeV/nucleon by using the post-accelerators REX-ISOLDE [348] and HIE-ISOLDE [349]. For efficient post-acceleration, charge breeding is necessary. But before injection to the charge breeder, the mass separated RIBs with original charge state 1 + are accumulated, cooled and bunched in a Penning trap. After extraction from the Penning trap, the ion bunches are injected to the charge breeder, where they are bombarded with electrons to further ionize them to higher charge states. This procedure is responsible for the main intensitiy losses of ISOL beams. For example, the trapping and charge breeding system at ISOLDE has efficiencies between 5% and 20%, depending mainly on the ion mass [350]. It reaches a maximum around A = 50 and decreases toward larger masses. One has also to consider that cooling and charge breeding needs time. Cooling times are typically longer than ∼ 10 ms. The breeding time varies between few milliseconds and 500 ms, also depending mainly on the ion mass. Therefore only ions with sufficiently long half-lives survive the procedure. Roughly, the half-life of the nuclei should at least be on the order of seconds to avoid additional losses due to radioactive decays before reaching the production target. Postaccelerated ISOL beams are of good quality, having small emittance and energy spread.

Deceleration of In-flight beams
Deceleration with degraders The simplest method to slow down relativistic RIBs is the application of degrader foils. For RIBs with energies of several 100 MeV/u, one needs degrader thicknesses of at least 10 g/cm 2 . Interesting is to see the impact of such thick foils on beam emittance and energy spread. Respective studies were performed for the low-energy branch of the future Super-FRS facility [351,352] at FAIR/GSI. Figure 40 shows the simulation Fig. 40 Simulation of a beam energy distribution, b beam spot size and c opening angle distribution after decelerating a (radioactive) oxygen beam with degrader foils from 190 MeV/nucleon to 6 MeV/nucleon [353]. For details see text results for the example of oxygen RIBs with an energy of 190 MeV/nucleon [353]. In this example, four thinner degraders instead of a thick one were installed in different places along the beam line through the Super-FRS and subsequent energy buncher. The degraders were assumed as aluminum wedges with thicknesses of 3.7 g/cm 2 , 1.8 g/cm 2 , 3.46 g/cm 2 and 0.63 g/cm 2 which decelerate the beam stepwise to 6 MeV/nucleon. Figure 40a shows the energy distribution of the oxygen ions after the last degrader. It is peaked around 6 MeV/nucleon with a full width at half maximum (FWHM) of 4 MeV/nucleon (ΔE/E = 67%). The excitation functions of hot fusion reactions leading to superheavy nuclei have relative widths of ΔE/E ≈ 20% at FWHM which means that only ∼ 10% of the beam intensity can be effectively used for the synthesis reaction. Another factor 10 of beam intensity is lost due to small angle scattering in the degraders in course of the deceleration process. Finally, the RIB intensity which is effectively available for the synthesis reaction is ∼ 1% of the original RIB intensity. Besides, the beam spot size and emittance are large: the beam diameter at the production target is 10 cm (!) and the opening angle is 100 mrad (Fig. 40b, c). The deceleration with degraders is a fast method, where the lower time limit is only given by the flight time of the RIBs through the separator, but it leads to bad beam quality.
Gas stopping and post-acceleration Another possibility to decelerate In-flight RIBs is their stopping in a gas catcher and successive re-acceleration after extraction. This method provides decelerated RIBs of good quality. A respective technique has so far been developed at the Facility for Rare Isotope Beams (FRIB) at the National Superconducting Cyclotron Laboratory (NSCL) of Michigan State University [354]. The bottlenecks in such systems are the stopping and extraction efficiencies of the ion catcher and the efficiency of the charge breeding which is performed at NSCL with an electron beam ion source. With state-ofthe-art ion catchers, stopping and extraction efficiencies of some 10% can be reached. The charge breeder efficiencies are on the same scale. At NSCL it is planned to replace in the future the gas-stopping cell by a gas-filled reverse-cyclotron [354,355].
Deceleration in a storage ring Another imaginable possibility to obtain Coulomb barrier beams with good quality could be the deceleration of In-flight RIBs in a storage ring (Fig. 39). Already existing storage rings for exotic beams like the ESR [356] at GSI Helmholtzzentrum or HIRFL-CSR [357] in Lanzhou would be predestined for this method. The GSI ESR, for example, can decelerate relativistic ions down to energies of 4 MeV/nucleon. During the deceleration process, beam cooling is necessary, which is performed by electron cooling. With this, the complete deceleration process takes about 10 s, which means that the method is suitable for ions with lifetimes larger than roughly 10 s. The bottlenecks in this method are beam losses during injection in the storage ring and during the deceleration and cooling process leading in worst cases to efficiencies on the 1% scale. But technical improvements are always possible.
If one carries the thought further, one could even imagine to let the decelerated RIBs circulate in the storage ring and let them traverse in each turn a thin target installed in the ring. At Coulomb barrier energies, the beam performs about 10 5 turns per second in a ring with circumference of about 100 m, like the ESR. This would effectively enhance the RIB intensity by the same factor 10 5 . There are so far no feasibility studies how realistic such a scenario could be. The critical point is the target, which leads to energy loss and deterioration of the beam emittance, followed by beam losses. Figure 41 shows the superheavy element region of the Chart of Nuclides with the presently known isotopes. Besides, it contains three areas of major interest which are still empty or not confirmed, respectively. One of these areas are the gaps between nuclei produced in cold and in hot fusion reactions. Filling these gaps would join the presently isolated region of relatively neutron-rich superheavy nuclei to the remaining, well established part of the nuclide chart. To fill these gaps, one does not even need RIBs because nearly all of them can be synthesized in fusion-evaporation reactions with stable beams, for example by using 248 Cm targets and the beams indicated on the right side of the chart. Another area of interest are the endpoint nuclei of the decay chains of superheavy isotopes from hot fusion reactions. Nuclei below Z = 112 in these chains were so far never produced directly and their direct synthesis would substantiate their correct assignment as well as the one of their mother nuclei. And finally, there is the "island of stability" with the predicted but still not confirmed new spherical shell closures in the superheavy element region. Concerning the proton shell closure, the model predictions diverge: Z = 114, 120 or 126 are suggested, but the location of the neutron shell is found by most theoretical models at N = 184 [358][359][360][361][362]. Figure 41 reveals that Fig. 41 Chart of Nuclides with the known isotopes of superheavy elements. The blue background represents shell correction energies (i.e. fission barriers) calculated in the macroscopic-microscopic model [358]. The chart contains three still empty, but very interesting regions: (i) the island of enhanced stability around Z = 114, N = 184, where spherical shell closures are expected; (ii) the endpoint nuclei of the decay chains of the most neutron-rich known isotopes, which were so far not produced directly and (iii) nuclei which fill the gaps between isotopes produced in hot and cold fusion reactions. To fill areas (i) and (ii), fusion reactions with radioactive ion beams or MNT reactions must be applied, while nuclei in area (iii) can be produced in fusion reactions with stable beams. The violet squares mark isotopes which can, for example, be produced in fusion reactions of 248 Cm targets with stable isotopes of the projectile nuclei indicated on the right side the most neutron-rich known isotopes are still seven neutrons away from N = 184. The previous section revealed that MNT reactions are realistic to produce several endpoint nuclei of decay chains from isotopes produced in hot fusion reactions, but they will very likely not allow one to synthesize nuclei on the island of stability. There remains the question if fusion reactions with RIBs would allow one to access it.

Expected cross-sections and yields
To get now an idea about the yields of superheavy nuclei which can be expected from fusion reactions with RIBs, we start with some basic considerations: All experiments on superheavy element synthesis with stable beams revealed so far that complete-fusion reactions with as asymmetric as possible projectile-target combinations lead to the largest evaporation residue cross-sections [104,105], because in this case the Coulomb repulsion between the reaction partners is minimized. There is currently no serious reason to assume that the situation would change if radioactive projectiles are used instead of stable ones. Principally there are, like in fusion reactions with stable beams, two possible approaches: (i) cold fusion reactions using Pb or Bi targets and (ii) hot fusion reactions using actinide targets. In cold fusion reactions, the evaporation residues are quite neutron-deficient. Also, the experiments with stable beams revealed so far that nuclei above Z = 112 are produced with significantly larger crosssections in hot fusion reactions. Therefore, we concentrate in the following representative examples on hot fusion reactions with 248 Cm targets, which is one of the heaviest and most neutron-rich available target material.
We first take a look on two representative model predictions of fusion-evaporation residue cross-sections. The first example in Fig. 42 shows DNS model cross-sections of Moscovium isotopes (Z = 115) from complete-fusion reactions of radioactive 47 K beams with 245,247,248 Cm targets [363]. The most neutron-rich fusion products are 291 Mc and 292 Mc, resulting from 4n and 3n evaporation channels. These are just the next neighbors of the already known isotope 290 Mc [364], and still seven neutrons far from the N = 184 shell. The maximum cross-section is about 6 pb, which is on the same order as cross-sections measured for Mc isotopes with stable projectiles. It has to be remarked that also the adiabatic model finds very similar cross-sections for superheavy nuclei in the same region [26]. The second example (Fig. 43) shows excitation functions of the neutron-rich rutherfordium isotopes 266,267 Rf produced in collisions of 22 O beams with 248 Cm targets. These calculations were performed with the adiabatic model which predicts maximum cross-sections of 50 nb and 5 nb, respectively [26]. The isotopes 266,267 Rf are the endpoint nuclei of the decay chains of 282 Nh [365] and 291 Lv [366] and were so far not produced directly.
The crucial parameter which decides about the experimental feasibility is the available RIB intensities. To get an impression, we show in Table 3 the expected average measurement times which would be necessary to produce one   Table 3 indicate some trends which are confirmed if one looks to further example reactions. We want to stress that the numbers which we gave above do not yet take into account the intensity losses of the RIBs during acceleration (or deceleration, respectively) to Coulomb barrier energies. The production yields of superheavy nuclei around the expected "island of stability" appear to be out of feasibility over the long term. But the synthesis of "light" superheavy nuclei with RIBs seems indeed realistic, particularly at future RIB facilities which will provide much higher beam intensities. By trend, inflight facilities seem best suited for such experiments because the intensities of the needed very light RIBs (roughly with A < 25) are much higher in In-flight production compared to ISOL. We want to stress that there is other work which considered this topic and arrived at very similar conclusions (e.g., [367,368]).

An "emergency solution"?
The question arises if there is an "emergency solution" to learn about possible new shell closures in the region of superheavy elements. It can be answered with "maybe". Namely, by studying QF and FF reactions, which have significantly larger cross-sections than fusion residues, a possible influence of shell closures could be explored. The presumption is that shell closures of the compound system might also be reavealed in the mass, angle and energy distributions of QF and FF fragments. So far, the only appropriate facility to perform such experiments is HIE-ISOLDE at CERN, which can provide since recently RIBs with energies up to 10 MeV/nucleon. Concerning the choice of projectile ion, one is also in such experiments bound to RIBs with intensities not lower than ∼ 10 6 /s in order to achieve statistical relevant data. This leaves 95 Rb as the most suitable projectile. In combination with 209 Bi targets, 95 Rb allows one to reach a compound system with proton and neutron shells at Z = 120, N = 184. An approved proposal for this experiment at HIE-ISOLDE is existing [369].

Existing and upcoming RIB facilities
There are activities in accelerator laboratories worldwide to upgrade existing RIB facilities or construct new ones to deliver high intensity exotic ion beams. The chart [370] in Fig. 44 gives an overview on presently existing ISOL and In-flight facilities as well as on facilities which are in planning or construction. In the following, we will pick out and describe a few of them to create an impression about the bandwidth of different technical approaches and potentials for future experiments, specifically in context with DIC at Coulomb barrier energies.
The "mother" of ISOL-type facilities is CERN ISOLDE which is in operation since 1967. Presently, ISOLDE is the only RIB facility which can provide intense exotic ion beams at Coulomb barrier energies up to 10 MeV/nucleon by using the post-accelerators REX-ISOLDE [348] and HIE ISOLDE [349] (Fig. 45). These energies enable nuclear reactions also in very heavy collision systems. Typical for ISOL facilities, the highest yields are obtained for heavier alkali metals like Rb or Cs, and noble gases like Xe, for a broad variety of their isotopes from neutron-deficient to neutron-rich. With presently available beam intensities and high-quality ISOLtype beams, ISOLDE provides a good potential for studying DIC and MNT reactions with exotic projectiles (see also Sect. 5.3.2).    [373] with its three experiment branches. The low-energy branch with the energy buncher spectrometer at its end is also foreseen to provide RIBs at Coulomb barrier energies for DIC studies. It is planned to decelerate the relativistic RIBs by using degraders A typical representative of an In-flight facility is the fragment separator (FRS) facility [372] at GSI/FAIR which is in operation since 30 years and uses primary beams with energies up to 1 GeV/nucleon. About 300 new isotopes were discovered at FRS so far. Presently, the new-generation RIB facility SuperFRS is under construction which is foreseen to deliver up to 100 times higher beam intensities compared to the present facility. Three branches for experiments are planned at the exit of the SuperFRS (see Fig. 46). The lowenergy branch, which will include an energy buncher, can also be used for DIC experiments at Coulomb barrier energies. At SuperFRS it is considered to decelerate the beams by using degraders. This leads, however, to fairly bad beam quality like discussed in Sect. 5.2.2. But there might be a way out which was so far not studied more closely, namely, to use a storage ring for RIB deceleration. For such studies, even the present FRS which is connected to the experimental storage ring (ESR) would be predestined.
The first powerful new-generation In-flight facility which is already in operation, is the Radioactive Isotope Beam Factory RIBF [374] at RIKEN Nishina Centre, Japan. Presently . After separation, the RIBs can either be used at full energy or stopped in an ion catcher. After extraction from the ion catcher one can use the beam directly for experiments at very low energy or reaccelerate it to Coulomb barrier energies it uses primary beams up to uranium with maximum energies of about 350 MeV/nucleon. Like the future SuperFRS, also RIBFs isotope separator BigRIPS will use superconducting magnets. RIBF is presently the worlds most powerful In-flight facility and provides the highest primary beam intensities.
The Facility for Rare Isotope Beams FRIB is also an Inflight type facility which is under construction at the National Superconducting Cyclotron Laboratory (NSCL) of Michigan State University [354]. There will be three branches for experiments at FRIB (Fig. 47). After In-flight separation, the RIBs can either be used at full energy at the "fast beam" branch. Or they can be stopped in an ion catcher. After extraction from the ion catcher, the beam can directly be used for experiments at very low energy. Or it can be reaccelerated to Coulomb barrier energies to use it for MNT experiments and DIC studies. As mentioned already in Sect. 5.2.2, it is planned to replace in the future the gas-stopping cell by a gas-filled reverse-cyclotron for beam deceleration.
Finally, there are RIB facilities under construction which will offer both, ISOL-type as well as In-flight production of exotic ion beams. One such representative is SPIRAL2 [375] at GANIL in Caen, France.

The quest for new detection techniques
Heavy MNT and fusion-evaporation products are emitted with low energies, therefore their identification via the universal E-ΔE-TOF method fails due to the pulse-height deficit. For these nuclei, the most sensitive and effective identification technique is α decay tagging. The signatures of an α decay chain are mostly so clear that a single decay chain is sufficient to pin down the isotope. But it is obvious that the method is not universal and restricts to α emitters with appropriate half-lives and unambiguous decay chains. For example, it is not applicable for neutron-rich β emitters below Pb if one considers to produce them in MNT reactions. All respective experiments used so far the emitted γ rays for identification which is not nearly as sensitive as α spectroscopy; cross-sections below 1 µb were not reached in any of these experiments. Also, many superheavy nuclei toward the neutron-rich side cannot be identified by α decays. Therefore, the success of extending the upper part of the nuclide chart depends crucially on the availability of sensitive and preferably universal detection techniques, regarding the expected small cross-sections. In the following we will consider possible techniques which are presently in discussion.

High-precision mass measurements
Reaction products can be identified if their mass is measured with sufficiently high resolution to distinguish different isobars. For this, mass resolving powers of m/Δm ≈ (10 5 − 10 6 ) are needed like shown in the example in Fig. 48. Penning traps or multiple reflection time-of-flight mass spectrometers (MR-TOF-MS) are appropriate devices to reach these resolutions. In MNT reactions a broad region of isotopes is populated which results usually in the appearance of several isobars for each A, containing also already known ones. Therefore, the mass parabola will comprise also known masses beside the new ones and the location of the new mass on the parabola will identify the nucleus, even if Z is not directly determined. Investigations of such a method are ongoing at the GSI fragment separator facility. The setup consists of a cryogenic stopping cell [376] and an MR-TOF-MS [377] which are installed behind the fragment separator. The idea is to install targets with thicknesses of several 10 µm inside the stopping cell which is filled with helium gas to stop the reaction products. This allows one to collect reaction products emitted to a wide angular range, which is particularly important for MNT products. The relativistic beams from FRS are decelerated by degraders, such that its mean energy is about 10 MeV/nucleon when it impinges the reaction target inside the ion catcher. The stopped reac- tion products are extracted and injected to the MR-TOF-MS which allows for the broadband detection of different isotopes. For a sufficiently precise mass measurement, about 10 events are necessary. The method allows one to detect nuclei with half-lives ≥ 10 ms. In order to avoid space charge effects in the stopping cell, the maximum beam intensity is limited. Estimates for uranium beams arrived at maximum intensities of 10 7 ions/s. According to simulations, the method requires presently reaction product cross-sections of 100 µb or more [378].

Laser resonance ionization
Another approach is the Z identification of MNT products by laser resonance ionization. A respective system which focuses on neutron-rich beta-emitters below Pb was set up at the KEK Isotope Separation System (KISS) at RIKEN [379,380]. The MNT products are stopped in an argon gas filled ion catcher. The ion catcher is doughnut-shaped to prevent the primary beam from entering the cell. Laser resonance ionization inside the ion catcher is used for Z selection. After extraction from the gas cell, the ions are passing a magnetic dipole field for A/q separation. The selected ions can then be guided to a detection system to perform spectroscopy studies. The method is applicable for ions with lifetime of 1 s or more. This is mainly determined by the extraction times from the gas catcher which are about 0.5 s. The overall efficiency of the present KISS setup is on the scale of 0.1%.

Combination of mass separation and decay tagging
A combination of mass separation and decay tagging in the region of heaviest elements was for the first time successfully applied at LBNL (see [381] and the references therein). Isotopes of element 115 were produced in fusion reactions of 48 Ca + 243 Am. After passing the gas-filled magnetic separator BGS, they were stopped in a radiofrequency gas catcher. After extraction, ion bunches were created in an RFQ trap and from there injected in a trochoid spectrometer. The separated ions were then implanted in a silicon strip detector, where their α decay chains were measured. With this method it was possible to measure directly the A of 288 Mc and its decay daughter 284 Nh. Ion losses occur mainly during stopping and extraction from the gas catcher. The transport efficiency from the exit of BGS to the focal plane detector is denoted in [381] with 14%. With this, the method is applicable for production cross-sections on the picobarn scale. Transportation times through the system are on the scale of several 10 ms.

Summary considerations
The present chart of nuclides reveals still large empty areas which need to be filled. To our opinion, the "conventional" fragmentation, fusion and fission reactions will continue to provide the highest potential for laboratory nucleosynthesis. Most promising are the new powerful exotic ion beam facilities which are presently emerging in many places. They are designed to deliver up to three orders of magnitude higher beam currents which will allow one to extend considerably the chart of nuclides in the region below uranium by further using fragmentation and fission reactions. Particularly the sill sparse area of neutron-rich r-process nuclei will profit from it.
More critical is the region of transuranium nuclei which is presently populated solely in complete-fusion reactions. Particularly their limitation toward the neutron-rich side cannot be overcome readily because fusion reactions with stable projectiles lead to rather neutron-deficient residues due to the bending of the stability line. Some new isotopes can be produced in charged particle evaporation channels in the de-excitation of CN [382][383][384][385]. Fusion reactions with RIBs would be the most elegant way to produce neutron-rich transuranium nuclei, but the available intensities of suitable RIBs are still far too low and would mostly require unrealistically long measurement times. Regarding this situation, also the new generation of RIB facilities will not provide a satisfactory solution. Interesting is that about 50 new superheavy isotopes between Z = (104-118) can still be produced in fusion reactions with stable projectiles and actinide targets and would fill the present gap between the regions of cold and hot fusion-evaporation residues.
Aiming for a solution to synthesize more neutron-rich superheavy nuclei, we arrive at MNT reactions. The kinematics of MNT reactions is less favorable than fusion kinematics because it leads to broad angular and energy distributions of the residual nuclei. This makes their efficient separation from background events difficult and leads necessarily to larger losses of reaction products compared to fusion reactions. Nevertheless, the so far accomplished experiments and theoretical predictions indicate that one can expect a variety of new isotopes from MNT reactions. Toward the neutronrich side, one can realistically expect new isotopes of elements with proton numbers Z ≤ 108 and neutron numbers next to the already known ones. This will allow for the direct synthesis of the endpoint nuclei of hot fusion decay chains. Rather unrealistic seems from the present point of view the synthesis of very heavy MNT products with Z > 108. A large potential of MNT reactions we expect, however, for the region of very neutron-deficient transuranium nuclei, even if it can also be populated in complete-fusion reactions with stable projectiles. Experiments revealed so far that the yields which can be expected from complete fusion and MNT are comparable in this area. However, in MNT reactions a wide span of different nuclides can be populated at the same experimental setting due to their broad excitation functions which makes MNT reactions in this case more efficient.
There remains the question how to reach the predicted center of the "island of stability" at Z = 114, 120-126 and N = 184. Nuclei with such large neutron numbers are from present state of knowledge neither reachable in MNT reactions nor in complete-fusion reactions with RIBs due to the tiny cross-sections. But RIBs might hold an "emergency solution". The RIB intensities are large enough to study quasifission and fusion-fission reactions in respective systems because these reactions have large cross-sections. The hope is that possible shell closures leave their fingerprint not only in fusion-evaporation residues but also in the mass, angle and energy distributions of quasi-fission and fusion-fission fragments. The challenge is to distinguish these signatures in the measured distributions. We will start respective studies in near future at the new HIE-ISOLDE facility at CERN.
Synthesis of new nuclides is usually connected with experiments at the limits of feasibility. Therefore, reliable model predictions are extremely important for the successful conductance of the experiments. Roughly, the different models can be divided into macroscopic, macroscopic-microscopic, and purely microscopic. Among them, the presently most wide-spread models are the macroscopic-microscopic models which can again be divided in two main branches. One is the class of DNS models, which use the diffusion master equation, diabatic internuclear potentials and the quantum nature of interacting nuclei. Models of the other class use Langevin-type equations of motion and almost adiabatic internuclear potentials. These models indicate some common trends and act as good guidelines for experiments.

Ten supplementary questions
1. How many new transuranium and superheavy nuclei can we expect in the near and mid-term future?
We estimate that MNT reactions can provide roughly 150 new transuranium isotopes in the neutron-deficient and neutron-rich area. Another 50 new superheavy nuclei can still be expected from fusion reactions with stable beams, resulting in a total of about 200 new isotopes. This is a good number compared to the 348 presently known transuranium nuclides.
2. Which collision systems seem more favorable to synthesize neutron-rich superheavy nuclei in MNT reactions, intermediate heavy projectiles on heavy target nuclei like 48 Ca + 248 Cm, or heavy projectiles on heavy target nuclei like 238 U + 248 Cm?
Experimental results reveal that the choice of projectile plays no big role as long as the envisaged MNT products are not too far from the target nucleus (see Sect. 4.2.2). For MNT products with up to 5 protons above the target, the crosssections obtained with 48 Ca and 238 U projectiles on 248 Cm targets were the same within error bars. But one must also consider that lighter beams like 48 Ca are usually provided with much higher intensities than 238 U beams, which will lead to, respectively, larger yields at the same cross-section. 3. Which are the most serious bottlenecks if we plan to use MNT reactions to produce neutron-rich superheavy nuclei?
We regard the lack of appropriate separation and detection techniques as most serious restriction. The wide angular and energy distributions of MNT products make In-flight separation inefficient and for isotope identification the sensitive α decay tagging method largely fails in the neutron-rich superheavy region. We have to deal with cross-sections on the picobarn scale resulting in event rates of one nucleus per week. Therefore we cannot afford large losses, and the sensitivity of the experimental setup must approximately meet the production cross-section. The most sensitive techniques which are presently investigated as possible alternatives to α decay tagging are so far applicable for microbarn crosssections and we still need to bridge six orders of magnitude to make them suitable for picobarns. 4. Are MNT reactions really better than fragmentation to synthesize N = 126 nuclei?
At the moment, we regard fragmentation reactions still as unrivalled to produce neutron-rich nuclei below Pb. It is true that first experiments confirmed model predictions according to which MNT cross-sections are comparable to, or even exceeding fragmentation cross-sections. But in projectile fragmentation reactions the experimental conditions are unequally better which means that: (i) 10000 times thicker targets lead to 10,000 times higher yields for the same crosssection and beam intensity, (ii) strongly forward directed kinematics enables very efficient In-flight separation, (iii) relativistic energies allow for isotope identification by E-ΔE-TOF measurement, sensitive enough to detect single isotopes independent of their decay properties. MNT products in this region can so far only be identified via their gamma decays which means serious restrictions: the method is only applicable if excited states in the daughter nucleus are populated with sufficient branching and gamma transitions in the daughter nucleus must be already well known in order to attribute them correctly. Nevertheless, it is still very important to collect more data on MNT reactions in this region. Particularly the observed trend that MNT cross-sections increase strongly with respect to fragmentation cross-sections toward lower Z of the MNT products is very interesting. In reactions of Xe + Pb, MNT cross-sections started to exceed fragmentation cross-sections by four orders of magnitude for residues with N = 126, Z ≤ 77. So far, these MNT cross-sections were deduced indirectly from the ones of projectile-like MNT products and it is now very important to verify this effect by directly identifying the target-like MNT products. 5. Are MNT reactions with RIBs an option to produce neutron-rich isotopes?
It is revealed experimentally and theoretically that more neutrons in the system shift the isotopic distributions of MNT products toward the neutron-rich side. Therefore, MNT reactions with neutron-rich RIBs should lead to more neutronrich MNT products. However, one is quite restricted in the choice of suitable exotic projectile nuclei from terms of intensity. The best suitable RIBs are heavy alkali elements like Rb or Cs and noble gases like Xe. They are available with good intensity for a broad region of their exotic isotopes. But the intensities of stable and relatively neutron-rich beams like 48 Ca are still some orders of magnitude larger. Finally, it is the product of beam intensity and cross-section, I beam × σ , which determines if the yield of a certain isotope is larger with an exotic or with a stable projectile. 6. Would it be better to use inverse kinematics for MNT reactions?
Apart from the advantages of the more forward focused kinematics in such reactions, one has to consider some technical aspects. For example, very heavy beams are usually delivered with much smaller intensities than lighter beams which leads to smaller yields at a given cross-section. Also, if one uses the heavy reaction partner as beam, one is limited to uranium as heaviest possible projectile. This is, for example, unfavorable if one wants to synthesize new superheavy nuclei which are well above uranium. In earlier experiments it was found that the cross-section of heavy MNT products decreases on average by one order of magnitude with the transfer of every proton more. 7. Which still empty region on the nuclide chart is most difficult to access? This is the predicted island of stability in the superheavy element region. Fusion reactions with stable projectiles do not provide enough neutrons. Fusion reactions with RIBs fail due to the by far too small RIB intensities. Transfer reactions appear not feasible due to the tiny cross-sections. In near and mid-term future one can expect only indirect information about possible enhanced stability in this area, for example by studying quasi-fission and fusion-fission reactions with neutron-rich RIBs (see Sect. 5.3.2).

Are MNT reactions rather an emergency solution or a rich source for new isotopes?
It depends on the region which we intend to populate. To produce new neutron-rich superheavy nuclei, MNT reactions are presently the only alternative -and to our opinion an emergency solution. The kinematically much more favorable fusion reactions are not applicable here due to the lack of neutrons in the system if one uses stable projectiles and due to the tiny intensities if one uses neutron-rich RIBs. The weak point in low-energy MNT reactions is their kinematics. Large angular distributions and energy spread makes their separation and distinction from background events very difficult.
To produce neutron-rich N = 126 nuclei below Pb, MNT reactions might even be not necessary because the very effective fragmentation reactions seem still the better alternative (see question 4).
Finally, there is the region of neutron-deficient transuranium isotopes which can principally also be accessed in complete-fusion reactions. Experimental data reveal that the yields become increasingly comparable in fusion and MNT reactions the more one moves toward the neutron-deficient side. At that point, an advantage of MNT reactions comes into play, namely, that they are not selective on few distinct reaction products like complete fusion but allow for the population of a broad region of nuclei in the same experiment. Here, indeed MNT reactions might be the better source for new isotopes. 9. Which are the best experimental setups for nucleosynthesis in MNT reactions? Table 4 summarizes the experimental techniques which were so far applied in searches for new isotopes in MNT reactions. The lowest cross-sections were so far reached with the radiochemical and the velocity filter method. In some sense they are complementary to each other. With radiochemical methods MNT products emitted to large opening angles can be collected while the velocity filter is selective on products emitted to forward angles. Another difference is that radiochemical methods need more time and are better suitable for nuclei with minimum half-lives of about 1 min. The In-flight technique is suitable for very short-lived nuclides, where the limitation is mainly given by the flight time of the nuclei through the separator which is typically on the 1 µs scale. An important and noteworthy point is that both techniques reached their small cross-section limits only for isotopes which can be identified by α decay tagging! 10. Is it reasonable to consider MNT reactions at intermediate energies?
Binary reaction products have been experimentally observed among the products of projectile fragmentation reactions at incident energies around and above the Fermi energy. The yields of the most neutron-rich isotopes of light nuclei tend toward the Q gg systematics, indicating MNT processes. Therefore, the question arises of how important the MNT binary process is in peripheral collisions at quite high bombarding energies. Experimental and theoretical investigations of this subject are highly desirable.