Issue 
EPJ Nuclear Sci. Technol.
Volume 4, 2018
Special Issue on 4th International Workshop on Nuclear Data Covariances, October 2–6, 2017, Aix en Provence, France – CW2017



Article Number  17  
Number of page(s)  5  
Section  Nuclear Data Library  
DOI  https://doi.org/10.1051/epjn/2018039  
Published online  05 September 2018 
https://doi.org/10.1051/epjn/2018039
Regular Article
Resonance parameter covariance representation: file32 versus file33
Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSNEXP, SNC,
92260
FontenayauxRoses France
^{*} email: luiz.leal@irsn.fr
Received:
1
December
2017
Received in final form:
25
January
2018
Accepted:
28
May
2018
Published online: 5 September 2018
Sensitivity and uncertainty analysis in error propagation studies are carried out based on nuclear data uncertainty information available in the basic nuclear data libraries such as ENDF, JEFF, JENDL and others. The uncertainty files (covariance matrices) are generally obtained from analysis of experimental data. In the resonance region, the computer code SAMMY is used for analyses of experimental data and generation of resonance parameters. In addition to resonance parameters evaluation, SAMMY also generates resonance parameter covariance matrices (RPCM). The intent of this paper is to discuss the use of the RPCM in contrast to the use of a cross section covariance matrix (CSCM) representation. For so, the resonance evaluation for ^{48}Ti in the resonance range from 10^{−5} eV to 400 keV is used. The RPCM is translated into two distinct CSCM by using a coarser and a finer group structures. The results are presented and discussions of the feasibility of using these covariance representations are depicted.
© L. Leal, published by EDP Sciences, 2018
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 objective of a resolved resonance parameters evaluation is to generate parameters which combined with a resonance formalism, the Rmatrix methodology for instance, lead to a reasonable representation of the experimental data such as the total, capture, fission cross sections, and etc. Along with the resonance parameters (RP), resonance parameter covariance matrices (RPCM) are also obtained. For practical applications the RP parameters are converted to the evaluated nuclear data file format (ENDF) in the socalled FILE2 representations whereas the RPCM is converted to the ENDF FILE32 representation. While the FILE32 representation is preferred, however situations may arise whereby the use of the cross section covariance matrix (CSCM) is recommended which consists of translating FILE32 into FILE33. For instance, if the RP evaluation results in a large set of RP a drawback may be connected to computer's disk storage of the RPCM. Nevertheless, the process of converting to FILE33 is carried out based on two assumptions: a) the CSCM are provided in the energymultigroup form; b) a weighting spectrum is needed for the crosssection averaging. The CSCM provides the variances of the cross section within a specified energy region and also the correlation between cross sections of nearby energy regions. The objective of this paper is to examine the use of the CSSM representation in lieu of the RPCM. The studies are performed by way of using the RPCM generated on the basis of a resolved resonance evaluation carried out for ^{48}Ti. A brief description of the resonance of ^{48}Ti will be introduced prior to the discussion on the use of the CSCM. The paper will examine whether the FILE32 and FILE33 conversion grants an alternative leading to comparable calculated group covariance. The adequacy of energy group choice for the CSCM representation in the nuclear data library will be investigated.
2 ^{48}Ti resolved resonance evaluation
In natural titanium the most abundant isotope is ^{48}Ti which has an abundance of 73.72% in weight. Additionally, ^{48}Ti has the largest thermal capture cross section among the natural titanium isotopes. In response to needs for criticality safety applications, revision and evaluation of nuclear systems in which titanium is present were requested. It was identified that existing titanium evaluation should be revised and a new evaluation should be put in place. Hence, an Rmatrix resonance evaluation based on the code SAMMY [1] was performed [2] for ^{48}Ti in the energy range from thermal to 400 keV. In addition, the uncertainty in the cross sections were derived from resonance parameter covariance obtained in the evaluation. Since no reaction other than total, scattering and capture exists below 400 keV, the resonance parameters are uniquely determined in the evaluation by the energy of the resonance, the neutron width, the capture width, and the resonance spin. Correspondingly to these parameters are the associated variances as well as their correlations. For instance, for two resonance parameters p_{j} and p_{k}, which can be energy or widths, their correlations leading to the RPCM is referred to as ⟨δp_{j}δp_{k}⟩. More o.0lication of the RPCM will be seen in section 3. As an example of the use the resonance parameter to reproduce the experimental cross section, the fitting of the differential total and capture cross sections in the energy range of 50 keV to 150 keV is displayed in Figure 1.
The dots are the experimental data whereas the solid line going through the data represents the evaluation. The fitting of the experimental data was done using the ReichMoore formalism included in the code SAMMY in conjunction with the Bayes' methodology. In addition to the crosssection fitting the uncertainty in the nuclear data can also be extracted from the evaluation process in which resonance parameter covariance are obtained. The data uncertainties are statistical and systematic. The latter is the predominant data uncertainty. The uncertainties information on the experimental data used in the resonance evaluation are as follow: The observed total and capture cross section systematic uncertainties are in between 1–4% and 2–10%, respectively. These uncertainties are reflected in the final RPCM and consequently in the cross sections. Figure 2 illustrates the uncertainties and covariance for the ^{48}Ti evaluated capture cross section in the energy region from 10^{−5} eV to 20 MeV. The covariance information is used in the analysis and design of nuclear systems. The SAMMY generated ^{48}Ti RPCM were converted into the ENDF format [4] on the basis of the LCOMP = 1 option. This option requires that the entire covariance data are given for the full RPCM covering the resolved resonance region. Thermal values and capture resonance integral for the ^{48}Ti evaluation are compared with that included in the ENDF/BVII.0. The results shown in Table 1 include the uncertainty in the cross sections that were calculated from the covariance derived from the resonance evaluation as indicated in the example shown in Figure 2. It should be pointed out that the most recent releases of ENDF (ENDF/BVIII.beta5) and JEFF (JEFF3.3T4) include the new ^{48}Ti resonance parameter evaluation.
Fig. 1 SAMMY fit of the experimental differential total and capture cross sections in the energy range 50 keV–150 keV [3]. 
Fig. 2 Uncertainty and correlation in the Ti capture cross section using 44energy groups. 
Cross section at thermal energy (0.0253 eV) and capture resonance integral.
3 Covariance generation
3.1 File32 into File33
For an energygroup with energy boundary E_{g} and E_{g}_{+1} the group cross section covariance for a reaction x as a function of the resonance parameter covariance is given as (1)
The covariance of , the fluxweighted group cross section for reaction x, is a function of the derivative of the cross sections with respect to the parameters p_{j} and p_{k}, that is, the cross section sensitivity to the RP, and the RPCM given as ⟨δp_{j}δp_{k}⟩ which are quantities stored in FILE32. On the other hand, if the CSCM listed in FILE33 is M(σ), the is obtained as (2) where (3)
Here φ_{g} is the energydependent neutron flux in the energy group g. The results of equations (1) and (2) are equivalent as long as the conversion of ⟨δp_{j}δp_{k}⟩ into M(σ) is properly done [5]. The FILE33 CSCM representation is intended to characterize the variances of the cross sections within a specified energy region, and the correlation between cross sections of adjacent energy regions. The choice of CSCM over the RPCM is expected to lead to a reduction in computer storage. However, a gain in computer use, meaning data storage and perhaps computer time, may not preserve the accuracy sought in practical applications such as, for instance, the establishing of criticality safety limits and their related uncertainties. The issues, pros and cons, in connection with the conversion from FILE32 to FILE33 representation is investigated in the following sections.
3.3 File32 to File33 Conversion for ^{48}Ti
The conversion from the RPCM to the CSCM representation has been investigated by means of using two energy group configurations, a coarser and a finer energy group structures. The finer group structure consists of an energy grid of 381 boundaries, i.e. 380 energy groups, whereas the coarse energy group is rather composed of 44energy groups chosen as a subset of the 380group structure. A graphically representation of these group structures is displayed in Figure 3. While the 380group representation is very detailed approaching to a single line in the picture (black line) the 44group (red line) clearly shows the effect of the coarse energy mesh. The conversions were carried out with the AMPX code system [6] using the modules, PRELL, PRILOSEC, PUFF, and COVCONV. The detailed information on how to use these AMPX modules can be found in the AMPX manual. A constant energydependent neutron flux, required for weighing the covariance as indicated in equations (2) and (3), were used in the conversion. It should be recalled that the results presented here are for the energy region 10^{−5} eV–400 keV corresponding to the resolved resonance evaluation of ^{48}Ti. For the 44group structure 35energy groups fall under 400 keV whereas 343energy groups have energies under 400 keV for the 380group structure. The exercise led, respectively, to two ENDF formatted crosssection libraries including the 44group and the 380group covariance generated in the conversion of the RPCM to the CSCM representation. The data libraries used in the studies presented in this work are identified as follows: a) the covariance library based solely on the resonance parameter representation will be referred to as the RPCM_LIB; b) the covariance library resulting from the conversion from the RPCM to the CSCM using the 44group structure is referred to as the CSCM_44LIB; c) lastly, the 380group library is named CSCM_380LIB. For the sake of information the computer storage size of the covariance portion of these libraries are listed in Table 2.
Fig. 3 Comparison of the 380 group (black line) and 44 group (red line) structures. 
Computer storage size for covariance libraries.
3.4 Tests and results
For practical applications the covariance information in the nuclear data libraries such as ENDF, JEFF, JENDL and others must be processed in a form suitable for use in sensitivity uncertainty (S/U) analysis. There are several tools that can make use of the processed libraries for error propagation studies. The most common processed form of the covariance data these codes use can take either the COVERX [7] or the BOXER format. The S/U code may have been built to accept both forms. In the study carried out here both forms of the processed covariance were generated.
The impact of the covariance representation in the resonance region of the ^{48}Ti evaluation was assessed by processing the libraries RPCM_LIB, CSCM_44LIB, and CSCM_380LIB on the basis of two energy group structures. The energy groups chosen are those listed in the SCALE code system, [6] namely the 27group and the 238group. The results were obtained by processing the covariance information in the libraries using the ERRORR module of the NJOY code system [8] with a constant energydependent neutron flux. The reported results are for the ^{48}Ti total and capture cross sections. The results are displayed in the graphical form. For each group structure (27group and 238group) three calculations were carried out corresponding to the three libraries. The results are overlaid for a better visualization of the impact of using the three covariance representations. The results for the total and capture cross sections are shown in Figures 4 and 5, respectively. The upper portion of each picture represents the corresponding cross section whereas the bottom part is the calculated uncertainty. It appears that for the 27group generated covariance the CSCM representation underestimate the uncertainty compared to the RPCM for the total cross section in the energy range where resonances exist. The same trend is also observed for the capture cross section but nevertheless the difference where the resonances start are very small. The results based on the 238goup are shown in Figures 6 and 7 , respectively. The two results for the total and the capture cross sections indicate that in the energy region where no resonance is present apparently the 44group covariance representation seems to be represented reasonably well. However, as the resonances arise the results based on the 27group representation deteriorates showing large differences in comparison with the RPCM calculations. On the other hand the CSCM based on the 380group is in very good agreement with the RPCM. Nonetheless the improvement on the results based on the 380group are based on a covariance matrix that requires a large computer storage size as indicated in Table 2 and therefore not offering any advantage over the resonance parameter representation. The results of these calculations lead to some inferences:
care must be taken when using fewgroup covariance representations. It may be appropriate for the energy range where the data are smooth and the resonances are not present as for instance in the high energy region;
for a detailed covariance results the use of the RPCM is recommended;
conversion of the RPCM to CSCM using a fine energygroup structure must be examined so as to assure that the computer allocation size is not overwhelming larger than that of the RPCM;
in the resonance region fewgroup representation is acceptable when a general overview of the impact of covariance results is sought.
Fig. 4 27group calculated covariance for the total cross section. 
Fig. 5 27group calculated covariance for the capture cross section. 
Fig. 6 238group calculated covariance for the total cross section. 
Fig. 7 238group calculated covariance for the capture cross section. 
4 Conclusions
This paper presents, at a certain extent, a discussion on the representation of the covariance data in the evaluated nuclear data libraries. The resonance parameter covariance for ^{48}Ti was used in the study. The impact of using the RPCM as opposed to the CSCM representation and the results of calculations based on these covariance representations are presented. Three covariance representations were used in the study, namely, one based on the resonance parameters, a coarser group covariance using 44group and a finer group covariance using 380energy groups. The results suggested that the use of the RPCM is preferred over the CSCM representation. However there are situation in which a finegroup group structure may be advised as for instance in the case where the RPCM lends itself to a huge covariance matrix representation. Another alternative could be to reduce the size of the RPCM by using a compact representation of the covariance or another preferred approach.
There are other issues that were not addressed in this paper such as the impact of the weighting spectrum in the covariance conversion from RPCM to CSCM. While this issue may impact the results, the conversion of the resonance parameter representation to cross section representation of the covariance together with the choice of the energy group structure is believed to be a more dominant subject. However, work is underway to understand the impact of the covariance weighting in the conversion process. Other subjects under examination are the impact of the cross section energy selfshielding and the temperature effect in the covariance. The results presented in this work were obtained for room temperature and infinite diluted cross sections. Also, no attempt has been made in this paper to examine the implication of the RPCM and CSCM representations in practical applications such as the magnitude of the propagated uncertainty to k_{eff} for system sensitive to the ^{48}Ti.
References
 N.M. Larson, Updated Users Guide for SAMMY: MultiLevel RMatrix Fits to Neutron Data Using Bayes' Equations, (ENDF364/R2, Oak Ridge National Laboratory, 2008), available at the Radiation Safety Information Computational Center (RSICC) as PSR158 [Google Scholar]
 L. Leal, K. Klaus, G. Arbanas, D. Wiarda, P. Koehler, A. Kahler, Resonance evaluation of Titanium48 including covariance for criticality safety applications, in International Conference on Nuclear Criticality Safety, Edinburgh, United Kingdom, 2011 [Google Scholar]
 K.H. Guber, P.E. Koehler, D. Wiarda, J.A. Harvey, Neutron CrossSection Measurements on Structural Materials at ORELA, in Proceedings of the International Conference on Nuclear Data for Science and Technology, Jeju Island, Korea , 2010 [Google Scholar]
 ENDF6 Formats Manual, Formats and Procedures for the Evaluated Nuclear Data Files ENDF/BVI and ENDF/BVII, edited by M. Herman, A. Trkov, ReportBNL903652009 Rev. 1 (National Nuclear Data Center, Upton, NY, 2010) [Google Scholar]
 L. Leal, D. Mueller, G. Arbanas, D. Wiarda, H. Derrien, Impact of the ^{235}U Covariance Data in Benchmark Calculations, in International Conference on the Physics of Reactors Nuclear Power: a sustainable resource CasinoKursaal Conference Center, Interlaken, Switzerland, 2008 [Google Scholar]
 B.T. Rearden, M.A. Jesse, et al., SCALE Code System,' ORNL/TM2005/39 Version 6.2.1, 2016 [Google Scholar]
 J.D. Drischler, The COVERX Service Module of the FORSS System, ORNL/TM7181, Oak Ridge National Laboratory, Oak Ridge, Tenn, 1980 [Google Scholar]
 R.E. MacFarlane, D.W. Muir, A.C. Kahler, The NJOY Nuclear Data Processing System, Version 2012, LAUR12 27079, Los Alamos National Laboratory, 2012 [Google Scholar]
Cite this article as: Luiz Leal, Resonance parameter covariance representation: file32 versus file33, EPJ Nuclear Sci. Technol. 4, 17 (2018)
All Tables
All Figures
Fig. 1 SAMMY fit of the experimental differential total and capture cross sections in the energy range 50 keV–150 keV [3]. 

In the text 
Fig. 2 Uncertainty and correlation in the Ti capture cross section using 44energy groups. 

In the text 
Fig. 3 Comparison of the 380 group (black line) and 44 group (red line) structures. 

In the text 
Fig. 4 27group calculated covariance for the total cross section. 

In the text 
Fig. 5 27group calculated covariance for the capture cross section. 

In the text 
Fig. 6 238group calculated covariance for the total cross section. 

In the text 
Fig. 7 238group calculated covariance for the capture cross section. 

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.