Issue 
EPJ Nuclear Sci. Technol.
Volume 2, 2016



Article Number  32  
Number of page(s)  16  
DOI  https://doi.org/10.1051/epjn/2016025  
Published online  23 August 2016 
https://doi.org/10.1051/epjn/2016025
Regular Article
Measurements of the effective cumulative fission yields of ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd for ^{235}U in the PHENIX fast reactor
^{1}
CEA, DEN Cadarache, 13108 Saint Paul Les Durance, France
^{2}
Centre d’Etudes Nucléaires Bordeaux Gradignan, CNRS/IN2P3, Univ. Bordeaux 1, 33175 Gradignan, France
^{3}
National Nuclear Laboratory, B170 Sellafield Works, Seascale, Cumbria CA20 1PG, UK
^{⁎} email: gilles.noguere@cea.fr
Received:
8
March
2016
Received in final form:
13
June
2016
Accepted:
20
June
2016
Published online: 23 August 2016
The effective Neodymium cumulative fission yields for ^{235}U have been measured in the fast reactor PHENIX relatively to the ^{235}U fission crosssection. The data were derived from isotoperatio measurements obtained in the frame of the PROFIL1, PROFIL2A and PROFIL2B programs. The interpretations of the experimental programs were performed with the ERANOS code in association with the Joint Evaluated Fission and Fusion library JEFF3.1.1. Final results for ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd were 5.61%, 3.70%, 2.83%, 1.64% and 0.66%, respectively. The relative uncertainties attached to each of the cumulative fission yields lie between 2.1% and 2.4%. The main source of uncertainty is due to the fluence scaling procedure (<2%). The uncertainties on the Neodymium capture crosssections provide a contribution lower than 1%. The energy dependence of the fission yields was studied with the GEF code from the thermal energy to 20 MeV. Neutron spectrum average corrections, deduced from GEF calculations, were applied to our effective fission yields with the aim of estimating fission yields at 400 keV and 500 keV, as given in the International Evaluated Nuclear Data Files (JEFF, ENDF/B and JENDL). The neutron spectrum average correction calculated for the PROFIL results remains lower than 1.5%.
© E. Privas et al., published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The present work reports experimental cumulative fission yields Y̅_{c}(^{A}Nd) that represent the total number of atoms of ^{A}Nd produced over all time after one fission of ^{235}U and averaged over the fastneutron spectrum of the PHENIX reactor. Such experimental results are needed for the computer simulation of reactors, fuel cycles and waste management [1,2]. Values and relative uncertainties from current evaluations of interest for this work are listed in Table 1. A good knowledge of these quantities enables estimation of the spent fuel inventories and resultant quantities such as decay heat, delayed neutron and gamma ray emissions. They are also important in understanding fuel performance, such as reactivity loss during reactor operation. Apart from burnup and safety calculation, the isotope inventories can be used to characterize a spent fuel. The percentage burnup of actinide in the fuel can be determined by analyzing the fission product inventories. One of the fission yields routinely used for reactor applications is ^{148}Nd. The burnup of fuel can be estimated using the stable nuclide ^{148}Nd, which is determined by destructive assay. The precursors of ^{148}Nd are all short lived and have small crosssections and thus the cumulative yield can be used to determine the number of fissions that has taken place in the fuel.
Previous interpretations of the PROFIL1, PROFIL2A and PROFIL2B integral experiments carried out in the PHENIX reactor (CEA Marcoule, France) with the ERANOS code [6,7] in association with the evaluated nuclear data library JEFF3.0 and JEFF3.1 showed some significant discrepancies for the cumulative fission yields of several Neodymium isotopes for the fast fission of ^{235}U [8,9]. The average calculatedtoexperimental ratios C/E associated to the Neodymium produced by fission in several samples are gathered in Figure 1. Conspicuously, there is a large drift with respect to the atomic mass number A from A = 143 to A = 150, far exceeding the claimed experimental uncertainties on fission yields, which lie in the range between 1% and 2.4%. Such a drift is not observed for the ^{239}Pu and ^{241}Pu samples.
A reanalysis of the PROFIL experiments was performed with the ERANOS code by using the nuclear data of the JEFF3.1.1 library [5,10]. A similar drift of the average C/E values related to the Nd isotopes measured in the ^{235}U samples was observed. This result confirms that the fast cumulative fission yields of the Nd isotopes recommended in the JEFF3.1.1 library need to be revised with improved experimental values. The present work details how experimental fission yields with realistic uncertainties can be extracted from the integral trends provided by the PROFIL experiments. The originality of our data reduction and uncertainty analysis relies on the use of the Generalized Perturbation Theory [11], implemented in the ERANOS code, and of sensitivity coefficients calculated by direct perturbations of the nuclear data and experimental correction parameters.
The effective cumulative fission yields, given for the PHENIX neutron spectrum, were derived from the C/E values with an uncertainty lying between 2.1% and 2.4%. Special care was taken in evaluating the uncertainties related to the use of a new ERANOS calculation scheme in association with an improved fluence scaling procedure, needed to normalize the experimental results to the ^{235}U fission crosssection. Biases caused by the Neodymium neutron capture reactions during irradiation were studied and quantified by direct perturbations of the capture crosssections of the ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd isotopes. Additional calculations performed with the GEF code [12,13] allowed accounting for the systematic behavior of the fission yields with the incident neutron energy. Final results are compared with experimental data reported in the literature and with evaluated data given at 400 keV or 500 keV in the JENDL, ENDF/B and JEFF nuclear data libraries.
Fig. 1 Average C/E ratios for the prediction of the Neodymium buildup in the PROFIL1 and PROFIL2 programs for ^{235}U, ^{239}Pu and ^{241}Pu obtained with the ERANOS code and the JEFF3.1 library [9]. 
2 Definition of the effective cumulative fission yields
A fission product is defined symbolically by the notation (A, Z, I) where A and Z are respectively the mass number and the atomic number, and I indicates the isomeric state. The ground state is denoted by I = 0 and I = 1, 2, … represents the 1st, 2nd, … isomeric states. If a fission product has no isomers, or if one is referring to the sum of yields of all its isomers, the notation (A, Z) is used. Using this terminology, we can distinguish the independent and cumulative fission yields. For a given incident neutron energy E, the independent fission yield Y_{post}(A, Z, I, E) is the number of atoms of a specific nuclide produced directly per 100 fission reactions, after prompt neutrons evaporation but before any radioactive decay. To take into account radioactive decay, cumulative fission yield Y_{c}(A, Z, I, E) is used. It represents the number of atoms of a specific nuclide produced directly and via decay of precursors per 100 fission reactions. The cumulative fission yields can be defined by the expression [14]:(1)where N is the dimension of the whole fission fragment inventory, i indicates a generic triplet (A_{i}, Z_{i}, I_{i}) and b_{ji}(E) is the branching ratio, which gives the probability that an isomer (A_{j}, Z_{j}, I_{j}) decays into (A_{i}, Z_{i}, I_{i}).
In the present work, we report effective cumulative fission yields Y̅_{c} measured in ^{235}U samples, which have been irradiated in the fast reactor PHENIX. They can be obtained from:(2)in which σ_{f}(E) is the ^{235}U fission crosssection and ϕ(E) stands for the fastneutron spectrum representative of the PROFIL experiment. The maximum energy E_{max} is set to 20 MeV.
As stable Neodymium isotopes (Z = 60) in the ground state deformation (I = 0) are considered in this work, the notation Y_{c}(A, Z, I, E) and Y̅_{c}(A,Z,I) will be replaced throughout the document by the explicit notation Y_{c}(^{A}Nd,E) and Y̅_{c}(^{A}Nd).
3 Description of the PROFIL1 and PROFIL2 programs
The PROFIL programs were carried out in the fast reactor PHENIX (CEA Marcoule, France) from 1974 to 1980. They are mainly used to provide nuclear data feedbacks on neutron crosssections ((n,γ), (n,2n), (n,f)) of several fission products and actinides. A detailed study of the integral results related to the cumulated fission yields of ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd for ^{235}U were never reported. The principle of the experimental setup is briefly described below. More information can be found in previous papers [8,9].
3.1 Principle of the experiments
The PROFIL experiments were designed to study the fast energy range of the neutron crosssections, called “continuum” in Figure 2. In the same figure, the ^{235}U capture and fission crosssections are compared to the fastneutron spectrum calculated by using the ECCO lattice code of the ERANOS system [6,7] with 1968 energy groups and a critical buckling. Fission yields for the main actinides such as ^{235}U, ^{238}U, ^{238}Pu, ^{239}Pu, ^{240}Pu, ^{241}Pu, ^{242}Pu and ^{241}Am are also available. Only results obtained for the fission of ^{235}U are presented in this paper.
The principle of the experiments consists of irradiating samples containing a small amount of pure isotope in a fastneutron flux. The two experiments were carried out four decades ago during the first cycles of the 250 MWe sodiumcooled fast reactor PHENIX (CEA Marcoule, France). A simplified radial view of the core is presented in Figure 3. The first phase (PROFIL1) was composed of a single experimental pin loaded in a standard PHENIX assembly placed at the center of the core. The pin contained 46 separate samples. The second phase (PROFIL2A and PROFIL2B) was a more ambitious program with two experimental pins, labelled A and B, placed in the inner core of PHENIX. Each pin contained 42 separate samples. Figure 3 also shows the pin location in their respective fuel subassemblies. A total of 130 samples were irradiated during four cycles for PROFIL1 and PROFIL2 (>10 months).
Samples loaded in the PROFIL pins were cylindrical containers with an outer diameter of 5.5 mm and a height of 8 mm. An example of a stainless steel double container is shown in Figure 3. For ^{235}U, a total of 6 and 7 regularly spaced samples were loaded in the PROFIL1 and in each PROFIL2 pins, respectively. As shown in the axial layout of the PROFIL1 pin (Fig. 4), the location of the ^{235}U samples was chosen to provide accurate information for fluence scaling issue. We were able to use the experimental results of the 6 samples loaded in the PROFIL1 pin, while only 3 samples (of the 7 samples) were analyzed in each PROFIL2 pins. In total, up to 12 integral trends were collected for determining the cumulative fission yield of ^{148}Nd and 9 integral trends for the other Neodymium isotopes.
Fig. 2 ^{235}U capture and fission crosssections available in the JEFF3.1.1 library compared to a neutron spectrum, in arbitrary unit, representative of the PROFIL experiment (red histogram). 
Fig. 3 Radial core geometry of the fast reactor PHENIX, location of the subassemblies containing the PROFIL pins, axial layout of the PROFIL1 pin and example of a stainless steel double container designed for the PROFIL irradiation. 
Fig. 4 Axial layout of the PROFIL1 pin with the position of the ^{235}U samples used in this work. The sample number 19 is located at z = 0 mm, which corresponds to the core midplane of the PHENIX reactor. 
3.2 Isotoperatio measurements
The experimental results used for the interpretation of the PROFIL experiments are the isotopic composition of the samples before and after the irradiation period. As indicated in Section 3.1, the time interval δt is of the order of 10 months. The variation, or equivalently the buildup, of the number of atoms provides integral information on the capture, fission and (n,2n) crosssections.
We focus the present work on the isotopic composition of the ^{235}U samples before and after the irradiation time δt. At the beginning of the experiment, the ^{235}U samples are uranium oxide samples composed of ^{235}U (≃90%), ^{234}U (≃1%) and ^{236}U (>0.5%). After the irradiation period, we observe a decrease of ≃20% of the ^{235}U content.
The determination of the sample composition was achieved by Inductively Coupled Plasma Mass Spectrometry (ICPMS) [15]. For each sample k, mass spectrometry measurements supplied high precision isotopic ratios whose relative uncertainties range from 0.1% to 0.2%. The variations between t_{0} and t_{0} + δt define the following experimental quantities:(3)
Experimentally, the isotopic variations , , and have to be combined with the ratio in order to deduce the isotopic variations of interest , , and . The corresponding values calculated with the ERANOS code are denoted by . Throughout the document, the C/E ratios have the generic form:(4)
In references [8,9], the quantities are assumed to be normally distributed and independent. As a consequence, the weighted average over the samples k of the calculatedtoexperimental ratios, formally defined as:(5)are characterized by unrealistic small relative uncertainties. In practice, correlations between the uncertainties of the different ICPMS results were not provided by the experimentalists. Mathematical solutions to calculate appropriate uncertainties were discussed in reference [16] in which the experimental uncertainty on each isotopic ratio is treated as a “statistical” contribution while the “systematic” part is driven by two fluence scaling parameters.
3.3 Fluence scaling issue
The interpretation of the PROFIL experiments is performed with the deterministic code ERANOS [6,7]. The calculation scheme relies on several hypotheses that are used to obtain a mean neutron flux over the irradiation period δt. The dispersion of the C/E ratios, observed for the ^{235}U samples, suggests that such a raw treatment overestimates the magnitude of the global fluence. The origins of the observed dispersions could be due to the approximations used in the calculation scheme but also to the influence of the surrounding materials, such as the steel containers, that could locally change the neutron spectrum. The strategy for solving the inconsistency between the calculations and the experimental results consists in normalizing the results utilizing the ^{235}U fission crosssection. For that purpose, two free parameters were introduced in the ERANOS calculations. The first parameter was introduced to shift the axial shape of the flux and the second parameter was introduced to normalize the global fluence level. Optimal parameter values are determined by using as reference the following fission indicator:(6)so that the related calculatedtoexperimental ratios satisfy the following constraint:(7)
For example, in the case of the PROFIL2A configuration, the sensitivity of this experimental ratio to the ^{235}U fission crosssection is close to unity (1.07), whereas the sensitivity is relatively small to the ^{238}U capture crosssection (0.14) and negligible to the ^{235,236}U capture crosssections (−0.04 and 0.01, respectively). The ^{235}U fission crosssection is one of the major reactions investigated within the “neutron crosssection standards” group of the International Atomic Energy Agency (IAEA) [17]. This crosssection is considered as a standard at thermal energy (25.3 meV) and from 0.15 MeV to 200 MeV with relative uncertainties lower than 1%. Below 0.15 MeV, crosssection values reported by IAEA are given as a recommendation only. As calculated in reference [16], the accuracy of the fluence scaling, due to the precisions of the ^{235}U fission crosssection, is close to 1.5% in average.
4 Integral results on Neodymium produced by fission
Integral results obtained from the interpretation of the PROFIL experiments with the ERANOS code are expressed in terms of calculatedtoexperimental ratios (Eq. (4)), normalized to the ^{235}U fission crosssection by using the fluence scaling operations shortly explained in Section 3.3. The main sources of uncertainties discussed in the following sections are the uncertainties related to (i) the calculation scheme, (ii) the calculation of the mean value (Eq. (5)), (iii) the calculation of the axial shape of the neutron flux, (iv) the fluence scaling and (v) the Nd capture crosssections.
4.1 Calculatedtoexperimental ratios obtained with JEFF3.1.1
The C/E values obtained with the ERANOS code by using the JEFF3.1 and JEFF3.1.1 libraries are listed in Table 2 for each ^{235}U sample. The two JEFF libraries share the same Nd cumulative fission yields [5]. The main difference between the two sets of results arises from the fluence scaling. In the previous calculations, the fluence scaling parameters (Sect. 3.3) were finetuned by trial and error. In the present work, their values were automatically adjusted so that the theoretical ratios of equation (6) agree with the experimental values within the limit of the uncertainties [16]. No significant differences are observed between the JEFF3.1 and JEFF3.1.1 results. The agreement lies between 0.1% and 0.5% for the PROFIL and PROFIL2 results, respectively. A contribution of 0.5% will be added in the final uncertainty in order to account for the possible biases due to the calculation scheme.
The average ratios ⟨C/E⟩ reported in Table 2 are displayed in Figure 5 as a function of the mass of the Neodymium isotopes. These results confirm the trends reported in the previous works (Fig. 1). A satisfactory agreement is observed between the two PROFIL2 pins (A and B). The discrepancies range from 0.5% (for the ^{145}Nd and ^{146}Nd buildup) to 1.3% (for the ^{148}Nd buildup). The main interesting result shown in Figure 5 is the sizeable difference between the PROFIL1 and PROFIL2 results in the case of the ^{150}Nd buildup. The origin of such a difference is difficult to explain with the available experimental information.
C/E ratios for the prediction of the Neodymium buildup, after fluence scaling. Results obtained in the present work with the JEFF3.1.1 library are compared with those obtained with the JEFF3.1 library [18]. The 1st and 2nd weighted mean values for PROFIL were calculated with and without the C/E results obtained for the ^{235}U sample number 19 located in the core midplane (z = 0 mm). The quoted uncertainties are the experimental uncertainties coming from the chemical analysis of the ^{235}U samples.
Fig. 5 Average ⟨C/E⟩ ratios for the prediction of the Neodymium buildup in the PROFIL1, PROFIL2A and PROFIL2B programs for ^{235}U obtained with the ERANOS code and the JEFF3.1.1 library. Individual results are listed in Table 2. 
4.2 Uncertainties on the average values
A closer inspection of the integral results has shown that the average values in the PROFIL1 experiment are slightly dependent on the C/E values reported for the ^{235}U sample number 19. This sample is located in the PROFIL1 pin at z = 0 mm, which corresponds to the core midplane of the PHENIX reactor (Fig. 4). As a complementary information, Table 3 reports the C/E values obtained for the ^{144}Nd + ^{144}Ce buildup. Results found for the ^{235}U sample number 19 are always significantly greater than the other ^{235}U samples. A similar trend with a lower amplitude can be observed in the plots of Figure 6. They display the C/E ratios related to the prediction of the ^{143}Nd, ^{144}Nd + ^{144}Ce, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd buildup in the ^{235}U samples of the PROFIL1 experiment as a function of the sample position in the pin. In each case, the highest value is reached at z = 0 mm, which makes questionable the accuracy attached to the sample number 19. The average ratios ⟨C/E⟩ calculated with and without the C/E results obtained for the ^{235}U sample number 19 are reported in Table 2. The differences between the two sets of average values range from 0.2% to 0.3%. In the following, the 2nd set of mean values reported in Table 2 (ignoring the sample 19 results) will be considered for the determination of the cumulative fission yields. An additional uncertainty of 0.3%, associated to this choice, will be taken into account.
Fig. 6 Comparison of the individual C/E results obtained with JEFF3.1.1 (open circle) and JEFF3.1 (black circle) for the prediction of the ^{143}Nd, ^{144}Nd + ^{144}Ce, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd buildup in the ^{235}U samples of the PROFIL1 experiment as a function of the sample position. The abscise z = 0 mm corresponds to the core midplane of the PHENIX reactor. The solid line is a parabolic curve fitted to the data acting as a eye guide. 
4.3 Uncertainties due to the axial shape of the neutron flux
A suspicious experimental trend can also be observed in Figure 6. In each plot, there is an apparent trend, fitted to a parabolic curve as an eye guide, with axial position. Such a behavior as a function of the sample position can be seen as a systematic trend between the middle and the extremity of the PROFIL pin, which is usually associated with the capability of our calculation scheme to correctly reproduce the axial shape of the neutron flux and the edge effects far away from the core midplane. This could be an explanation why the C/E results for the sample number 19 (located in the core midplane of the PHENIX reactor) are systematically higher (see Sect. 4.2). Full 4D MonteCarlo simulations for correcting this trend were not investigated in the present work. Instead, we introduce an additional source of uncertainty of 1% related to the axial shape of the neutron flux.
4.4 Uncertainties due to the fluence scaling
The recent work on the PROFIL programs reported in reference [16] shows that the uncertainties on the calculated isotopic ratio are dominated by the accuracy of the neutron fluence scaling close to Δφ/φ = 1.5% on average. The corresponding uncertainty attached to the average ⟨C/E⟩ values can be estimated with a sensitivity analysis.
The sensitivity S(φ) of the calculated isotopic ratio to the fluence was established by direct perturbation of the 2nd scaling parameter, which normalizes the global fluence level (Sect. 3.3). Values of S(φ) are reported in Table 4 and result from the difference of two ERANOS calculations. The quoted uncertainties on the Neodymium ratios are the product of the sensitivity times the uncertainty associated to the fluence scaling:(8)
This sensitivity study demonstrates that no sizeable dependences of S(φ) are observed between the different isotopic ratios. Therefore, an incorrect fluence scaling cannot explain the increasing bias with the Nd mass number shown in Figures 1 and 5.
4.5 Uncertainties due to the Nd capture crosssections
The final source of uncertainties investigated in this work is the contribution of the uncertainties of the Nd capture crosssections to the calculated isotopic ratio.
Figure 7 compares the Nd capture crosssection available in the JEFF3.1.1 library with a neutron spectrum representative of the PROFIL experiments. For the capture process, the fastenergy range of interest in the neutronic calculations lies between 1 keV and 1 MeV. The sensitivity is maximum around 30 keV. This energy also corresponds to the temperature relevant for the neutron capture nucleosynthesis. Improved compilation of (n,γ) reactions at kT = 30 keV has been updated in order to provide a set of recommended Maxwellian Averaged Neutron Cross Section between 1 keV and 100 keV. Two sets of Nd neutron crosssections published in 2000 and 2006, respectively are reported in Table 5. The recommended values are similar with slightly different uncertainties. We can observe large differences with the JEFF3.1.1 neutron capture crosssections ranging from 4% to 20%, far exceeding the uncertainties recommended in references [19,20]. The bias of 9% and 14% calculated for the ^{143}Nd(n,γ) and ^{145}Nd(n,γ) reactions, respectively, are consistent with the overestimation of 13% and 19% reported in the previous PROFIL interpretations [8,9].
The sensitivity S(σ_{γ}) of the calculated Neodymium buildup in the ^{235}U samples to the Nd capture crosssections are listed in Table 6. They were calculated with the perturbation formalism of the ERANOS code. The mathematical framework is detailed in reference [11]. Low sensitivity coefficients are obtained because the Nd capture crosssections in the fastenergy range remain lower than 1 barn. The negative impact of the ^{A}Nd capture crosssection on the ratio ^{A}Nd/^{235}U is associated to a positive impact of similar magnitude on the ratio ^{A+1}Nd^{235}U. Due to published data not included for the ^{150}Nd(n,γ) reaction at 30 keV in the Atlas of Neutron Resonances [20], the recommendation of Bao et al. [19] were used to estimate the uncertainties reported in the last two columns of Table 6. The uncertainties attached to the Nd capture crosssections and the bias between the JEFF3.1.1 library and the recommended values are treated separately as follows:(9)and(10)
The present results confirm the conclusions reported in reference [22], in which the author suggests that a global uncertainty lower than 1% could have resulted from omitting to make this correction for the Nd capture crosssections.
Fig. 7 Nd capture crosssections available in the JEFF3.1.1 library and recommended in reference [19] at 30 keV compared to a neutron spectrum, in arbitrary units, representative of the PROFIL experiments (red histogram). 
Sensitivity S(σ_{γ}) in %/% of the calculated Neodymium buildup in the ^{235}U samples to the Nd capture crosssections. The last two columns report the relative uncertainties calculated with equations (9) and (10) by using the relative uncertainties reported in reference [19] and ratios given in Table 5.
5 Results and discussions
The experimental results provided by the PROFIL experiments are related to effective cumulated fission yields Y̅_{c}(^{A}Nd), which are associated to a given neutron spectrum. Fastreactor fission yields reported in the literature could have variations due to the neutron energy dependence of the fission yields. The GEF code [12,13] was used to study for this energy dependence and to extract preliminary trends on cumulated fission yields Y_{c}(^{A}Nd,E) at E = 400 keV and E = 500 keV.
5.1 Effective cumulative fission yields
The average calculatedtoexperimental ratios ⟨C/E⟩ listed in Table 2 for PROFIL1, PROFIL2A and PROFIL2B (calculated with JEFF3.1.1) are reported in Table 7 together with the sources of uncertainties discussed in Sections 4.1–4.5. Their quadratic sum leads to overall relative uncertainties ranging from 2% to 2.5%, mainly dominated by the contribution of the fluence scaling uncertainty.
The uncertainty analysis, reported in Section 4, indicates that the main sources of uncertainties cannot alone explain the large drift on the calculatedtoexperimental ratios with respect to the Nd mass number observed in Figures 1 and 5. As a consequence, we can expect that the observed discrepancies between the calculated and experimental values are mainly due to the Nd cumulative fission yields used in the ERANOS calculations. This statement has been verified with a sensitivity analysis.
The sensitivity coefficients S(Y̅_{c}) of the average ratios ⟨C/E⟩ to the effective Nd fission yields were estimated by a direct perturbation analysis. Results reported in Table 7 were calculated from the difference of two ERANOS calculations. They indicate that the S(Y̅_{c}) values are of similar magnitude and close to 0.98, confirming the strong sensitivity of the calculatedtoexperimental ratios to the Nd cumulative fission yields. Since the ERANOS calculations were performed by using the cumulative fission yields of the JEFF3.1.1 library, experimental values Y̅_{c}(^{A}Nd) can be determined as follows:(11)where ΔY̅_{c}(^{A}Nd) is related to the average calculatedtoexperimental ratios ⟨C/E⟩ via the relationship:(12)with the condition:(13)
By introducing equations (12) and (13) in equation (11), we obtain:(14)
In the ERANOS calculations, we have to account for that the cumulative fission yields in the JEFF3.1.1 library are given at three incident neutron energies (25.3 meV, 400 KeV and 14 MeV). In absence of a fine description of the energy dependence of the cumulative fission yields, we always assume that:(15)where Y_{cJ}(^{A}Nd,E) represents the fastcumulative fission yield of JEFF3.1.1, reported in the 3rd column of Table 1. The final experimental values Y̅_{c}(^{A}Nd), calculated with equation (14), are reported in Table 7. The mean values ⟨Y̅_{c}(^{A}Nd)⟩, averaged over the PROFIL1, PROFIL2A and PROFIL2B results, are given below:
The obtained uncertainties were automatically calculated with the uncertainty propagation capabilities of the CONRAD code [21] by using the standard law:(16)in which represents the sensitivity matrix of ⟨Y̅_{c}(^{A}Nd)⟩ to the average ratios ⟨C/E⟩ and D is the covariance matrix between the average ratios ⟨C/E⟩. The covariance matrix D can be separated in two parts to account for the statistical and systematic uncertainties listed in Table 7. The covariance matrix for the statistical part is a diagonal matrix, while the covariance matrix for the systematic part contains correlations close to unity, mainly because of the fluence scaling uncertainty.
Summary of the weighted mean values ⟨C/E⟩, given in Table 2, and of the statistical and systematic uncertainties obtained for the PROFIL1 (1), PROFIL2A (2A) and PROFIL2B (2B) experiments with the JEFF3.1.1 library. The effective fission yields Y̅_{c}(^{A}Nd) were calculated with equation (14) and ⟨Y̅_{c}(^{A}Nd)⟩ are the average values calculated over the three PROFIL experiments.
5.2 Comparison with experimental data
Few sets of data are reported in the literature for fastreactor conditions. The main experimental data of interest were measured in the DOUNREAY [22], EBRI [23], RAPSODIE [24,25], EBRII [26] and PROTEUS [27] fast reactors in the 70s. Figure 8 compares the PROFIL results with a selected set of data used in the evaluation procedure of the fission yields for the JEFF library [14]. For the ^{145}Nd and ^{146}Nd fission yields, our results agree with the data within the limit of the uncertainties. For the ^{148}Nd and ^{150}Nd fission yields, a larger spread between the data can be observed. We have focused our attention on the singular trends of the data reported for ^{148}Nd.
In Figure 8, the 1st values in each plot (open circle) are the experimental fission yield of Davies [22] measured in the DOUNREAY reactor (Caithness, UK). The Davies’ data are systematically higher. For ^{148}Nd, the cumulative fission yield used in the evaluation procedure of JEFF is 0.0175 ± 0.0005. The low uncertainty of 3% has the consequence of increasing its weight on the evaluated yield value.
Robin et al. [24] also report a cumulated fission yield for ^{148}Nd of 0.017 with a low uncertainty of 1.8% at 2σ, in which systematic uncertainties are not included. The Robin's data come from the TACO experiment performed in the fastneutron critical mockup reactor RAPSODIE located in the CEA of Cadarache (France). TACO provided useful technical knowledge for the design of the PROFIL programs in a power reactor. A dedicated analysis of the TACO results were performed by Koch [25] for establishing systematics of fastcumulative fission yields. The analysis of two ^{235}U samples lead to a cumulated fission yield for ^{148}Nd of 0.01665, which is 2% lower than the Robin's value.
If realistic uncertainties are taken into account during the evaluation procedure, a cumulated fission yield for ^{148}Nd close to 0.0168 is expected. This value is fully consistent with the experimental result of 0.0168 ± 0.0002 reported by Maeck in 1975 [26]. Unfortunately, in 1981, Maeck et al. published revised and updated fastreactor fission yields and indicated that these data superseded earlier values from EBRII published in 1975 [28], making questionable any average values that include results coming from EBRII. The discussion of the uncertainties, proposed in this work, demonstrates the difficulty for assessing correct mean values for the evaluated nuclear data libraries.
Fig. 8 Effective Nd cumulative fission yields obtained in this work (Tab. 7) compared to data used in the evaluation procedure of the JEFF3.1.1 Fission Yield library. 
5.3 Comparison with evaluated data at 400 keV and 500 keV
For fastreactor applications, evaluated nuclear data libraries recommend cumulative fission yields at 400 keV (JEFF library) and 500 keV (ENDF/B and JENDL libraries). Those for ^{235}U are reported in Table 1. In Figure 9, our effective values for ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd are compared with values recommended in the JEFF, ENDF/B and JENDL libraries. For JEFF3.1.1, the decreasing trend with the Nd mass number explains the behavior of the calculatedtoexperimental ratios shown in Figures 1 and 5. The cumulative fission yield from the PROFIL experiments are in better agreement with JENDL4.0 and ENDF/BVII.1. A systematic bias close to 2% is observed, meaning that the relative uncertainties of 0.5% quoted in the Japanese and US libraries are underestimated.
In the evaluated libraries, the cumulative fission yields for fast neutrons are given at 400 keV or 500 keV, while any effect of neutron energy on fission yields is accounted for during the evaluation procedure. In reference [26], Maeck indicates that this effect can be significant for different positions in a same reactor, without providing an order of magnitude estimate of this effect for the fission yields measured in the EBRII reactor. In the present work, the variation of fission yields with the neutron energy was investigated with the GEF code. The code is able to provide systematics for fission yields as a function of the energy for a large number of fissile systems. Preliminary results obtained for the ^{235}U(n,f) reaction are shown in Figure 10. The results from GEF were normalized at the thermal energy by using the cumulative fission yields of JEFF3.1.1. The normalization factors range from 1.08 (for ^{143}Nd) to 0.6 (for ^{150}Nd). The comparison with a fastneutron spectrum representative of the PROFIL experiments shows the smooth variation of the fission yields in the energy range of interest for fastreactor applications (around 400 keV and 500 keV). Table 8 reports the normalized cumulative fission yields calculated with GEF at two energies and the effective cumulative fission yields calculated with equation (2). Then, the effect of the energy dependence of the fission yields can be estimated with the ratio:
In the case of the PROFIL experiments, such a neutron spectrum average correction lies between 0.5% and 1.4%. For each Nd isotope, we also observe that the obtained correction varies slowly with the neutron energy, between 400 keV and 500 keV.
The calculated corrections Δ_{ϕ}(E) reported in Table 8 were applied on the effective fission yields ⟨Y̅_{c}(^{A}Nd)⟩ given in Section 5.1 as follow:(18)in order to provide an estimate of Y_{c}(^{A}Nd,E) at E = 400 keV:and at E = 500 keV:
The results are reported in Figure 9 (black circle). We observe a slight reduction of the dispersion of the Nd values in the case of the comparison with the JENDL and ENDF/B fission yields. The variation of the fission yields with the neutron energy could also contribute to explain the PROFIL trends with the Nd mass number obtained with the JEFF library (Figs. 1 and 5). The encouraging results from the GEF code (version 1.7, 2013) are still preliminary and have to be considered with care.
Fig. 9 Comparison of the Nd cumulative fission yields obtained in this work and reported in the JEFF, JENDL and ENDF/B evaluated nuclear data libraries. The open circles represent the effective cumulative fission yields deduced from the PROFIL experiments. The black circles are the cumulative fission yields at 400 keV (for the JEFF library) and 500 keV (for the ENDF/B and JENDL libraries) also deduced from the PROFIL results, but they account for a neutron spectrum average correction calculated with the GEF code (Eq. (17)). 
Fig. 10 Energy dependence of the Nd cumulative fission yields for ^{235}U calculated with the GEF code. The GEF results have been normalized at the thermal energy (25.3 meV) with the thermal fission yield recommended in the JEFF3.1.1 library. They are compared to a neutron spectrum, in arbitrary units, representative of the PROFIL experiments. 
Results calculated with the GEF code after normalization to the thermal fission yields of JEFF3.1.1. Columns (a) and (b) give the cumulative fission yields calculated at 400 keV and 500 keV. The column (c) contains the fission yields from GEF weighted by the uranium fission rate calculated for the ^{235}U samples in the PROFIL experiments. The last two columns give the neutron spectrum average correction defined as the ratios of the fission yields at 400 keV and 500 keV to the effective fission yield.
6 Conclusion
Fast effective Neodymium cumulative fission yields with realistic uncertainties were deduced from the PROFIL experiments. Final results for ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd were 0.0561, 0.0370, 0.0283, 0.0164 and 0.0066 respectively. The comparison of our effective values with the fission yields recommended in the US and Japanese libraries (ENDF/BVII.1 and JENDL4) confirms that the fission yields available in the JEFF3.1.1 library need to be revised considering this new information. We also found a remaining systematic bias of 2% with the US and Japanese libraries. This bias could indicate that the cumulative fission yield for ^{148}Nd is overestimated in the international libraries. Such a result is difficult to confirm with the PROFIL experiments because the final uncertainty attached to each of the effective fission yields lies between 2.1% and 2.4%.
The final uncertainty includes the contribution of the fluence scaling uncertainty, which is close to 2%. This wellknown source of uncertainty mainly depends on the accuracy of the ^{235}U(n,f) reaction. The contribution of the uncertainties on the Nd capture crosssections was also quantified. The magnitude of such a contribution is low (<1%), confirming statements found in previous works. We also report for the first time the impact of the energy dependence of the cumulative fission yields that smoothly decrease in the fastenergy range. The GEF code was used to provide neutron spectrum average corrections ranging from 0.5% to 1.4%. According to these calculations, the neutron spectrum average correction seems to be rather small in the case of the PROFIL experiments.
The GEF calculations have provided encouraging results for future evaluation works. However, this preliminary result has to be confirmed by using an improved version of the code over the seven other fissile systems measured during the PROFIL experiments (^{238}U, ^{238}Pu, ^{239}Pu, ^{240}Pu, ^{241}Pu, ^{242}Pu and ^{241}Am). Such a work is in progress within the frame of the JEFF project.
References
 IAEA, Fission product nuclear data, Technical document IAEA213, 1977 [Google Scholar]
 OECD, Fission product nuclear data, Technical document NEA/NSC/DOC(92)9, 1992 [Google Scholar]
 M.B. Chadwick et al., ENDF/BVII.1 nuclear data for science and technology: cross sections, covariances, fission product yields and decay data, Nucl. Data Sheets 112, 2887 (2011) [CrossRef] [Google Scholar]
 K. Shibata et al., JENDL4.0: a new library for nuclear science and engineering, J. Nucl. Sci. Technol. 48, 1 (2011) [Google Scholar]
 M.A. Kellett et al., The JEFF3.1/3.1.1 radioactive decay data and fission yields sublibraries, Technical document, JEFF report 20, 2009 [Google Scholar]
 G. Rimpault et al., The ERANOS code and data system for fast reactor neutronic analyses, in Proc. Int. Conf. PHYSOR 2002, Seoul, Korea (2002) [Google Scholar]
 J.M. Ruggieri et al., ERANOS 2.1: the international code system for GENIV fast reactor analysis, in Proc. Int. Congress on Advances in Nuclear Power Plants, ICAPP06, Reno, USA (2006) [Google Scholar]
 J. Tommasi et al., Analysis of sample irradiation experiments in PHENIX for JEFF3.0 nuclear data validation, Nucl. Sci. Eng. 154, 119 (2006) [Google Scholar]
 J. Tommasi et al., Analysis of the PROFIL and PROFIL2 sample irradiation experiments in PHENIX for JEFF3.1 nuclear data validation, Nucl. Sci. Eng. 160, 232 (2008) [Google Scholar]
 A. Santamarina et al., The JEFF3.1.1 Nuclear Data Library, Technical document, JEFF report 22, 2009 [Google Scholar]
 A. Gandini, J. Nucl. Energy 21, 755 (1967) [CrossRef] [Google Scholar]
 K.H. Schmidt et al., General Description of Fission Observables, JEFF Report 24, NEA Data Bank, 2014 [Google Scholar]
 K.H. Schmidt et al., General description of fission observables: GEF model code, Nucl. Data Sheets 131, 107 (2016) [CrossRef] [Google Scholar]
 R.W. Mills, Fission product yield evaluation, PhD thesis, University of Birmingham, 1995 [Google Scholar]
 A.R. Date, A.L. Gray, Analyst 108, 159 (1983) [CrossRef] [Google Scholar]
 E. Privas et al., The use of nuclear data as Nuisance parameters in the integral data assimilation of the PROFIL experiments, Nucl. Sci. Eng. 182, 377 (2016) [Google Scholar]
 A.D. Carlson et al., International evaluation of neutron cross section standards, Nucl. Data Sheets 110, 3215 (2009) [CrossRef] [Google Scholar]
 J. Tommasi, private communication, 2007 [Google Scholar]
 Z.Y. Bao et al., Neutron cross sections for nucleosynthesis studies, At. Data Nucl. Data Tables 76, 70 (2000) [Google Scholar]
 S.F. Mughabghab, Atlas of Neutron Resonances , 5th ed. (Elsevier, Amsterdam, 2006) [Google Scholar]
 P. Archier et al., CONRAD evaluation code: development status and perspectives, Nucl. Data Sheets 118, 488 (2014) [CrossRef] [Google Scholar]
 W. Davies, Absolute measurements of fission yields for ^{235}U and ^{239}Pu in the Dounreay fast reactor, Radio. Acta 12, 174 (1969) [Google Scholar]
 F.L. Lisman et al., Fission yields of over 40 stable and longlived fission products for thermal neutron fissioned ^{233}U, ^{235}U, ^{239}Pu and ^{241}Pu and fast reactor fissioned ^{235}U and ^{239}Pu, Nucl. Sci. Eng. 42, 191 (1970) [CrossRef] [Google Scholar]
 M. Robin et al., The importance of fission product nuclear data in burnup determination, in Proc. Int. Conf. on Chemical and Nuclear Data, Canterbury, UK (1973) [Google Scholar]
 L. Koch, Systematics of fast cumulative fission yields, Radio. Acta 29, 61 (1981) [Google Scholar]
 W.J. Maeck, Fast reactor fission yields for ^{233}U, ^{235}U, ^{238}U and ^{239}Pu, and recommendations for the determination of burnup on FBR mixed oxide fuels, in Proc. Int. Conf. on Nuclear Cross Sections and Technology, Washington, USA (1975) [Google Scholar]
 M. Rajagopalan et al., Mass yields in the fission of Uranium235 and Plutonium239 in the neutron spectrum of a gascooled fast reactor, Nucl. Sci. Eng. 58, 414 (1975) [CrossRef] [Google Scholar]
 W.J. Maeck et al., Revised EBRII fastreactor fission yields for ^{233}U, ^{235}U and ^{238}U, Exxon Nuclear Idaho Company, report ENICO1091, 1981 [Google Scholar]
