Study on deuteron formation mechanism in nucleon-induced reactions

Ya-Jun He · Chen-Chen Guo· Jun Su · Long Zhu· Zhen-Dong An

Abstract The mechanism of deuteron formation in neutron-induced reactions is studied within the framework of the isospin-dependent quantum molecular dynamics model,using the GEMINI code. The influence of the n + p →d reaction channel is investigated by analyzing the deuteron production cross sections in the neutron-induced reactions 12C(n,d), 16O(n,d), and 28Si(n,d), with incident energies of 20-100 MeV. By including the n + p →d reaction channel when modeling the collision, the deuteron production cross sections increase, optimizing the cross-section results and bringing them closer to the experimental data values. This indicates that the n + p →d reaction channel is an important mechanism for enhancing deuteron production.

Keywords Cluster mechanism · Deuteron formation cross section · Nucleon-induced reactions

1 Introduction

Spallation reactions, where high-energy, light particles collide with a heavy target, resulting in a lighter remnant nucleus and the ejection of numerous light particles, play an important role in a wide range of applications [1]. Ever since cosmic ray spallation reactions were studied in the 1930s [2], they have attracted much attention from researchers in applied and fundamental fields [3]. This interest only increased with the proposition of an optimum neutron resource, as well as the application of nuclear waste transmutation in an accelerator-driven system. Over the years, there have been extensive developments of applications entailing spallation reactions, including material physics [4], nuclear waste disposal [5-7], particle physics and nuclear physics [8], rare isotope production,accelerator radiation protection [9], cancer hadron therapy[11], and cosmic rays in the atmosphere [10].

A nucleon-induced spallation reaction is a two-stage process. First, the hard nucleon-nucleon (NN) collisions take place, dissipating the incident energy and resulting in a nuclei with high excitation energy. In the second stage,de-excitation takes place via evaporation and fission processes[12,13].In addition to fission and fusion,spallation is an effective process of releasing neutrons from nuclei;unlike the other two reactions,it is an endothermic process.To sustain a spallation reaction, a beam of high-energy particles,usually protons,must be supplied and focused on a heavy target. Spallation can be an important source of neutrons whose flux can easily be controlled by the drive beam [14]. However,it is hard to obtain the large quantity of reaction data necessary, in consideration of the low feasibility and high costs.Therefore,theoretical studies are imperative to improve our understanding of these reactions.There are many spallation models and codes currently used, including the particle and heavy-ion transport code system (PHITS) [15], the Boltzmann-Uehling-Uhlenbeck(BUU) transport model [16], the quantum molecular dynamics (QMD) model [17-19], the cascade model(JAM1.0) [20], the intranuclear cascade model (INCL)4[21], ISABEL [22, 23], and the Monte Carlo code GEANT4 [24, 25]. These models and codes are used in conjunction with various statistical decay models, such as the generalized evaporation model (GEM) [26, 27],GEMINI [28], and the abrasion-ablation model (ABLA)[29]. More recently, a Bayesian neural network (BNN)method has been proposed for spallation reactions [30].Nonetheless, the discrepancies between observables in theoretical calculations and experimental data are still relatively large,proving that the theoretical models are not well founded. One of the problems the theoretical models struggle with is being able to reproduce the production cross sections of the light cluster well.Attempts have been made to study the formation mechanism of clusters,such as deuteron, to improve upon the modeling. A surface coalescence mechanism has been used in INCL, Jaeri QMD(JQMD), and improved QMD (ImQMD) models to improve the description of the light complex particles emitted in various nucleon-induced reactions [21, 31, 32].The pre-equilibrium light complex particle emission has been described reasonably well [18, 19]. An investigation into the formation mechanism of deuteron clusters is conducive to solving the current controversial clustering problem.

In this work, we study neutron-induced reactions by using the isospin-dependent QMD model (IQMD) model.The GEMINI model is applied to simulate the pre-fragment decays. The version of the IQMD model used is IQMDBNU.The framework of the IQMD+GEMINI model has been applied to study nuclear reactions in an energy region from 10 MeV to GeV [33-42]. By using these models, we attempt to study the formation mechanism of deuteron clusters.The reaction channel of n + p →d is added into the model. In Sect. 2, a brief introduction of the IQMD model and cluster mechanism is made.Section 3 describes and analyzes the subsequent calculation results, with some brief discussion about the results.In Sect. 4,we summarize our findings, provide some concluding remarks, and discuss potential future directions that our work can take.

