Search

The dynamics of $\eta$ meson produced in proton-induced nuclear reactions via the decay of N$^{\ast}(1535)$ has been investigated within the Lanzhou quantum molecular dynamics transport model (LQMD). The in-medium modifications of the $\eta$ production in dense nuclear matter are included in the model, in which an attractive $\eta$-nucleon potential is implemented. The impact of the $\eta$ optical potential on the $\eta$ dynamics is investigated. It is found that the attractive potential leads to the reduction of high-momentum (kinetic energy) production from the spectra of momentum distributions and inclusive cross sections and increasing the reabsorption process by surrounding nucleons., Comment: 5 pages, 4 figures. arXiv admin note: text overlap with arXiv:1509.04792

Within the framework of the Lanzhou quantum molecular dynamics (LQMD) transport model, pion dynamics in heavy-ion collisions near threshold energies and the emission of preequilibrium particles (nucleons and light complex fragments) have been investigated. A density, momentum and isospin dependent pion-nucleon potential based on the $\Delta$-hole model is implemented in the transport approach, which slightly leads to the increase of the $\pi^{-}/\pi^{+}$ ratio, but reduces the total pion yields. It is found that a bump structure of the $\pi^{-}/\pi^{+}$ ratio in the kinetic energy spectra appears at the pion energy close to the $\Delta$(1232) resonance region. The yield ratios of neutrons to protons from the squeeze-out particles perpendicular to the reaction plane are sensitive to the stiffness of nuclear symmetry energy, in particular at the high-momentum (kinetic energy) tails., Comment: 8 pages, 9 figures, submitted EPJA. arXiv admin note: text overlap with arXiv:1509.04792

A statistical approach based on the Weisskopf evaporation theory has been developed to describe the de-excitation process of highly excited heavy and superheavy nuclei, in particular for the proton-rich nuclei. The excited nucleus is cooled by evaporating $\gamma$-ray, light particles (neutrons, protons, $\alpha$ etc) in competition with the binary fission, in which the structure effects (shell correction, fission barrier, particle separation energy) contribute to the processes. The formation of residual nuclei is evaluated via sequential emission of possible particles above the separation energies. The available data of fusion-evaporation excitation functions in the $^{28}$Si+$^{198}$Pt reaction can be reproduced nicely well within the approach., Comment: 10 pages, 4 figures

Transport simulations are very valuable for extracting physics information from heavy-ion collision experiments. With the emergence of many different transport codes in recent years, it becomes important to estimate their robustness in extracting physics information from experiments. We report on the results of a transport code comparison project. 18 commonly used transport codes were included in this comparison: 9 Boltzmann-Uehling-Uhlenbeck-type codes and 9 Quantum-Molecular-Dynamics-type codes. These codes have been required to simulate Au+Au collisions using the same physics input for mean fields and for in-medium nucleon-nucleon cross sections, as well as the same initialization set-up, the impact parameter, and other calculational parameters at 100 and 400 AMeV incident energy. Among the codes we compare one-body observables such as rapidity and transverse flow distributions. We also monitor non-observables such as the initialization of the internal states of colliding nuclei and their stability, the collision rates and the Pauli blocking. We find that not completely identical initializations constitute partly for different evolutions. Different strategies to determine the collision probabilities, and to enforce the Pauli blocking, also produce considerably different results. There is a substantial spread in the predictions for the observables, which is much smaller at the higher incident energy. We quantify the uncertainties in the collective flow resulting from the simulation alone as about $30\%$ at 100 AMeV and $13\%$ at 400 AMeV, respectively. We propose further steps within the code comparison project to test the different aspects of transport simulations in a box calculation of infinite nuclear matter. This should, in particular, improve the robustness of transport model predictions at lower incident energies where abundant amounts of data are available., Comment: 26 pages, 13 figures, a manuscript on the heavy-ion collision part of transport code comparison project