Cite this article as: Edwin Privas, Gilles Noguere, Jean Tommasi, Cyrille De Saint Jean, KarlHeinz Schmidt, Robert Mills, Measurements of the effective cumulative fission yields of ^{143}Nd, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd for ^{235}U in the PHENIX fast reactor, EPJ Nuclear Sci. Technol. 2, 32 (2016)
All Tables
C/E ratios for the prediction of the Neodymium buildup, after fluence scaling. Results obtained in the present work with the JEFF3.1.1 library are compared with those obtained with the JEFF3.1 library [18]. The 1st and 2nd weighted mean values for PROFIL were calculated with and without the C/E results obtained for the ^{235}U sample number 19 located in the core midplane (z = 0 mm). The quoted uncertainties are the experimental uncertainties coming from the chemical analysis of the ^{235}U samples.
Sensitivity S(σ_{γ}) in %/% of the calculated Neodymium buildup in the ^{235}U samples to the Nd capture crosssections. The last two columns report the relative uncertainties calculated with equations (9) and (10) by using the relative uncertainties reported in reference [19] and ratios given in Table 5.
Summary of the weighted mean values ⟨C/E⟩, given in Table 2, and of the statistical and systematic uncertainties obtained for the PROFIL1 (1), PROFIL2A (2A) and PROFIL2B (2B) experiments with the JEFF3.1.1 library. The effective fission yields Y̅_{c}(^{A}Nd) were calculated with equation (14) and ⟨Y̅_{c}(^{A}Nd)⟩ are the average values calculated over the three PROFIL experiments.
Results calculated with the GEF code after normalization to the thermal fission yields of JEFF3.1.1. Columns (a) and (b) give the cumulative fission yields calculated at 400 keV and 500 keV. The column (c) contains the fission yields from GEF weighted by the uranium fission rate calculated for the ^{235}U samples in the PROFIL experiments. The last two columns give the neutron spectrum average correction defined as the ratios of the fission yields at 400 keV and 500 keV to the effective fission yield.
All Figures
Fig. 1 Average C/E ratios for the prediction of the Neodymium buildup in the PROFIL1 and PROFIL2 programs for ^{235}U, ^{239}Pu and ^{241}Pu obtained with the ERANOS code and the JEFF3.1 library [9]. 