2 Model

2.1 IQMD-BNU

The IQMD model is based on the same principles as the QMD model [43-45]. By considering the evolution of nucleons in mean field potentials and the binary NN collisions,the IQMD model has been widely and successfully used for the analysis of heavy-ion collisions (HICs) at incident energies below 2 GeV/nucleon.

The IQMD-BNU model was introduced and compared with other model versions in a transport-code-comparison project [46, 47]. In a QMD-type model [48], each nucleon is represented by a Gaussian wave packet [49]

where φkiis the wavefunction of the ith particle at a state ki(pi=¯hki). φkiis chosen as the Gaussian wave packet to avoid the negative values of the phase-space distribution(fN). By the Wigner transformation of the N-body wavefunction, the phase-space density is given by

where ri0= ri〈 〉 and pi0= pi〈 〉 are the centroids of wave packets in coordinate and momentum space, respectively.The time evolution of the nucleons in the system is considered under self-consistent conditions.The mean field is governed by the Hamiltonian equations of motion,

The density can be calculated by taking the integral of the phase-space density [Eq. (4)] over momentum space. The nuclear potential Vn includes a two- and three-body Skyrme interaction term,as well as a symmetry potential term.The nuclear potential energy density of the asymmetric nuclear matter, with density ρ and asymmetry δ, is given by

where ρ0=0.16 fm-3is the normal density. The constant parameters are taken to have values of α=-356 MeV,β=303 MeV,γ=7/6，Csp =38.06 MeV,and γi=0.75.These values provide a compressibility of 200 MeV at saturation density for isospin symmetric nuclear matter.

The binary NN collisions are also included in our model.Specifically, the NN elastic scattering cross section can be expressed by the following formula,

2.2 GEMINI

The statistical code GEMINI,proposed by Charity in the 1980s [53], can be utilized to model the decay of a compound nucleus in fusion reactions and excited fragments in HICs. The GEMINI code is applied to the process of deexcitation of hot fragments [54, 55].

The IQMD code will stop running whenever the excitation energies of the two heaviest pre-fragments are less than a specified value, Estop =3 MeV/A. The charge number, mass number, and excitation energy of each prefragment are output by the IQMD code. If the excitation energy is greater than zero, the modeling of the sequential decay of the pre-fragment will be performed by the statistical code GEMINI.This statistical code not only allows light-particle evaporation and symmetric fission, but all possible binary decay modes.Subsequently,a Monte Carlo technique is used to track all the decay chains until the resulting products are unable to decay further. Part of the necessary decay width data comes from the Hauser-Feshbach formalism of the evaporation of light particles.

In summary, if the excited energy of the pre-fragments is lower than 3 MeV/A, the statistical decay model GEMINI is preferred. We start the simulation of GEMINI after a simulation time of 1000 fm/c in the IQMD model,also using the PSDC method.We assume that the evolution is stable enough until 1000 fm/c. For further details on GEMINI, refer to Ref. [53].

2.3 IQMD + GEMINI model updates

where σd is the cross section of the n + p →d reaction,and fmed d is the in-medium factor.The cross section of the n + p →d reaction has been measured in various experiments,with the incident energy ranging from 0.01 to 100 MeV [56-65]. The experimental data are depicted as circles in Fig.1.By fitting the data,an empirical formula of σd as a function of the incident energy is obtained,where ρ is the density in the center of the colliding neutron and proton, and ρ0is the normal density. This in-medium factor value means that the deuteron can only be produced in the low-density region, i.e., the deuteron is generated when a neutron and a proton are close enough, and the medium density is less than 0.2ρ0. When a deuteron evolves to a medium density greater than 0.2ρ0, it will break up into a neutron and a proton.

It should be noted that the gamma rays were ignored in the n+p →d+γ decay channel in this work, since the beam energy is much larger than the energy. As a weak composite system, the propagation of the deuterons is approximately treated the same as the propagation of a neutron and proton under a nucleonic mean field potential.For the deuteron break-up mechanism, we adopt a phenomenological method that causes the deuteron to decay when its medium density is less than 0.2ρ0. A long-term deuteron production mechanism in the QMD model is still an open question, which requires further study.

3 Results and discussion