The neck dynamics in Fermi-energy heavy-ion collisions, to probe the nuclear symmetry energy in the domain of sub-saturation densities, is investigated within an isospin dependent transport model. The single and double ratios of neutron/proton from free nucleons and light clusters (complex particles) in the isotopic reactions are analyzed systematically. Isospin effects of particles produced from the neck fragmentations are explored, which are constrained within the midrapidities ($|y/y_{proj}|<$0.3) and azimuthal angles (70$^{o}\sim$110$^{o}$, 250$^{o}\sim$290$^{o}$) in semiperipheral nuclear collisions. It is found that the ratios of the energetic isospin particles strongly depend on the stiffness of nuclear symmetry energy and the effects increase with softening the symmetry energy, which would be a nice probe for extracting the symmetry energy below the normal density in experimentally. A flat structure appears at the tail spectra from the double ratio distributions. The neutron to proton ratio of light intermediate mass fragments (IMFs) with charged number Z$\leq$8 is related to the density dependence of symmetry energy with less sensitivity in comparison to the isospin ratios of nucleons and light particles., Comment: 6 pages, 6 figures

Formation mechanism of fragments with strangeness in collisions of antiprotons on nuclei has been investigated within the Lanzhou quantum molecular dynamics (LQMD) transport approach combined with a statistical model (GEMINI) for describing the decays of excited fragments. Production of strange particles in the antiproton induced nuclear reactions is modeled within the LQMD model, in which all possible reaction channels such as elastic scattering, annihilation, charge exchange and inelastic scattering in antibaryon-baryon, baryon-baryon and meson-baryon collisions have been included. A coalescence approach is developed for constructing hyperfragments in phase space after de-excitation of nucleonic fragments. The combined approach could describe the production of fragments in low-energy antiproton induced reactions. Hyperfragments are formed within the narrower rapidities and lower kinetic energies. It has advantage to produce heavier hyperfragments and hypernuclides with strangeness s=-2 (double-$\Lambda$ fragments) and s=1 ($\overline{\Lambda}$-fragments) in antiproton induced reactions., Comment: 6 pages, 5 figures

Dynamics of pseudoscalar mesons ($\pi$, $\eta$, $K$ and $\overline{K}$ ) and hyperons ($\Lambda$ and $\Sigma$) produced in heavy-ion collisions near threshold energies has been investigated within the Lanzhou quantum molecular dynamics (LQMD) transport model. The in-medium modifications on particle production in dense nuclear matter are implemented in the model through corrections on the elementary cross sections and by inclusion of the meson (hyperon)-nucleon potentials, in which the isospin effects are considered. It is found that the transportation of particles are influenced with the in-medium corrections. The total number of pions is reduced with an isospin, density and momentum dependent pion-nucleon potential. However, the ratios of charged pions is enhanced with inclusion of the potential. The production of eta in the domain of mid-rapidities and high momenta is sensitive to the $\eta$-nucleon potential, but weakly depends on symmetry energy. The attractive antikaon-nucleon potential enhances the subthreshold $\overline{K}$ production and also influences the structure of phase-space distributions. Dynamics of etas, kaons, antikaons and hyperons is also influenced by the pion potential because of collisions between pions and nucleons (resonances). The impacts of mean-field potentials on particle dynamics are investigated, such as the phase-space distributions from rapidity and transverse momentum spectra, inclusive invariant spectra, collective flows etc., Comment: 18 pages, 13 figures. arXiv admin note: text overlap with arXiv:1310.8072

Reaction dynamics in collisions of antiprotons on nuclei is investigated within the Lanzhou quantum molecular dynamics model. The reaction channels of elastic scattering, annihilation, charge exchange and inelastic collisions of antiprotons on nucleons have been included in the model. Dynamics on particle production, in particular pions, kaons, antikaons and hyperons, is investigated in collisions of $\overline{p}$ on $^{12}$C, $^{20}$Ne, $^{40}$Ca and $^{181}$Ta from a low to high incident momenta. It is found that the annihilations of $\overline{p}$ on nucleons are of importance on the dynamics of particle production in phase space. Hyperons are mainly produced via meson induced reactions on nucleons and strangeness exchange collisions, which lead to the delayed emission in antiproton-nucleus collisions., Comment: 8 pages, 4 figures. arXiv admin note: substantial text overlap with arXiv:1405.1997

The dynamics of nucleon transfer processes in heavy-ion collisions is investigated within the dinuclear system model. The production cross sections of nuclei in the reactions $^{136}$Xe+$^{208}$Pb and $^{238}$U+$^{248}$Cm are calculated and the calculations are in good agreement with the experimental data. The transfer cross sections for the $^{58}$Ni+$^{208}$Pb reaction are calculated and compared with the experimental data. We predict the production cross sections of neutron-rich nuclei $^{165-168}$Eu, $^{169-173}$Tb, $^{173-178}$Ho, and $^{181-185}$Yb based on the reaction $^{176}$Yb+$^{238}$U. It can be seen that the production cross sections of the neutron-rich nuclei $^{165}$Eu, $^{169}$Tb, $^{173}$Ho, and $^{181}$Yb are 2.84 $\mu$b, 6.90 $\mu$b, 46.24 $\mu$b, and 53.61 $\mu$b, respectively, which could be synthesized in experiment.