In the text 
Fig. 2 ^{235}U capture and fission crosssections available in the JEFF3.1.1 library compared to a neutron spectrum, in arbitrary unit, representative of the PROFIL experiment (red histogram). 

In the text 
Fig. 3 Radial core geometry of the fast reactor PHENIX, location of the subassemblies containing the PROFIL pins, axial layout of the PROFIL1 pin and example of a stainless steel double container designed for the PROFIL irradiation. 

In the text 
Fig. 4 Axial layout of the PROFIL1 pin with the position of the ^{235}U samples used in this work. The sample number 19 is located at z = 0 mm, which corresponds to the core midplane of the PHENIX reactor. 

In the text 
Fig. 5 Average ⟨C/E⟩ ratios for the prediction of the Neodymium buildup in the PROFIL1, PROFIL2A and PROFIL2B programs for ^{235}U obtained with the ERANOS code and the JEFF3.1.1 library. Individual results are listed in Table 2. 

In the text 
Fig. 6 Comparison of the individual C/E results obtained with JEFF3.1.1 (open circle) and JEFF3.1 (black circle) for the prediction of the ^{143}Nd, ^{144}Nd + ^{144}Ce, ^{145}Nd, ^{146}Nd, ^{148}Nd and ^{150}Nd buildup in the ^{235}U samples of the PROFIL1 experiment as a function of the sample position. The abscise z = 0 mm corresponds to the core midplane of the PHENIX reactor. The solid line is a parabolic curve fitted to the data acting as a eye guide. 

In the text 
Fig. 7 Nd capture crosssections available in the JEFF3.1.1 library and recommended in reference [19] at 30 keV compared to a neutron spectrum, in arbitrary units, representative of the PROFIL experiments (red histogram). 

In the text 
Fig. 8 Effective Nd cumulative fission yields obtained in this work (Tab. 7) compared to data used in the evaluation procedure of the JEFF3.1.1 Fission Yield library. 

In the text 
Fig. 9 Comparison of the Nd cumulative fission yields obtained in this work and reported in the JEFF, JENDL and ENDF/B evaluated nuclear data libraries. The open circles represent the effective cumulative fission yields deduced from the PROFIL experiments. The black circles are the cumulative fission yields at 400 keV (for the JEFF library) and 500 keV (for the ENDF/B and JENDL libraries) also deduced from the PROFIL results, but they account for a neutron spectrum average correction calculated with the GEF code (Eq. (17)). 

In the text 
Fig. 10 Energy dependence of the Nd cumulative fission yields for ^{235}U calculated with the GEF code. The GEF results have been normalized at the thermal energy (25.3 meV) with the thermal fission yield recommended in the JEFF3.1.1 library. They are compared to a neutron spectrum, in arbitrary units, representative of the PROFIL experiments. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.