The excited energy and charge distribution of pre-fragments for the n+28Si reaction at 60 MeV are displayed in Fig.2.The percentage of pre-fragments is represented by p.The excitation energy of the fragments can be calculated using,

To study the formation of a deuteron cluster in detail,the formation time,collision energy, density at its position of formation, and excitation energy are shown in Fig. 3.The highest percentage of deuteron clusters are formed in the time frame 0-100 fm/c (Fig. 3a). The deuteron cluster is formed at this time because the incoming neutron captures a proton in the target,in the low-density region.From 100 to 1000 fm/c, the deuteron is formed from the excited nucleus in the low-density area during the reaction. It should be noted that the neutrons that form deuterons are not necessarily incident neutrons. Figure 3b shows the relative kinetic energy of the neutron and proton in the excitation system when the deuteron cluster is formed. As the relative kinetic energy increases from 0 to 10 MeV,the probability of deuteron generation increases gradually.After 10 MeV, this probability gets smaller. At 60 MeV,the probability that a deuteron is formed is lower than one percent. This is because the cross section of deuterons decreased gradually above 10 MeV, as shown in Fig. 1, in turn causing the probability of deuteron production to also decrease gradually. Figure 3c shows the ratio between density and normal density, at the position where the deuteron cluster formed. The deuteron cluster prefers to form in a density ratio of 0.3-0.5.When the ratio is greater than 0.6, even though the deuteron cluster will form, the cluster will soon disappear owing to interactions with other nucleons. When the ratio is below 0.3, the density is too small for a deuteron cluster to form.Figure 3d displays the distribution of the total energy of the deuteron cluster.The reason for negative total energies is a result of the choice to ignore the gamma photons.

Figure 4 shows the deuteron cluster formation cross section in28Si(n,d),16O(n,d), and12C(n,d) reactions at beam energies of 20-100 MeV.The dotted line is the result of the deuteron cluster produced by fragmentation without the introduction of n + p →d reaction channel.The solid line is the result of the deuteron cluster produced by fragmentation and the neutron-proton collision, after the addition of this reaction channel. The hollow circles represent experimental data [66-68]. With the introduction of the n + p →d reaction channel into the modeling, the cross-section results are significantly optimized for each set of reactions, bringing the results much closer to the experimental data values. In Fig. 4a, the prediction of the cross section for28Si(n,d) reactions has the same trend asthe experimental data at 20-50 MeV. By calculating the ratio of neutron-proton collisions forming a deuteron cluster, it is found that the contribution rate is greater than 52.6%. By the same method, the contribution rate of neutron-proton collisions for deuteron cluster formation is 36.5% for16O(n,d) reactions (Fig. 4b) and is 47.3% for12C(n,d)reactions(Fig.4c).Hence,the addition of the n +p →d reaction channel is very significant when modeling the deuteron cluster cross section. It is necessary to place greater focus on neutron-proton collisions when studying the effect on the formation and evolution of deuteron clusters.

4 Summary

A combination of the IQMD model and the GEMINI code was applied to study neutron-induced reactions.It was shown that the deuteron production cross sections calculated by the standard version of the IQMD+GEMINI model are much smaller than the data. In an attempt to optimize the modeling of the deuteron production cross sections,the n + p →d reaction channel was added to the IQMD model. The empirical formula of the cross section and its in-medium factor of the n + p →d reaction were used. When the deuteron evolved to a medium density greater than 0.2ρ0, it decayed into a neutron and a proton by the reaction d+N →n+p+N,where N is a nucleon.The contribution of neutron-proton collisions forming deuteron clusters was calculated. It was found that the contribution rate was greater than 52.6% for28Si(n,d)reactions, 36.5% for16O(n,d) reactions, and 47.3% for12C(n,d) reactions. The deuteron production cross sections in the neutron-induced reactions12C(n,d),16O(n,d), and28Si(n,d) with incident energy 20-100 MeV were calculated and compared to the experimental data.The deuteron production cross sections increased having incorporated the n + p →d reaction channel in the model and bring the values closer to the experimental data. The results indicated that the deuteron cluster is generated in a low-density region via the neutron-proton collision.In future work,the deuteron production mechanism in the QMD still requires further study, but our results hint that the n + p →d reaction channel should be considered in the transport model when studying deuteron clusters in HICs.