Issue 
EPJ Nuclear Sci. Technol.
Volume 8, 2022



Article Number  3  
Number of page(s)  7  
DOI  https://doi.org/10.1051/epjn/2021028  
Published online  17 January 2022 
https://doi.org/10.1051/epjn/2021028
Regular Article
Impact of H in H_{2}O thermal scattering data on criticality calculation: uncertainty and adjustment
^{1}
Reactor Physics and Thermal hydraulic Laboratory, Paul Scherrer Institut,
Villigen, Switzerland
^{2}
International Atomic Energy Agency, Vienna International Centre,
PO Box 100,
1400
Vienna, Austria
^{*} email: dimitrialexandre.rochman@psi.ch
Received:
20
September
2021
Received in final form:
29
November
2021
Accepted:
1
December
2021
Published online: 17 January 2022
In this paper, the impact of the thermal scattering data for H in H_{2}0 is estimated on criticality benchmarks, based on the variations of the CAB model parameters. The Total Monte Carlo method for uncertainty propagation is applied for 63 k_{eff} criticality cases, sensitive to H in H_{2}0. It is found that their impact is of a few tenth of pcm, up to 300 pcm maximum, and showing highly nonlinear distributions. In a second step, an adjustment is proposed for these thermal scattering data, leading to a better agreement between calculated and experimental k_{eff} values, following an increase of scattering contribution. This work falls into the global approach of combining advanced theoretical modelling of nuclear data, followed by possible adjustment in order to improve the performances of a nuclear data library.
© D. Rochman et al., Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://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 evaluation of neutroninduced nuclear data, e.g. cross sections, fission yields, spectra, is a continuous process, mainly because the nuclear reaction theory is still improving, new measurements are frequently available, and because there is a need from the user community. Examples of recently evaluated nuclear data are the Nuclear Energy Agency JEFF3.3 library [1], the US ENDF/BVIII.0 library [2], the Japanese JENDL4.0 library [3] and finally the TENDL libraries [4]. The need from users might concern better estimations of current reactor quantities (boron letdown curves, power maps) [5–7], spent fuel quantities (decay heat, source terms) [8–12], or quantities for advanced systems [13–15]. A commonlyused neutronics parameter is the neutron multiplication factor, or k_{eff}, which plays for instance a crucial role in spent fuel management and in burnup credit [16–18]. The validation of existing simulation codes and the assessment of the uncertainties from various sources is a necessity, often imposed by national and international regulations (see Refs. [19–21] and references within). In this context, the impact of the nuclear data uncertainties needs to be evaluated in a regular time manner, with respect to the release of new nuclear data libraries and their associated covariance information.
This is precisely the first goal of this paper: providing uncertainties for a number of criticalitysafety benchmarks [22], due to the H in H_{2}O thermal scattering data (also called Thermal Scattering Laws, or TSL). These TSL are from the latest evaluations based on the CAB model [23], and included in the ENDF/BVIII.0 and TENDL2021 libraries. Regarding their uncertainties, reference [24] is used as a reference with a number of modifications (see Sect. 2). The uncertainty propagation method is the Total Monte Carlo method [25], leading to k_{eff} uncertainties (later called Δk_{eff}), corresponding to one standard deviation of the observed probability density functions (pdf), as well as higher moments. Similar works on uncertainty propagation and covariance generation for TSL can be found in references [26,27]. In the first reference, the TMC method (and fast TMC) was applied, whereas the second reference defines TSL covariance matrices, with an application to UO_{2} and MOX fuels.
The second goal of this paper is to propose adjusted TSL, or at least to demonstrate the feasibility of the adjustment, based on the result of the TMC method. Such possibility was already explored almost ten years ago and presented in reference [28] for the TSL based on a different evaluation model. The method, called the Petten method, simply proposes the best performing sampled (random) TSL, based on the uncertainty propagation for criticality benchmarks, as the new recommended evaluation. Such method was also applied to Cu and ^{239}Pu [29,30]. In this approach, the best performing sampled TSL is the one leading to the smallest χ^{2}. Other methods of adjustments have been proposed since, apart from the traditional GLLS method (Generalized Linear Least Squares), based on nonlinear Monte Carlo parameter search, see for instance references [31–36], the Unified Monte Carlo method [37] and references [38,39] for a comparison and review between GLLS and the Monte Carlobased methods. Naturally, based on the first part of this study (uncertainty propagation), we are calculating χ^{2} values for each sampled TSL, and then proposing the best performing one as an example of adjustment.
Details of the study are provided in the next sections: first a recall on the sampling method of the TSL, followed by uncertainty values for a number of criticality benchmarks. Additionally, other moments are obtained, such as the skewness and kurtosis, indicating a nonlinear propagation of uncertainties from the TSL to k_{eff}. Finally, the last section will analyze the χ^{2} distribution and draw conclusions and a potential adjustment of the TSL.
2 Variation of thermal scattering data
The same approach and sampled TSL files were presented and used in reference [40]. It is basically based on the sampling of the parameters of the CAB model [23]. The nominal (not perturbed) parameters of this model can be found in reference [41] and as mentioned, such TSL were included in some of the modern nuclear data libraries. The variations of the model parameters are derived from reference [42], with a number of modification, following the recommendations of the authors. Such variations are as follows:
σ_{s} (elastic cross section): ± 0.2%
Δ (scaling factor): ± 10%
C (diffusion coefficient): ± 0.5%
ω_{t} (translational weight): ± 15%
ω_{c} (continuous spectrum weight), ω_{1} and ω_{2} (1^{st} and 2^{nd} oscillator weights: correlated with ω_{t} based on equation (1)
E_{1} (first oscillator energy): ± 5%, and
E_{2} (second oscillator energy): ± 30%.
For the oscillator weights ω_{c}, ω_{1} and ω_{2}, a full correlation is considered with the following equation: (1)
Such approach follows the recommendation of reference [24]. All variations are a Normal distributions. Uncertainties correspond to one standard deviation (1σ). As in the case of reference [40], 3000 sampled TSL are produced, and processed with NJOY into ACE files [43].
3 Results
The impact of the H in H_{2}O thermal scattering data are here assessed on criticalitysafety benchmarks. Such benchmarks are obtained from the ICSBEP handbook [22]. Such handbook contains a large number of benchmarks (simulation input files and experimental quantities such as k_{eff}), but only some of them are sensitive to TSL. Those benchmarks are mostly thermal ones (the neutron fluence is mainly in the thermal energy region), but not only. In the following, a total of 63 benchmarks are considered, among which only one is considered as fast (called pmf111) and one intermediate (usi11); the other ones are thermal (61 benchmarks). The nomenclature used here is similar to the handbook: all benchmark names consist of three letters: the first one denotes the type of fuel (“l” for lowenriched ^{235}U, “h” for highenriched ^{235}U, “u” for ^{233}U fuel, “m” for a mixture of ^{235}U and ^{239}Pu, and “p” for ^{239}Pu fuel), the second one defines the physical form (“s” for solution and “c” for compound), and finally the last letter defines the neutron spectrum (“t” for thermal, “i” for intermediate and “f” for fast). All benchmark names are continued with an additional information (after a dash), indicating its subnumber, referring to the ICSBEP cases.
The benchmark calculations are performed with MCNP6; apart from the sampled TSL ACE files, the other nuclear data files correspond to the beta version of the JEFF4.0 library. The method of uncertainty propagation is the Total Monte Carlo (or TMC) [25]; it was already presented in many publications and it is summarized here as follows: it consists in repeating i times the same calculation, each time with a different sampled quantity of interest as input. In other words, each random file differs slightly from the previous one and averaging over all of them would return a credible standard deviation and covariance matrix. The averaging however takes place after the applied calculations. In our case, the quantity of interest corresponds to the TSL data, as explained in Section 2. In this work, i is equal to 3000 and it allows to achieve an acceptable convergence for various moments of the observed k_{eff} probability density functions. In all uncertainty values, the statistical uncertainties from the MCNP6 is already removed, assuming that the observed variance is the sum of the effect of the TSL and statistics. Finally, together with the uncertainty Δ k_{eff} (or one standard deviation), the quantity MAD (for Median Absolute Deviation) is also provided. This is another measure of the width of a distribution, less affected by extreme values (both the standard deviation and MAD values will significantly differ for skewed distributions; for a Normal distribution, the MAD is equal to 0.67 of the standard deviation). The definition of the MAD is Median(x_{i}), with x_{i} the observed realization for the sample i, and the average over all i. Results for the average, uncertainty, skewness and kurtosis for each considered benchmark are presented in Table 1.
Uncertainties (standard deviation and MAD), skewness and kurtosis on criticality benchmarks due to thermal scattering data of H in H_{2}0 obtained from MCNP6 calculations. The base library (except the thermal scattering data) is JEFF4.0β. The statistical uncertainties are already removed.
3.1 Uncertainties and other moments
Compared to uncertainties calculated in reference [28], the present values are in general smaller. This is not surprising as the model parameters in reference [28] were modified in an adhoc way, being one of the first trials of the TMC method on TSL. As observed, the maximum uncertainties are less than 300 pcm (and in general less than 100 pcm). It is worth noticing that for skewed distributions with a strong kurtosis, as the majority of cases presented here, one standard deviation does not refer to 68% of the population: for lct8.5, 68.3% of the cases are within the mean ± 1σ, but for lst2.1, 73% are within the mean ± 1σ. One can also compare the present results for the lct and pst benchmarks with the values presented in reference [27]. In this reference, uncertainties due to TSL for the MISTRAL1 experiment, similar to the lct benchmarks, range from 70 to 150 pcm. Such uncertainties were obtained by a direct perturbation approach, based on different model uncertainties, and are of the same order of magnitude compared to the present values. In the case of the MOX cases, reference [27] reports uncertainties between 110 and 200 pcm for the MISTRAL2 system, which can be compared to the present values for the pst benchmarks. In this case too, uncertainties are relatively similar.
Generally, the uncertainties (as well as skewness and kurtosis) are similar within types of benchmarks: for instance all lct benchmarks have negative skewnesses, all lct7 cases present uncertainties close to the same values. Examples of some pdf are presented in Figure 1. These distributions are sensibly different, with strong different skewnesses (for the pst and uct cases), zero skewness (the lct case), or zero kurtosis (corresponding to a flat distribution, for mct2). Given that the model parameters for the TSL production were sampled from Normal distributions, such a variety of k_{eff} pdf indicates non linear cases of uncertainty propagation.
Similarities in pdf within families of benchmarks are certainly related to the fact that these benchmarks share similar geometry, enrichment, and possibly reflectors. This can be observed for the majority of cases, but with a number of exceptions, as for the uct1 benchmarks and the cases 3 and 7: uncertainties strongly differ, and skewnesses are far from zero with opposite signs. Both benchmarks are also strongly anticorrelated, with a Pearson correlation factor of −0.8. These two cases nevertheless present very different physical systems (e.g. rectangular and hexagonal cores with different reflectors).
The values from the lct benchmarks (being clusters of UO_{2} rods in water) can be compared to the k_{∞} uncertainties from reference [40] (being for a UO_{2} pincell, 3.5% enriched, burned up to 100 MWd/kgU; same sampled TSL files were used), where the maximum uncertainty was about 30 pcm. In the case of the lct benchmarks, the maximum is about 80 pcm. These two values differ, but are also relatively close, given that a direct comparison is not possible (different k_{eff} values, different enrichments, and different geometry).
Based on these random k_{eff} values (3000for each benchmark), the adjustment step can be performed, as presented in the next section.
Fig. 1 Probability density function for a number of criticality benchmarks, obtained by randomly changing the H in H_{2}O thermal scattering data. 
3.2 Adjustment
The previous results were obtained based on 3000 sampled TSL; they all differ and lead to different benchmark k_{eff} values. Let’s i denotes a specific sampled TSL, and n the benchmark case (n goes from 1 to N = 63). For each i, a χ^{2} value can be calculated as follows: (2)
k_{i,n,calc} and k_{n,exp} correspond to the calculated and experimental k_{eff} values for the benchmark n, respectively. For instance, for i = 0, which corresponds to the beta version of the future JEFF4.0 library, called JEFF4T0, . By considering all i sampled cases, the obtained χ^{2} distribution can be calculated and is presented in Figure 2. The shape of this distribution is typical of χ^{2} values, as seen in references [28–30]. As it can be easily understood, some sampled TSL do not perform well (high χ^{2}), whereas some of them are performing better than the original evaluated TSL. Based on such distribution, the Petten method is excessively simple: select the sampled TSL with the lowest χ^{2}. Naturally, as the TSL have a limited effect on the benchmark calculated values, the low tail of the χ^{2} distribution does not reach χ^{2} = 0, but as observed, noticeable improvements in terms of C/E are possible. Out of the 3000 sampled TSL, the one with the smallest χ^{2} is number 1871. The calculated k_{eff} from this specific TSL are presented in Figure 3, compared with the nominal (unperturbed) TSL evaluation. The changes on k_{eff} values from the nominal calculations and from the TSL file number 1871 are not that striking, but the χ^{2} values were definitively reduced. One has to keep in mind that the TSL induces low uncertainties on k_{eff}, meaning that their impact is nevertheless limited. One can also repeat this work and double the uncertainties on the CAB parameters; this way, it is likely to find an even better solution. Naturally, this would imply that the TSL are “responsible” for more than they should be, given that other nuclear data can also be adjusted. Users who are not able or prepared to do the full uncertainty quantification with all 3000 random files, can then at least use random file number 1871 to get the best result.
This specific TSL run number 1871 was obtained with different CAB model parameters, see Table 2 (parameters ω_{c}, ω_{1} and ω_{2} are fully correlated with ω_{t} using Eq. (1)). The parameters for the file number 1871 are one example of adjustments with integral benchmarks. Table 2 is also showing parameters for random files performing almost as good as number 1871. As observed, their parameters are sensibly different, and it is difficult to extract a pattern. The E_{2} for run 1871 might not be suitable, as such value is not valid for the second oscillator, which was closer to 400 meV. To avoid such result, a rejection condition could be used, based on water physics. This indicates that the adjustment solely based on this number of integral benchmarks might not be enough to obtain a unique solution, either because there is not a unique solution, or because the number of benchmarks is not enough. Additionally, the χ^{2} definition can take into account differential data (measured cross sections), thus leading to a compromise between integral and differential agreements. More advanced method can also be applied, based on Bayesian Monte Carlo, or methods borrowed from Machine Learning, e.g. separating integral and differential data between a training set and a validation set. Such more advanced methods can be applied in future work.
Fig. 2 Distribution of the χ^{2} per degree offreedom considering sampled (random) H in H_{2}O thermal scattering data. The vertical black line indicates the position of χ^{2}∕N = 1. 
Fig. 3 Calculated over Experimental values for the considered criticality benchmarks. The case “Best” corresponds to a selected H in H_{2}O thermal scattering data case. 
Nominal (unperturbed) CAB parameters (file number 0) and sampled parameters leading to the minimum χ^{2} (file number 1871).
4 Conclusion
This work presents an estimation of the impact of the thermal scattering data for H in H_{2}O (or TSL) on k_{eff} of criticality benchmarks. The method used is the sampling of the CAB model parameters, leading to sampled thermal scattering data in the form of ACER files, and finally to sampled k_{eff} values. It was found that the impact of such TSL is modest compared to other nuclear data (less than 300 pcm). It was also found that the effect of TSL is highly non linear. Based on these results, a selection of “best performing” sampled TSL is proposed, using a simple χ^{2} definition. In conclusion, this paper is an application of the TMC and Petten methods towards nuclear data adjustment, more precisely thermal scattering data adjustment. If used in combination of traditional nuclear data evaluation methods, it can help to reach a better agreement between simulated and experimental integral values.
Author contribution statement
The authors have equally contributed to the various steps of technical work and to the preparation of the paper.
Acknowledgements
We would like to thank J.I. Marquez Damian from the European Spallation Source and G. Noguere from the CEA Cadarache for their support and efforts in helping for the production of random TSL data. This work would not have been possible without their help.
It was also partly funded by the European Union’s Horizon 2020 Research and Innovation Programme with the SANDA project, under grant agreement No 847552.
References
 A.J.M. Plompen, et al., The joint evaluated fission and fusion nuclear data library, JEFF3.3, Eur. Phys. J. A 56, 181 (2020) [Google Scholar]
 D.A. Brown et al., ENDF/BVIII.0: The 8^{th} Major Release of the Nuclear Reaction Data Library with CIELOproject Cross Sections, New Standards and Thermal Scattering Data, Nucl. Data Sheets 148, 1 (2018) [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]
 A. Koning, D. Rochman, J. Sublet, N. Dzysiuk, M. Fleming, S. van der Marck, TENDL: complete nuclear data library for innovative nuclear science and technology, Nucl. Data Sheets 155, 1 (2019) [Google Scholar]
 A. Aures, W. Bernnat, F. Bostelmann, J. Bousquet, B. KrzykaczHausmann, A. Pautz, K. Velkov, W. Zwermann, Reactor simulations with nuclear data uncertainties, proceedings of the ANS Best Estimate Plus Uncertainty International Conference (BEPU 2018), BEPU2018KNE7, Real Collegio, Lucca, Italy, May 1319, 2018 [Google Scholar]
 D. Rochman, A. Dokhane, A. Vasiliev, H. Ferroukhi, M. Hursin, Nuclear data uncertainties for core parameters based on Swiss BWR operated cycles, Ann. Nucl. Energy 148, 107727 (2020) [CrossRef] [Google Scholar]
 O. Cabellos, E. Castro, C. Ahnet, C. Holgado, Propagation of nuclear data uncertainties for PWR core analysis, Nucl. Eng. Technol. 46, 299 (2014) [CrossRef] [Google Scholar]
 I.C. Gauld, M.L. Williams, F. MichelSendis, J.S. Martinez, Integral nuclear data validation using experimental spent nuclear fuel compositions, Nucl. Eng. Technol. 49, 1226 (2017) [CrossRef] [Google Scholar]
 M. Frankl, M. Hursin, D. Rochman, A. Vasiliev, H. Ferroukhi, Nuclear data uncertainty quantification in criticality safety evaluations for spent nuclear fuel geological disposal, Appl. Sci. 11, 6499 (2021) [CrossRef] [Google Scholar]
 D. Rochman, A. Vasiliev, A. Dokhane, H. Ferroukhi, Uncertainties for Swiss LWR spent nuclear fuels due to nuclear data, Eur. Phys. J. Nuclear Sciences & Technologies 4, 6 (2018) [CrossRef] [EDP Sciences] [Google Scholar]
 G. Ilas, H. Liljenfeldt, Decay heat uncertainty for BWR used fuel due to modeling and nuclear data uncertainties, Nucl. Eng. Des. 319, 176 (2017) [CrossRef] [Google Scholar]
 A. Rintala, Evaluating the Effect of Decay and Fission Yield Data Uncertainty on BWR Spent Nuclear Fuel Source Term, in Proceedings of the 29^{th} international conference on nuclear energy for the new Europe, NENE 2020, Sept. 7–10, Portoroz, Slovenia (2020) [Google Scholar]
 I. Trivedi, J. Hou, G. Grasso, K. Ivanov, F. Franceschini, Nuclear data uncertainty quantification and propagation for safety analysis of leadcooled fast reactors, Sci. Technol. Nucl. Instal. 3961095 (2020) [Google Scholar]
 L.A. Bernstein, D.A. Brown, A.J. Koning, B.T. Rearden, C.E. Romano, A.A. Sonzogni, A.S. Voyles, W. Younes, Our future nuclear data needs, Annu. Rev. Nucl. Part. Sci. 69, 109 (2019) [CrossRef] [Google Scholar]
 F. Bostelmann, G. Ilas, W.A. Wieselquist, Key Nuclear Data Impacting Reactivity in Advanced Reactors, Oak Ridge National Laboratory, Report ORNL/TM2020/1557 (2020) [CrossRef] [Google Scholar]
 W.A. Wieselquist, K.S. Kim, G. Ilas, I.C. Gault, Comparison of burnup credit uncertainty quantification methods, in Proceedings of the ANS NCSD 2013  Criticality Safety in the Modern Era: Raising the Bar, Wilmington, NC, September 29October 3, 2013, on CDROM, American Nuclear Society, LaGrange Park, IL (2013) [Google Scholar]
 A. Vasiliev, J. Herrero, M. Pecchia, D. Rochman, H. Ferroukhi, S. Caruso, Preliminary assessment of criticality safety constraints for Swiss spent nuclear fuel loading in disposal canisters, Materials 12, 494 (2019) [CrossRef] [Google Scholar]
 I.C. Gauld, U. Mertyurek, Validation of BWR spent nuclear fuel isotopic predictions with applications to burnup credit, Nucl. Eng. Des. 345, 110 (2019) [CrossRef] [Google Scholar]
 Safety assessment and verification for nuclear power plants, IAEA Safety Standards Series NSG1.2 (International Atomic Energy Agency, Vienna, Austria, 2009) [Google Scholar]
 Best estimate safety analysis for nuclear power plants: Uncertainty evaluation, IAEA Safety Report Series No. 52 (International Atomic Energy Agency, Vienna, Austria, 2008) [Google Scholar]
 E.J. Bonano, J.E. Meacham, G.J. Appel, Radioactive waste management: it’s not all science and engineering, in Proceedings of the conference on Management of Spent Fuel from Nuclear Power Reactors: Learning from the Past, Enabling the Future, International Atomic Energy Agency, Vienna, Austria, 2020, on CDROM, IAEACN272/65, ID65 (2020) [Google Scholar]
 J.B. Briggs Ed., International Handbook of evaluated Criticality Safety, Benchmark Experiments, NEA/NSC/DOC(95)03/I (Organization for Economic Cooperation and Development, Nuclear Energy Agency, 2004) [Google Scholar]
 J.I. Marquez Damian, J.R. Granada, D.C. Malaspina, CAB models for water: a new evaluation of the thermal neutron scattering laws for light and heavy water in ENDF6 format, Ann. Nucl. Energy 65, 280 (2014) [CrossRef] [Google Scholar]
 G. Noguere, J.P. Scotta, C. De Saint Jean, P. Archier, Covariance matrices of the hydrogen neutron cross sections bound in light water for the JEFF3.1.1 neutron library, Ann. Nucl. Energy 104, 132 (2017) [CrossRef] [Google Scholar]
 A.J. Koning, D. Rochman, Towards sustainable nuclear energy: putting nuclear physics to work, Ann. Nucl. Energy 35, 2024 (2008) [CrossRef] [Google Scholar]
 L. Maul, Thermal Scattering Law Uncertainties and Propagation into Small Thermal Fission Reactors, PhD Thesis, University of New South Wales, Sydney, Australia (2018) [Google Scholar]
 J.P. Scotta, G. Noguère, J.I. Marquez Damian, Generation of the ^{1}H in H_{2}O neutron thermal scattering law covariance matrix of the CAB model, Eur. Phys. J. Nuclear Sciences & Technologies 4, 32 (2018) [CrossRef] [EDP Sciences] [Google Scholar]
 D. Rochman, A. Koning, Improving the H in H_{2}O thermal scattering data using the Petten method, Nucl. Sci. Eng. 172, 287 (2012) [CrossRef] [Google Scholar]
 D. Rochman, A. Koning, Evaluation and adjustement of the neutroninduced reactions of ^{63,65}Cu, Nucl. Sci. Eng. 170, 265 (2012) [CrossRef] [Google Scholar]
 D. Rochman, A. Koning, How to randomly evaluate nuclear data: a new method applied to ^{239}Pu, Nucl. Sci. Eng. 169, 68 (2011) [CrossRef] [Google Scholar]
 E. Bauge, M. Dupuis, S. Hilaire, S. Peru, A.J. Koning, D. Rochman, S. Goriely, Connecting the dots, or nuclear data in the age of supercomputing, Nucl. Data Sheets 118, 32 (2014) [CrossRef] [Google Scholar]
 D. Rochman, E. Bauge, A. Vasiliev, H. Ferroukhi, Correlation nusigmachi in the fast neutron range via integral information, Eur. J. Phys. Nuclear Sciences & Technologies 3, 14 (2017) [CrossRef] [EDP Sciences] [Google Scholar]
 D. Rochman, E. Bauge, A. Vasiliev, H. Ferroukhi, S. Pelloni, J.Ch. Koning, Sublet, Monte Carlo nuclear data adjustment via integral information, Eur. Phys. J. Plus 133, 537 (2018) [CrossRef] [Google Scholar]
 D. Rochman, A. Vasiliev, H. Ferroukhi, S. Pelloni, E. Bauge, A.J. Koning, Correlation nusigma for UPu in the thermal and resonance neutron range via integral information, Eur. Phys. J. Plus 134, 453 (2019) [Google Scholar]
 D. Kumar, S.B. Alam, H. Sjöstrand, J.M. Palau, C. De Saint Jean, Nuclear data adjustment using Bayesian inference, diagnostics for model fit and influence of model parameters, Eur. Phys. J. Web Conf. 239, 13003 (2020) [CrossRef] [EDP Sciences] [Google Scholar]
 C. De Saint Jean, P. Archier, E. Privas, G. Noguere, B. Habert, P. Tamagno, Evaluation of neutroninduced cross sections and their related covariances with physical constraints, Nucl. Data Sheets 148, 383 (2018) [CrossRef] [Google Scholar]
 C.W. Chapman, Thermal neutron scattering evaluation framework, PhD Thesis, Georgia Institute of Technology, USA, August 2017 [Google Scholar]
 D. Siefman, M. Hursin, D. Rochman, S. Pelloni, A. Pautz, Stochastic vs. sensitivitybased integral parameter and nuclear data adjustments, Eur. Phys. J. Plus 133, 429 (2018) [Google Scholar]
 E. Ivanov, C. de Saint Jean, V. Sobes, Nuclear data assimilation, scientific basis and current status, Eur. J. Phys. NuclearSciences & Technologies 7, 9 (2021) [Google Scholar]
 D. Rochman, M. Hursin, A. Vasiliev, H. Ferroukhi, Impact of H in H_{2}O thermalscattering data on depletion calculation: k_{eff}, nuclide inventory and decay heat, submitted to Eur. Phys. J. Nuclear Sciences & Technologies (2021) [Google Scholar]
 tslHinH2O, J.I. Marquez Damian, J.R. Granada, D. Roubtsov, https://github.com/marquezj/tslHinH2O, September 2020 [Google Scholar]
 J.P. Scotta, G. Noguere, D. Bernard, J.I. Marquez Damian, A. Santamarina, Impact of the thermal scattering law of H in H_{2}O on the isothermal temperature reactivity coefficients for UOX and MOX fuel lattices in cold operating conditions, Eur. Phys. J. Nuclear Sciences & Technologies 2, 28 (2016) [CrossRef] [EDP Sciences] [Google Scholar]
 R.E. MacFarlane, A.C. Kahler, Methods for processing ENDF/BVII with NJOY, Nucl. Data Sheets 111, 2739 (2010) [CrossRef] [Google Scholar]
Cite this article as: D. Rochman, A. Vasiliev, H. Ferroukhi, A. Koning, J.Ch. Sublet, Impact of H in H_{2}O thermal scattering data on criticality calculation: uncertainty and adjustment, EPJ Nuclear Sci. Technol. 8, 3 (2022)
All Tables
Uncertainties (standard deviation and MAD), skewness and kurtosis on criticality benchmarks due to thermal scattering data of H in H_{2}0 obtained from MCNP6 calculations. The base library (except the thermal scattering data) is JEFF4.0β. The statistical uncertainties are already removed.
Nominal (unperturbed) CAB parameters (file number 0) and sampled parameters leading to the minimum χ^{2} (file number 1871).
All Figures
Fig. 1 Probability density function for a number of criticality benchmarks, obtained by randomly changing the H in H_{2}O thermal scattering data. 

In the text 
Fig. 2 Distribution of the χ^{2} per degree offreedom considering sampled (random) H in H_{2}O thermal scattering data. The vertical black line indicates the position of χ^{2}∕N = 1. 

In the text 
Fig. 3 Calculated over Experimental values for the considered criticality benchmarks. The case “Best” corresponds to a selected H in H_{2}O thermal scattering data case. 

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.