Dynamics of strange particles produced in proton induced nuclear reactions near threshold energies has been investigated within the Lanzhou quantum molecular dynamics (LQMD) transport model. The in-medium modifications on particle production in dense nuclear matter are considered through corrections on the elementary cross sections via the effective mass and the mean-field potentials. It is found that the attractive antikaon-nucleon potential enhances the subthreshold $\bar{K}$ production and also influences the structure of inclusive spectra. The strangeness production is strongly suppressed in proton induced reactions in comparison to heavy-ion collisions. The measured K$^{-}$/K$^{+}$ ratios in the $^{12}$C+$^{12}$C and p+$^{12}$C collisions from KaoS collaboration have been well explained with inclusion of the in-medium potentials., Comment: 7 pages, 5 figures

The quantum molecular dynamics model has been improved to investigate the reaction dynamics induced by antiprotons. The reaction channels of elastic scattering, annihilation, charge exchange and inelastic collisions have been included in the model. Dynamics on particle production, in particular pions, kaons, antikaons and hyperons, is investigated in collisions of $\overline{p}$ on $^{12}$C, $^{20}$Ne, $^{40}$Ca, $^{112}$Sn, $^{181}$Ta, $^{197}$Au and $^{238}$U from a low to high incident momentum. The rapidity and momentum distributions of $\pi^{+}$ and protons from the LEAR measurements can be well reproduced. The impacts of system size and incident momentum on particle emissions are investigated from the inclusive spectra, transverse momentum and rapidity distributions. It is found that the annihilations of $\overline{p}$ on nucleons are of importance on the particle production. Hyperons are mainly produced via meson induced reactions on nucleons and strangeness exchange collisions when the incident momentum is below the threshold value of annihilation reaction. A higher nuclear temperature is obtained from the kaon emission, but it has a lower value for hyperon production., Comment: 11 pages, 8 figures

Strangeness (K$^{0,+}$, $\Lambda$ and $\Sigma^{-,0,+}$) production in heavy-ion collisions near threshold energies has been investigated within the Lanzhou quantum molecular dynamics (LQMD) transport model. The kaon (anti-kaon)-nucleon and hyperon-nucleon potentials in dense nuclear matter are implemented in the model, which are computed from the chiral perturbation approach and the relativistic mean-field model, respectively. It is found that the in-medium potentials change the structure of transverse flow, and also affect the rapidity distributions and the inclusive spectra for strangeness production. The local temperature of the fire ball extracted from the kaon spectra of inclusive invariant cross sections is influenced by the kaon-nucleon potential. The stiffness of nuclear symmetry energy and the kaon-nucleon potential by distinguishing isospin effect play a significant role on the ratio of K$^{0}$/K$^{+}$, in particular at the subthreshold energies. The ratio of $\Sigma^{-}/\Sigma^{+}$ depends on the high-density symmetry energy, in which the $\Sigma$-nucleon potential has a neglectable contribution on the isospin ratio., Comment: 25 pages, 13 figures. arXiv admin note: text overlap with arXiv:1304.2897

Subthreshold kaon (K$^{0}$ and K$^{+}$) production in neutron-rich nuclear reactions to probe the kaon-nucleon interaction in nuclear medium and to extract the isospin dependent part of the nuclear equation of state at high-baryon densities, is investigated within an isospin and momentum dependent transport model. A repulsive kaon-nucleon potential is implemented in the model through fitting the flow data and inclusive spectra in heavy-ion collisions, which enhances the energetic kaon emission squeezed out in the reaction zone and leads to a variation of the high-momentum spectrum of the K$^{0}$/K$^{+}$ yields. It is found that the stiffness of nuclear symmetry energy plays a significant role on the isospin ratio with decreasing the incident energy and a hard symmetry energy has a larger value of the K$^{0}$/K$^{+}$ ratio in the domain of subthreshold energies., Comment: 8 pages, 6 figures. arXiv admin note: text overlap with arXiv:1111.6863

