Pre-neutron fragment mass yields for 235U(n, f ) and 239Pu(n, f ) reactions at incident energies from thermal up to 20 MeV

Figures(8) / Tables(3)

Get Citation
Fanglei Zou, Xiaojun Sun, Kai Zhang, Hongfei Chen, Jie Yan, Junlong Tian and Yunyi Cui. Pre-neutron fragment mass yields for 235U(n, f) and 239Pu(n, f) reactions at incident energies from thermal up to 20 MeV[J]. Chinese Physics C. doi: 10.1088/1674-1137/acb910
Fanglei Zou, Xiaojun Sun, Kai Zhang, Hongfei Chen, Jie Yan, Junlong Tian and Yunyi Cui. Pre-neutron fragment mass yields for 235U(n, f) and 239Pu(n, f) reactions at incident energies from thermal up to 20 MeV[J]. Chinese Physics C.  doi: 10.1088/1674-1137/acb910 shu
Milestone
Received: 2022-02-19
Article Metric

Article Views(1380)
PDF Downloads(32)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Pre-neutron fragment mass yields for 235U(n, f ) and 239Pu(n, f ) reactions at incident energies from thermal up to 20 MeV

  • 1. College of Physics, Guangxi Normal University, Guilin 541004, China
  • 2. Guangxi Key Laboratory of Nuclear Physics and Nuclear Technology, Guilin 541004, China
  • 3. Institute of Nuclear physics and chemistry, Mianyang 421000, China

