_{1}

^{*}

In this work, new cross section formulae for (n, p) and (n, d) reactions of He-3 have been investigated as a function of incident neutron energy. The new Empirical formulae were produced by using the least squares method to the experimental cross sections data, which were taken from Experimental Nuclear Reaction Data EXFOR Database Version of 2021. Several functions were examined to choose the best one that fits the data. Statistical hypothesis testing was used to insure how well the suggested equations fit the set of data. Two statistical indicators were used for each case for goodness-of-fit. Very high compatibility was found between the empirical values and the experimental data for both reactions.

The nucleus helium-3 (He-3) which proposed for the first time by the Mark Oliphant in 1934 [

Cross section data may be needed during calculations of nuclear systems, medical applications, nuclear industrial studies, or other nuclear fields. Some nuclear database libraries such as the Evaluated Nuclear Data File (ENDF), Japanese Evaluated Nuclear Data Library (JENDL), TALYS-based Evaluated Nuclear Data Library (TENDL), and Joint Evaluated Fission and Fusion File (JEFF) are available for (n, d) reactions but theses libraries do not cover all the required ranges of cross sections [

“Curve Fitting” term indicates to the mathematical method that one can find the formula which can be represented as: y = f ( x ) that can be the best fit with the corresponding y-values and the set of observed x-values [

For the mathematical method, the creation of the best-fit formulas for the observed set of points was to reduce the sum of squares of the residuals of the points from the curve. This method is called the least square method [

By assuming that the observed data are in the form of Cartesian pairs ( x i , y i ) , 1 ≤ i ≤ N , and the suggested empirical function Y ( X i ) is assumed to be a straight line with the equation:

Y ( X i ) = a ∗ X i + b (1)

The i^{th} error Er(i) of the estimation data is then equaling to:

E r ( i ) = Y ( X i ) − a − b X i (2)

Then, total squared error “E” can be defined as:

E = ∑ i = 1 N ( y i − a − b X i ) 2 (3)

where N is the number of the observed pairs of data and E’s subscript refers to the solution’s order.

The lower the value of E, the closer the proposed solution is to the observed data. For calculation of a andb constants, E must be very small. By partial derivation of Equation (3) with respect to a and b two times and making them equal to zero, the constants a and b can be deduced.

∂ ∂ a E = ∂ ∂ b E = 0 (4)

∂ ∂ a ( E ) = ∂ ∂ a ∑ i = 1 N ( y i − a − b X i ) 2 = 0 (5)

∂ 2 ∂ b ( E ) = − 2 ∗ ∑ i = 1 N y i − a − b X i = 0 (6)

So

∑ i = 1 N y i = a ∗ N + b ∗ ∑ i = 1 N E r ( i ) (7)

∑ i = 1 N y i ∗ E r ( i ) = a ∗ ∑ i = 1 N E r ( i ) + b ∗ ∑ i = 1 N E r ( i ) 2 (8)

Parameters a and b are found from the solution of Equations (7) and (8). For general case, least-squares method is applied to reach the best solutions by finding the best possible values for a group of parameters.

Parametric fitting includes discovering parameters for at least one show that you fit to information. The information is thought to be measurable in nature and is separated into two segments: deterministic and arbitrary. The deterministic part is given by a parametric model and the arbitrary segment is frequently depicted as blunder related with the data. We can say that the observed data is equal to the sum of both parametric model and the error.

A goodness-of-fit test or indicator is performed to statistically examine whether the estimated formula gives a good description of the observed data. Nash-Sutcliffe efficiency indicator NSE and root mean square error indicator RMSE are two of the most popular indicators for goodness of fit [

If “N” is the number of observed points, then RMSE is defined as:

RMSE = 1 N ∑ 1 N ( y i − Y i ) 2 (9)

Its value ranges from zero to infinity, where zero value indicates the best fit.

The second goodness-of-fit indicator is the dimensionless indicator NSE and can be calculated as [

NSE = 1 − [ ∑ 1 N ( y i − Y i ) 2 / ∑ 1 N ( y i − Y mean ) 2 ] (10)

where y mean represented the mean of the estimated or fitted values of the cross section.

NSE ranges from minus infinity to 1, where unity NSE indicates the best fit.

During calculation of the lifetime of neutron [

The suggested empirical formulae for (n, p) and (n, d) reactions of He-3 as a function of only incident neutron energy depend on fitting of the observed data by using parametric fitting method. In this method, many different formulae were tried until finding the best parameters of the suggested equation that give the best fit of the observed data. Parametric fitting helped us to find parameters for many pre-suggested formulae and only those parameters that gave physical significance were approved. There are no theories to dictate the suggested formulae. Power functions, exponential functions, polynomial function sgaussian functions, and Fourier functions, with different degree, were used to reach the best fit.

Observed cross section data for ^{3}He (n, p) ^{3}H and ^{3}He (n, d) ^{2}H described above from the experimental nuclear reaction data EXFOR database of version of 2021 [

1) He-3(n, d) H-2

σ ( E ) = a ∗ e b ∗ E + c ∗ E 2.1 + d ∗ E 0.5 + e (11)

wherea = −786, b = −0.6132, c = −0.1607, d = 36.7, and e = −20.13.

2) He-3(n, p) H-3

σ ( E ) = a E + b ∗ E + c (12)

where a = 0.8458, b = 1.189, c = -1.174, and E ≤ 1 MeV.

σ ( E ) = a 1 e ( E − b 1 c 1 ) 2 + a 2 e ( E − b 2 c 2 ) 2 + a 3 e ( E − b 3 c 3 ) 2 + a 4 e ( E − b 4 c 4 ) 2 + a 5 e ( E − b 5 c 5 ) 2 + a 6 e ( E − b 6 c 6 ) 2 (13)

where a_{1} = 0.1342, b_{1} = 1.611, c_{1} = 0.679, a_{2} = 0.05731, b_{2} = 2.468, c_{2} = 0.7579, a_{3} = 0.6842, b_{3} = 1.899, c_{3} = 3.229, a_{4} = -0.462, b_{4} = 2.276, c_{4} = 3.161, a_{5} = 0.08333, b_{5} = 13.5, c_{5} = 18.76, a_{6} = 1.947,b_{6} = -18.07, c_{6} = 16.14, and 1 ≤ E ≤ 20 MeV.

Equation (13) gives the new empirical cross section (mb) formula as a function of energy (MeV) of the incident neutrons for (n, d), while both Equations (14) and (15) give the new empirical cross section formulae for (n, p) according to two regions of incident neutrons energy.

Also, the discrepancies between the observed data and the calculated were included in

Δ ( i ) = σ observed ( E i ) − σ calculated ( E i ) σ observed ( E i ) (14)

By applying Equation (9) and Equation (10) to the observed and computed data, the RMSE and NSE values for the three equations are calculated and given in

The values of RMSE and NSE of the curve fitting of the cross-sections of the two interactions showed a very high accuracy between the experimental values and the experimental values, given the values of the most important statistical

Nuclear Reaction | Goodness of fit indicators | |
---|---|---|

RMSE Equation (9) | NSE Equation (10) | |

He-3 (n, d) H-2 | 1.16 | 0.9981 |

He-3 (n, p) H-3 (E ≤ 1 MeV) | 0.4329 | 1 |

He-3 (n, p) H-3 (1 ≤ E ≤ 20 MeV) | 0.001036 | 1 |

indicator NSE, which are very close to the integer one in Equation (11) and equal to the integer one in the two Equations (12) and (13). Also, the values of the other indicator RMSE are ideal for Equations (12) and (13) and are considered acceptable for Equation (11).

The first step for evaluating the goodness of the new empirical formulae is the visual examination of the fitted curve displayed in couples with the observed data. And it is cleared from both

The authors declare no conflicts of interest regarding the publication of this paper.

Alramady, A.M. (2021) Empirical Neutron Cross Section Formulae for (n, p) and (n, d) Reactions of He-3 Target. Journal of Applied Mathematics and Physics, 9, 1434-1441. https://doi.org/10.4236/jamp.2021.97096