A recent updated version of the Lanzhou quantum molecular dynamics is reviewed, in which the momentum dependence of symmetry potential and the effective mass splitting of proton and neutron in nuclear medium are included in the model. The in-medium nucleon-nucleon (NN) elastic cross sections are evaluated with the scaling approach according to the effective mass and their influence on collective flows in heavy-ion collisions are discussed. The inelastic cross sections, in particular for the process of the channel N$\Delta\rightarrow$NN, are parameterized in accordance with the available experimental data. It is found that the in-medium cross sections play a significant role on isospin emissions and result in a flat distribution for transverse flows and elliptic flows of free nucleons comparing with the in-vacuum ones. The rapidity distribution of the difference between neutron and proton transverse flows is sensitive to the stiffness of nuclear symmetry energy as a promising observable, which can not be influenced by the in-medium effect and collision centrality. Furthermore, the elliptic flow of free energetic nucleons in the case of the mass splitting of $m_{n}^{\ast}>m_{p}^{\ast}$ is also related to the symmetry energy. However, the pion flows weakly depend on the symmetry energy and the mass splitting., Comment: 14 pages, 19 figures. arXiv admin note: substantial text overlap with arXiv:1106.3797

Transverse emission of preequilibrium nucleons, light clusters (complex particles) and charged pions from the isotopic $^{112,124}$Sn+$^{112,124}$Sn reactions at a beam energy of 400\emph{A} MeV, to extract the high-density behavior of nuclear symmetry energy, are investigated within an isospin and momentum dependent transport model. Specifically, the double ratios of neutron/proton, triton/helium-3 and $\pi^{-}/\pi^{+}$ in the squeeze-out domain are analyzed systematically, which have the advantage of reducing the influence of the Coulomb force and less systematic errors. It is found that the transverse momentum distribution of isospin ratios strongly depend on the stiffness of nuclear symmetry energy, which would be a nice observable to extract the high-density symmetry energy. The collision centrality and the mass splitting of neutron and proton in nuclear medium play a significant role on the distribution structure of the ratios, but does not change the influence of symmetry energy on the spectrum., Comment: 5 figures, 13 pages

Within the framework of an isospin and momentum dependent transport model, the emissions of isospin particles (nucleons and light clusters) squeezed out in heavy-ion collisions are investigated as probes of the poorly known symmetry energy at high baryon density. Two different mass splittings of neutrons and protons in nuclear medium as $m_{n}^{\ast}>m_{p}^{\ast}$ and $m_{n}^{\ast}

A Skyrme-type momentum-dependent nucleon-nucleon force distinguishing isospin effect is parameterized and further implemented in the Lanzhou Quantum Molecular Dynamics (LQMD) model for the first time, which leads to a splitting of nucleon effective mass in nuclear matter. Based on the isospin- and momentum-dependent transport model, we investigate the influence of momentum-dependent symmetry potential on several isospin-sensitive observables in heavy-ion collisions. It is found that symmetry potentials with and without the momentum dependence but corresponding to the same density dependence of the symmetry energy result in different distributions of the observables. The mid-rapidity neutron/proton ratios at high transverse momenta and the excitation functions of the total $\pi^{-}/\pi^{+}$ and $K^{0}/K^{+}$ yields are particularly sensitive to the momentum dependence of the symmetry potential., Comment: 12 pages, 5 figures

We present a review of the recent progress of theoretical models on the description of the formation of superheavy nuclei in collisions of heavy systems. Two sorts of reactions that are the fusion-evaporation mechanism and the massive damped collisions to produce superheavy nuclei are discussed. Problems and further improvements of the capture of colliding partners, the formation of compound nucleus and the de-excitation process are pointed out. Possible combinations in the synthesis of the gap of the cold fusion and $^{48}$Ca induced reactions are proposed by the calculations based on the dinuclear system model and also compared with other approaches. The synthesis of neutron-rich heavy isotopes near sub-shell closure N=162 via transfer reactions in the damped collisions of two actinides and the influence of shell closure on the production of heavy isotopes are investigated. Prospective possibility to reach superheavy nuclei near N=184 via neutron-rich radioactive beams of high intensity in the future is discussed., Comment: 22 pages, 17 figures