Abstract: Pre-neutron fragment mass yields in the vicinity of the thermal neutron energy are highly important for applications because of the larger fission cross sections of the 235U(n, f) and 239Pu(n, f) reactions. In this paper, pre-neutron fragment mass yields at incident energies from thermal up to 20 MeV are systematically studied using an empirical fission potential (EFP) model, the potential parameters of which are obtained from the measured data. The energy dependences of the peaks and valleys of the pre-neutron fragment mass yields are described by exponential and linear functions for the 235U(n, f) and 239Pu(n, f) reactions, respectively. The energy dependences of the evaporation neutrons, which play a crucial role in the reasonable description of pre-neutron fragment mass yields, are also obtained from the fission cross sections. The pre-neutron fragment mass yields in this study are not only consistent with the results of previous studies in regions of several MeVs but also agree well with existing measured data at incident energies from thermal up to 20 MeV. The results show that the feasibility of this EFP model is verified in this extended energy region.

    HTML

    I.   INTRODUCTION
    • Since the discovery of nuclear fission in the late 1930s [1], finding a method of accurately describing the nuclear fission process has become one of the biggest challenges in theoretical research. The fission process not only involves the collective motion of a large number of nucleons but is also affected by various structural effects. Moreover, the complexity and richness of fission have been revealed theoretically and experimentally over more than 80 years [24]. Nuclear fission data have many practical applications, including nuclear safety measures, accelerator technology, homeland security, medicine, nuclear energy, nuclear reactor waste transmutation [57], and the r-process of neutron star mergers [8]. A precise description of the pre-neutron fragment mass yield, which is one of the most important quantities for neutron-induced fission, is of great importance in understanding the fission mechanism and practical applications [9].

      Several important theories and models have been developed to understand the fission mechanism or quantitatively calculate fragment mass yields [10, 11]. Calculating fission observables based on the time-dependent generating coordinate method and Hartree-Fock-Bogoliubov (HFB) microscopic model is difficult and time-consuming when considering the effective interaction of the finite force range [12]. Although the macro-dynamic model has been used to calculate the mass distribution of fission fragments, no accurate calculation methods could be obtained for several parameters that are sensitive to the viscosity coefficient and mass tensor of the fission process [13].

      The GEF model, a semi-empirical phenomenological model, combines physical concepts from quantum mechanics and statistical mechanics with specific experimental information to adjust a set of suitable parameters for different fission systems. It can avoid the limitation of microscopic calculation, simplify the calculation process, and give reliable predictions for a variety of fission observables for a large number of fissile nuclei [14].

      The pre-neutron fragment mass yields of actinide nuclei induced by neutrons in a region of several MeVs have been theoretically studied based on the empirical fission potential (EFP) model proposed in our previous study [15]. Immediately, the EFP model was extended up to 60 MeV for the 232Th(n, f) and 238U(n, f) reactions to quantitatively describe the pre-neutron fragment mass yields [9, 16]. However, there are still no publications describing pre-neutron fragment mass yields in the vicinity of thermal neutron energy. It is worth mentioning that several quantities after scission are highly important for applications because of the larger cross sections in the low incident energy region. Thus, the feasibility of this EFP model at the incident energies from thermal up to 20 MeV should be verified. In particular, the neutron-induced 235U and 239Pu (as important fuel nuclides) reactions have no threshold. Unfortunately, there is little experimental data on pre-neutron fragment mass yields at incident energies from thermal up to 20 MeV for these two reactions. Therefore, pre-neutron fragment mass yields should be theoretically recommended at incident energies from thermal up to the 20 MeV region.

      This paper is organized as follows: In Sec. II, the EFP model is briefly introduced. The calculated results of the pre-neutron fragment mass yields for the 232Th(n, f) and 238U(n, f) reactions are shown and compared with the results of the GEF model in Sec. III. Finally, a brief summary is given in Sec. IV.

    II.   THEORETICAL MODEL
    • Figure 8.  (color online) Predicted pre-neutron fragment mass yields for the 239Pu(n, f) reaction at an incident energy of 14 MeV.

      Based on the bimodal characteristic of actinide nuclear fission, the EFP model is proposed to describe the pre-neutron fragment mass yield [9, 15, 16]. In terms of the EFP model, the pre-neutron fragment mass yield can be expressed as

      Y(A)=CeU(A),

      (1)

      where C is the normalization constant, and the variable A denotes the mass number of the primary fragment.

      The empirical fission potential U(A) is expressed as

      U(A)={u1(AA1)2Aau0(AA0)2+RaAb,u2(AA2)2Ab

      (2)

      where the parameters a and b=(A0a)(A0A1)A2A0+A0 are the smooth connection points, A1 and A2 are the positions of the light and heavy fragment peaks of the pre-neutron fragment mass yields, respectively, and A0 denotes the corresponding position at the symmetric fission point. The potential parameters u0, u1, and u2 are expressed as

      u0=R(A0a)(A0A1),u1=R(A0A1)(aA1),u2=R(A2A0)(A2b).

      (3)

      The pre-neutron fragment mass yield of binary fission should be normalized to 200%. Therefore, the normalization constant C can be analytically expressed as

      C=200%0exp[U(A)]dA=200%I0+I1+I2,

      (4)

      with

      I0=πeR2u0{erfi[(A0a)u0]+erfi[(bA0)u0]},I1=π2u1{erf[(aA1)u1]+erf(A1u1)},I2=π2u2{1+erf[(A2b)u2]}.

      (5)

      With increasing incident neutron energy, the valley values of the pre-neutron fragment mass yields obviously increase, in contrast with the peak values. That is, the peak-to-valley ratios decrease with increasing incident neutron energy. The parameter R in Eq. (2) can easily be obtained as

      R=lnY(A1)Y(A0),

      (6)

      where Y(A1) and Y(A0) represent the peak and valley values of the pre-neutron fragment mass yields, respectively. Thus, all of the model parameters are functions of a,A0,A1,A2, and R. In particular, the parameter a can be obtained by solving Eq. (1) and only defining A=A2. Therefore, there are no artificially adjustable parameters.

      With the increase in incident neutron energy, the excitation energy of the compound nucleus increases, and a few neutrons (known as pre-neutrons, that is, evaporation neutrons) are evaporated before scission. The number of pre-neutrons can be obtained from the corresponding multi-chance fission cross sections. The relationship between the positions of the peak and valley and the mass number AFN of the fissile nucleus is expressed as

      {A1=AFNA2,A0=(A1+A2)/2,A2=Ags2˜n(En),

      (7)

      where Ags2 140 denotes the peak position of the heavy fragment mass yields for neutron-induced actinide fission at low incident energies [15], ˜n(En) is the number of pre-neutrons corresponding to the different multi-chance fission channels, and En denotes the incident energy of the neutron [15]. When the number of pre-neutrons is zero (that is, ˜n(En)=0), the mass number of the compound nucleus ACN is that of the fissile nucleus, that is, AFN=ACN.

    III.   RESULTS AND ANALYSIS

      A.   Pre-neutron fragment mass yields for the 235U(n, f) reaction

    • With increasing incident neutron energy, the excitation energy of the compound nucleus increases; hence, a few pre-neutrons evaporate before scission. The number of pre-neutrons can be obtained from the corresponding multi-chance fission cross sections. Therefore, the fission cross sections of the 235U(n, f) reaction are investigated. Figure 1 shows the experimental and evaluated fission cross sections of the 235U(n, f) reaction. The experimental data are obtained from Refs. [1719]. The evaluated values of ENDF/B-VIII, indicated by the pink dotted line, are recommended as standard cross sections for this reaction [20]. There is an obvious step phenomenon, and every step corresponds to the number of pre-neutrons. From Fig. 1, it can be established that the number ˜n(En) of pre-neutrons can be roughly expressed as follows:

      Figure 1.  (color online) Fission cross section of the 235U(n, f) reaction for incident neutron energies up to 20 MeV. The experimental data are obtained from Refs. [17] (circles), [18] (triangles), and [19] (squares). The pink dotted line represents the evaluated values of ENDF/B-VIII [20]. The perpendicular dashed lines label the incident energy regions corresponding to different multi-chance fission channels, such as (n, f), (n, nf), and (n, 2nf).

      ˜n(En)={0,En<5.5 MeV,1,5.5En<11.5 MeV,2,11.5En<20 MeV.

      (8)

      The experimental values of Y(A0) and Y(A1) [2124], which denote the valleys and peaks of the pre-neutron fragment mass yields for the 235U(n, f) reaction, are fitted to the functions of the incident neutron energies,

      Y(A0)=α0+β0Exp(En/λ0),Y(A1)=α1+β1Exp(En/λ1),

      (9)

      where the parameters α0, β0, λ0, α1, β1, and λ1 are finally determined by the experimental data of Y(A1) and Y(A0) and listed in Table 1. The exponential energy dependences of the peaks Y(A1) and valleys Y(A0) of the pre-neutron fragment mass yields for the 235U(n, f) reaction are shown in Fig. 2 as red solid lines. The experimental data are obtained from Refs. [21] (rhombus), [22] (squares), [23] (circles), and [24] (triangles). Furthermore, Eq. (9) approximately equals the linear results in the incident energy region 0.5 <En< 5.5 MeV, as shown in Ref. [15].

      αiβiλi
      235U(n, f)i=0−0.2150.218−10.754
      i=12.0524.62043.161
      239Pu(n, f)i=20.0390.0240
      i=36.060−0.1350

      Table 1.  Parameters of the exponential dependences in Eq. (9) and Eq. (11) for this study.

      Figure 2.  (color online) Peak Y(A1) and valley Y(A0) of the pre-neutron fragment mass yields for the 235U(n, f) reaction as a function of incident neutron energy. The experimental data are obtained from Refs. [21] (rhombus), [22] (squares), [23] (circles), and [24] (triangles). The red solid lines represent the results of this study.

      Using Eqs. (1)–(9) and Table 1, the calculated pre-neutron fragment mass yields for the 235U(n, f) reaction are shown by the red solid lines in Fig. 3 and Fig. 4 at different incident energies. We can see that the theoretical results of this study agree well with the experimental data, which are obtained from Refs. [2126]. In particular, Fig. 4 shows the results at the thermal neutron energy, which have not been previously published. The results of the GEF model [14] are also simultaneously compared in Fig. 3 and Fig. 4 as blue dash lines. The root mean square deviations at different incident energies of this study and the GEF model are listed in Table 2. As an example, the potential parameters at different incident energies (En= 0.0253 eV, 6 MeV, and 14 MeV) are listed in Table 3. This indicates that the EFP model can reasonably describe existing data on pre-neutron fragment mass yields for the 235U(n, f) reaction at incident energies from thermal up to 20 MeV.

      Figure 3.  (color online) Pre-neutron fragment mass yields for the 235U(n, f) reaction at incident energies of 5.5 MeV (a), 6 MeV (b), 7.04 MeV (c), and 15.5 MeV (d). The scattered symbols denote the experimental data, which are taken from Refs. [21] (rhombus), [22] (squares), [23] (circles), and [24] (triangles). The blue dashed and red solid curves denote the calculated results of the GEF model[14] and this study, respectively.

      Figure 4.  (color online) Same as Fig. 3 but at thermal neutron energy. The scattered symbols denote the experimental data, which are taken from Refs. [24] (triangles), [25] (circles), and [26] (squares).

      En/MeV235U(n, f)239Pu(n, f)Model
      2.53 ×1080.5420.554GEF
      0.5360.376This study
      0.720.585GEF
      0.299This study
      1.720.503GEF
      0.324This study
      2.720.457GEF
      0.310This study
      4.480.540GEF
      0.267This study
      5.050.546GEF
      0.355This study
      60.397GEF
      0.239This study
      7.040.499GEF
      0.330This study
      15.50.338GEF
      0.179This study

      Table 2.  Root mean square deviation at different incident energies compared with existing experimental data and the theoretical results.

      ReactionEn/MeVabI0I1I2R
      235U(n, f)2.53×108117.99118.010.00014.20714.2077.529
      6111.64124.360.53116.18116.1813.596
      14109.71126.292.55217.24617.2462.216
      239Pu(n, f)2.53×108119.99120.010.00015.78915.7895.031
      6119.99120.010.00018.37118.3713.325
      14119.99120.010.00221.53121.5312.373

      Table 3.  Potential parameters adopted in this study.

    • B.   Pre-neutron fragment mass yields for the 239Pu(n, f) reaction

    • Besides the 235U(n, f) reaction, the EFP model was also applied to the 239Pu(n, f) reaction. First, the fission cross sections of the 239Pu(n, f) reaction are investigated, as shown in Fig. 5. The experimental data are obtained from Refs. [2729], and the pink dotted line represents the evaluation results of ENDF/B-VIII [20], which are recommended as standard cross sections. There is also an obvious step phenomenon, and every step corresponds to the number of pre-neutrons. From Fig. 5, it can be established that the number ˜n(En) of pre-neutrons can be roughly expressed as follows:

      Figure 5.  (color online) Same as Fig. 1 but for the 239Pu(n, f) reaction. The experimental data are obtained from Refs. [27] (circles), [28] (squares), and [29] (triangles).

      ˜n(En)={0,En<5.6 MeV,1,5.6En<12 MeV,2,12En<20 MeV.

      (10)

      Based on monoenergetic experimental data from Ref. [30], the values of the valleys Y(A0) and peaks Y(A1) of the pre-neutron fragment mass yields for the 239Pu(n, f) reaction are fitted as linear functions of the incident neutron energies [15],

      Y(A0)=α2+β2En,Y(A1)=α3+β3En,

      (11)

      where the parameters α2, β2, α3, and β3 are finally determined by the experimental data of Y(A1) and Y(A0) and listed in Table 1. Because of scarce measured data, the linear energy dependences of the peaks Y(A1) and valleys Y(A0) of the pre-neutron fragment mass yields for the 239Pu(n, f) reaction are shown in Fig. 6 as red solid lines. The experimental data are obtained from Refs. [30]. Eq. (11) and its parameters can reproduce the results of Ref. [15]. It is expected that there are some measurements at higher incident energies to verify these linear dependences in the future.

      Figure 6.  (color online) Same as Fig. 2 but for the 239Pu (n, f) reaction. The experimental data are obtained from Ref. [30] (squares).

      According to Eqs. (1)–(7), (10), (11), and Table 1, the calculated pre-neutron fragment mass yields for the 239Pu(n, f) reaction can reproduce our previous results [15] at different incident energies (En= 0.72, 1.72, 2.72, and 4.48 MeV). Furthermore, Fig. 7 shows the calculated results at the thermal neutron energy, which have not previously been published. We can see that the theoretical results of this study agree well with the experimental data, which are obtained from Refs. [30] (triangles), [31] (circles), and [32] (squares). Simultaneously, the results of the GEF model [14] are also compared in Fig. 7 as blue dash lines. The root mean square deviations at different incident energies in this study and the GEF model are listed in Table 2. As an example, the potential parameters at different incident energies (En= 0.0253 eV, 6 MeV, and 14 MeV) are listed in Table 3. This indicates that the EFP model can also reasonably describe the pre-neutron fragment mass yields for the 239Pu(n, f) reaction at incident energies from thermal up to 20 MeV, as well as those of the 235U(n, f) reaction. Figure 8 also gives the predicted results of both this study and the GEF model at an incident energy of 14 MeV.

      Figure 7.  (color online) Pre-neutron fragment mass yields at thermal neutron energies for the 239Pu(n, f) reaction. The scattered symbols denote the experimental data, which are taken from Refs. [30] (triangles), [31] (circles), and [32] (squares). The blue dashed and red solid curves denote the calculated results of the GEF model and this study, respectively.

    IV.   CONCLUSION
    • In this paper, the pre-neutron fragment mass yields for the 235U(n, f) and 239Pu(n, f) reactions at incident energies from thermal up to 20 MeV are calculated on the basis of the EFP model, the potential parameters of which are uniquely obtained from the existing experimental values of the peaks and valleys. The energy dependences of the peaks and valleys are reasonably described by explicit exponential and linear functions for the 235U(n, f) and 239Pu(n, f) reactions, respectively. The energy dependences of the number of evaporation neutrons are also derived from multi-chance fission cross sections, which play a crucial role in the reasonable description of pre-neutron fragment mass yields. The calculated results of this study can not only reproduce our previous results but also agree well with existing data at incident energies from thermal up to 20 MeV. This extended energy region is exactly what nuclear data libraries need.

      Furthermore, the calculated results of the GEF model are also compared. The root mean square deviations of this study and the GEF model at different incident energies are listed in Table 2. As shown, the EFP model in this study can feasibly describe pre-neutron fragment yields quantitatively. The results show that the feasibility of this EFP model is verified in this extended energy region. Therefore, the EFP model is expected to be a tool used to set up the nuclear database files of pre-neutron fragment mass yields for the 235U(n, f) and 239Pu(n, f) reactions at incident energies from thermal up to 20 MeV.

    ACKNOWLEDGEMENTS
    • We thank Profs. Ning Wang, Li Ou, Min Liu, Dr. Yuting Rong, and anonymous referees for valuable suggestions.

Reference (32)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return