The dynamics of pions and strange particles in heavy-ion collisions in the region of 1A GeV energies is investigated by the lanzhou quantum molecular dynamics model for probing the nuclear equation of state at supra-saturation densities. The total multiplicities and the ratios obtained in $^{197}$Au+$^{197}$Au over $^{12}$C+$^{12}$C systems are calculated for selected Skyrme parameters SkP, SLy6, Ska and SIII, which correspond to different modulus of incompressibility of symmetric nuclear matter and different cases of the stiffness of symmetry energy. A decrease trend of the excitation functions of the ratios for strange particle production with increasing incident energy was observed. The available data of K$^{+}$ production measured by KaoS collaboration are described well with the parameter SkP, which results in a soft equation of state. The conclusions can not be modified by an in-medium kaon-nucleon potential., Comment: 7 pages, 6 figures

Within the framework of the improved isospin dependent quantum molecular dynamics (ImIQMD) model, the dynamics of strangeness (K$^{0,+}$, $\Lambda$ and $\Sigma^{-,0,+}$) production in heavy-ion collisions near threshold energies is investigated systematically, in which the strange particles are considered to be mainly produced by the inelastic collisions of baryon-baryon and pion-baryon. The collisions in the region of supra-saturation densities of the dense baryonic matter formed in heavy-ion collisions dominate the yields of strangeness production. The total multiplicities as functions of incident energies and collision centralities are calculated with the Skyrme parameter SLy6. The excitation function of strangeness production is analyzed and also compared with the KaoS data for the K$^{+}$ production in the reactions $^{12}$C+$^{12}$C and $^{197}$Au+$^{197}$Au., Comment: 4 pages, 4 figures

Within the framework of the improved isospin dependent quantum molecular dynamics (ImIQMD) model, the dynamics of pion emission in heavy-ion collisions in the region 1 A GeV as a probe of nuclear symmetry energy at supra-saturation densities is investigated systematically. The total pion multiplicities and the $\pi^{-}/\pi^{+}$ yields are calculated for selected Skyrme parameters SkP, SLy6, Ska and SIII, and also for the cases of different stiffness of symmetry energy with the parameter SLy6. The influence of Coulomb potential, symmetry energy and in-medium pion potential on the pion production is investigated and compared each other by analyzing the distributions of transverse momentum and longitudinal rapidity and also the excitation functions of the total pion and the $\pi^{-}/\pi^{+}$ ratio. The directed flow, elliptic flow and polar angle distributions are calculated at the cases of different collision centralities and symmetry energies., Comment: 12 pages, 14 figures

The dynamics of transfer reactions in collisions of two very heavy nuclei $^{238}$U+$^{238}$U is studied within the dinuclear system (DNS) model. Collisions of two actinide nuclei form a super heavy composite system during a very short time, in which a large number of charge and mass transfers may take place. Such reactions have been investigated experimentally as an alternative way for the production of heavy and superheavy nuclei. The role of collision orientation in the production cross sections of heavy nuclides is analyzed systematically. Calculations show that the cross sections decrease drastically with increasing the charged numbers of heavy fragments. The transfer mechanism is favorable to synthesize heavy neutron-rich isotopes, such as nuclei around the subclosure at N=162 from No (Z=102) to Db (Z=105)., Comment: 4 pages, 4 figures

Within the framework of the dinuclear system (DNS) model, production cross sections of new superheavy nuclei with charged numbers Z=108-114 are analyzed systematically. Possible combinations based on the actinide nuclides $^{238}$U, $^{244}$Pu and $^{248,250}$Cm with the optimal excitation energies and evaporation channels are pointed out to synthesize new isotopes which lie between the nuclides produced in the cold fusion and the $^{48}$Ca induced fusion reactions experimentally, which are feasible to be constructed experimentally. It is found that the production cross sections of superheavy nuclei decrease drastically with the charged numbers of compound nuclei. Larger mass asymmetries of the entrance channels enhance the cross sections in 2n-5n channels., Comment: 5 pages, 4 figures

The enhancement of the sub-barrier fusion cross sections was explained as the lowering of the dynamical fusion barriers within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model. The numbers of nucleon transfer in the neck region are appreciably dependent on the incident energies, but strongly on the reaction systems. A comparison of the neck dynamics is performed for the symmetric reactions $^{58}$Ni+$^{58}$Ni and $^{64}$Ni+$^{64}$Ni at energies in the vicinity of the Coulomb barrier. An increase of the ratios of neutron to proton in the neck region at initial collision stage is observed and obvious for neutron-rich systems, which can reduce the interaction potential of two colliding nuclei. The distribution of the dynamical fusion barriers and the fusion excitation functions are calculated and compared them with the available experimental data., Comment: 9 pages, 5 figures

Within the framework of the improved isospin dependent quantum molecular dynamics (ImIQMD) model, the pion emission in heavy-ion collisions in the region 1 A GeV is investigated systematically, in which the pion is considered to be mainly produced by the decay of resonances $\triangle$(1232) and N*(1440). The in-medium dependence and Coulomb effects of the pion production are included in the calculation. Total pion multiplicity and $\pi^{-}/\pi^{+}$ yields are calculated for the reaction $^{197}$Au+$^{197}$Au in central collisions for selected Skyrme parameters SkP, SLy6, Ska, SIII and compared them with the measured data by the FOPI collaboration., Comment: 9 pages, 5 figures

Within the framework of the improved isospin dependent quantum molecular dynamics (ImIQMD) model, the emission of pion in heavy-ion collisions in the region 1 A GeV as a probe of nuclear symmetry energy at supra-saturation densities is investigated systematically, in which the pion is considered to be mainly produced by the decay of resonances $\triangle$(1232) and N*(1440). The total pion multiplicities and the $\pi^{-}/\pi^{+}$ yields are calculated for selected Skyrme parameters SkP, SLy6, Ska and SIII, and also for the cases of different stiffness of symmetry energy with the parameter SLy6. Preliminary results compared with the measured data by the FOPI collaboration favor a hard symmetry energy of the potential term proportional to $(\rho/\rho_{0})^{\gamma_{s}}$ with $\gamma_{s}=2$., Comment: 12 pages, 5 figures

Within the framework of the dinuclear system (DNS) model, the production cross sections of superheavy nuclei Hs (Z=108) and Z=112 combined with different reaction systems are analyzed systematically. It is found that the mass asymmetries and the reaction Q values of the combinations play a very important role on the formation cross sections of the evaporation residues. Both methods by solving the master equations along the mass asymmetry degree of freedom (1D) and along the proton and the neutron degrees of freedom (2D) are compared each other and with the available experimental results., Comment: 10 pages, 5 figures

Within the framework of a dinuclear system (DNS) model, the evaporation-residue excitation functions and the quasi-fission mass yields in the $^{48}$Ca induced fusion reactions are investigated systematically and compared with available experimental data. Maximal production cross sections of superheavy nuclei based on stable actinide targets are obtained. Isotopic trends in the production of the superheavy elements Z=110, 112-118 based on the actinide isotopic targets are analyzed systematically. Optimal evaporation channels and combinations as well as the corresponding excitation energies are proposed. The possible factors that influencing the isotopic dependence of the production cross sections are analyzed. The formation of the superheavy nuclei based on the isotopes U with different projectiles are also investigated and calculated., Comment: 26 pages, 15 figures

The shell correction is proposed in the improved isospin dependent quantum molecular dynamics (ImIQMD) model, which plays an important role in heavy-ion fusion reactions near Coulomb barrier. By using the ImIQMD model, the static and dynamical fusion barriers, dynamical barrier distribution in the fusion reactions are analyzed systematically. The fusion and capture excitation functions for a series of reaction systems are calculated and compared with experimental data. It is found that the fusion cross sections for neutron-rich systems increase obviously, and the strong shell effects of two colliding nuclei result in a decrease of the fusion cross sections at the sub-barrier energies. The lowering of the dynamical fusion barriers favors the enhancement of the sub-barrier fusion cross sections, which is related to the nucleon transfer and the neck formation in the fusion reactions., Comment: 20 pages, 12 figures

Within the concept of the dinuclear system (DNS), a dynamical model is proposed for describing the formation of superheavy nuclei in complete fusion reactions by incorporating the coupling of the relative motion to the nucleon transfer process. The capture of two heavy colliding nuclei, the formation of the compound nucleus and the de-excitation process are calculated by using an empirical coupled channel model, solving a master equation numerically and applying statistical theory, respectively. Evaporation residue excitation functions in cold fusion reactions are investigated systematically and compared with available experimental data. Maximal production cross sections of superheavy nuclei in cold fusion reactions with stable neutron-rich projectiles are obtained. Isotopic trends in the production of the superheavy elements Z=110, 112, 114, 116, 118 and 120 are analyzed systematically. Optimal combinations and the corresponding excitation energies are proposed., Comment: 18 pages, 8 figures