https://doi.org/10.1051/epjn/2017009
Regular Article
Correlation ν̅_{p} − σ − χ in the fast neutron range via integral information
^{1}
Laboratory for Reactor Physics Systems Behaviour, Paul Scherrer Institut, Villigen, Switzerland
^{2}
CEA, DAM, DIF, 91297 Arpajon cedex, France
^{⁎} email: dimitrialexandre.rochman@psi.ch
Received:
18
January
2017
Received in final form:
20
March
2017
Accepted:
20
March
2017
Published online: 15 May 2017
This paper presents a Bayesian approach based on integral experiments to create correlations which do not appear with differential data. Some quantities such as the fission cross section (σ), neutron multiplicity (ν̅_{p}), neutron spectra (χ), etc. are usually neither modeled together nor measured in coincidence, thus there is no correlation matrices in evaluated nuclear data libraries. One can nevertheless use the information from integral experiments such as fast criticalitysafety benchmarks to correlate such quantities for possible inclusion in nuclear data libraries. A simple Bayesian set of equations is presented with random nuclear data, similarly to the usual methods applied with differential data. An example for ^{239}Pu is proposed.
© D. Rochman et al., published by EDP Sciences, 2017
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
Covariance information for basic nuclear data quantities is nowadays necessary for various types of nuclear applications. Many simulation tools are capable of using such matrices to propagate nuclear data uncertainties on final quantities, with either perturbation theories [1–6], or Monte Carlo sampling [3,7–13]. These results can for instance be used for the review procedure of new facilities, or during the safety assessment of new reactor core designs [14]. The origin of such covariance information can be found in various nuclear data libraries, such as ENDF/BVII [15,16], JENDL4.0 [17] or TENDL [18], and the production techniques can widely vary, see for instance references [19,20,22–26].
All the methods applied to experimental observables (measured cross sections, angular and energy distributions) provide nominal values, with uncertainties and correlations, for instance in terms of covariance matrices between incident neutron energies for a specific reaction channel, e.g. fission, or capture. With the help of nuclear reaction models, crosschannel covariance matrices can be obtained, for instance between the fission and the capture cross sections. Alternatively, the methods solely based on model parameters, such as in references [19,20], also lead to similar matrices.
One of the drawbacks of the existing approaches is that no correlations are proposed between quantities which are not simultaneously measured, or calculated within a unique theoretical model. For instance, there is currently no measurements providing in coincidence a fission rate (related to the fission cross section σ(n,f)) and the number of emitted prompt neutron per fission (socalled nubar, or ν̅_{p}). It can nevertheless be noted that some efforts are spent to obtain such data, as presented in reference [21]. In the following σ(n,f) and ν̅_{p} will simply be noted σ and ν, respectively. Similarly, the reaction models separately calculate cross sections and prompt neutron emitted spectra (χ) as independent quantities. Therefore, no theoretical correlation matrix can be obtained between σ, ν and χ. One may argue that if quantities are not measured or calculated together, there is no correlation between them. Nonetheless correlations and uncertainties are valid within the applied method and can therefore be calculated in a welldefined framework. With the example of σ − ν − χ, if there are no differential evidence of such correlation (at the level of the observables, i.e. cross sections and emitted neutrons), there are nevertheless indications from integral measurements. It is also important to realize that missing correlation matrices in evaluated files correspond to zero correlations. Such correlations can affect uncertainties for integral quantities, as uncorrelated nuclear data can be either independently sampled (for Monte Carlo uncertainty propagation), or can lead to simplified covariance matrices (for sensitivity & perturbation methods). Moreover, these zero correlation matrices are not visible and users are often not aware of the possible impact of these hidden correlations.
This paper proposes to evaluate correlation matrices between σ, ν and χ in the fast neutron range using integral information based on a simple Bayesian method, as presented in references [20,25]. As mentioned, such correlation matrices are not provided in nuclear data libraries, as they are in most of the cases based on differential measurements. But the use of integral benchmarks, such as fast criticalitysafety benchmarks [27] can provide such values, as presented in Section 2. Examples for the ^{239}Pu fission cross section (σ), its ν and χ will be presented in Section 3 for socalled fast benchmarks. Finally, an overview of the correlation matrices will be presented in Section 4, with possible other applications.
2 Correlation from integral benchmarks
As mentioned in the introduction, correlations between some quantities are not provided in libraries, unfortunately giving the presumption to the nuclear data users that they do not exist. One of the most important example concerns the fission cross section and ν for a given nucleus. Assuming that such quantities are estimated with an uncertainty and that design quantities are fixed, there is a small probability that their true values are equal to their mean estimates plus one standard deviation, for both of them at the same time. It is however more likely that from our simulation capabilities, their true values are anticorrelated: an estimated high fission cross section is probably to be compensated by a low ν (since the source term in the Boltzmann equation takes the form of σ × ν). A similar argument can be applied for the fission cross section and the hardness of χ (how high is the mean outgoing neutron energy, which drives the energy distribution of the fission neutrons).
This is the central point of this work: using numerical simulations based on these three quantities to extract their correlation matrices (e.g. with the Boltzmann transport equation). Derived integral quantities are numerous and easily calculated, and they can be compared to experimental integral data. In the present context, we will limit ourselves to criticality benchmarks with the calculation of the multiplication factor k_{eff}, for instance calculated with the neutron transport code MCNP [28]. In first approximations, k_{eff} is proportional to the product of ν and σ, thus these two quantities are anticorrelated for a fixed k_{eff} value.
The general proposed method is as follows. Within the assumptions that ν, σ and χ are not perfectly known and can be represented by independent probability distributions (prior assumption), it is possible to generate random realizations, and then to calculate random k_{eff} values for a specific system definition (geometry, content, etc.). Such random k_{eff} values can then be compared to a specific measurement of k_{eff} with its uncertainty (called k_{exp} ± Δk), using for instance a simplified chi2 definition. The correlation between ν and σ (and χ) can then be obtained using weights and the usual covariance definitions (leading to posterior probabilities for the 3 quantities of interest).
In the following, an example is given for ν and σ. Any of these two quantities can also be replaced by χ.
Independent probability distributions (prior pdf): the calculation of these pdf is based on the sampling of model parameters (as in the TMC and BFMC methods [19,20]). Model parameters are sampled a large number of times (with the index i = 1, … ,n) to generate full cross sections and other nuclear data quantities for ^{239}Pu from 0 to 20 MeV. The T6 system [18] is used to generate socalled random ENDF6 and ACE files, containing all necessary random nuclear data. Such random files for ^{239}Pu can be found in the TENDL2014 library [29]. As mentioned, in the T6 system, the ν values (as a function of the incident neutron energy) are independently calculated from σ, for all incident energies.
The systems are fast criticality benchmarks. The advantage of these benchmarks is that their k_{eff} is highly sensitive to the fast neutron range, thus avoiding the complications linked to the resolved and unresolved resonance range. MCNP6 is used to calculate k_{eff,i} for the specific random ACE file i for ^{239}Pu.
The comparison between n random calculated and experimental k_{eff,i = 1 …, n} is performed with the chi2 Q_{i} values and associated weights w_{i} (here, chi2 is called Q_{i} to differentiate it from the neutron spectra χ):(1)(2)Such formulation can easily be linked to the usual Bayesian likelihood.

Such weights are then assigned to the corresponding ν_{i} and σ_{i} which lead to k_{eff,i}. The final quantity for a specific benchmark consists of a matrix containing [ν_{i}, σ_{i}, w_{i}] for i = 1, … ,n. In the following, n = 10 000. Figure 1 presents the distribution of ν and σ with distinct colors being proportional to w_{i}: red for small w_{i} and black for w_{i} ≃ 1. The w_{i} are obtained using the pmf1 benchmark (or Jezebel). A simple spline function gives an idea of the correlation between ν and σ for their specific energies.
Fig. 1 Example of correlations between ν and σ for two different incident neutron energies for the fission cross section. The incident neutron energy for ν is 900 keV in both cases. The considered criticality benchmark is pmf1. 10 000 samples are presented; the distinct colors are proportional to w_{i}: red for small w_{i} and black for w_{i} ≃ 1.
Each of the random ENDF6 files in this work contains unique nuclear data realizations based on the sampling of the model parameters. They have already been used and characterized in previous work, see for instance references [2,3,26]. All possible nuclear data are varied, such as cross sections, angular distributions, emitted spectra and multiplicities. Therefore, the observed spread of k_{eff} based is the reflection of all varied quantities. Still, if one of them (e.g. the fission cross section at a specific energy) has a dominant impact on k_{eff}, there will be a correlation with k_{eff}, independently of the variations of the other variables. Such correlation can be observed as long as the nuclear data of interest has a significant impact, somewhat not overshadowed by the action of the other variations. Such correlation factors are already used in reactor physics to order important input factors and quantify their impacts. They are directly related to the importance factors R^{2}, see for instance references [5,30]. The R^{2} (or the correlation ρ) are very similar to the sensitivity vectors obtained in perturbation methods, but they also take into account the internal energyenergy correlation (for instance the fission cross section energyenergy correlation). Thus, correlations between input variables and the output might be found in a specific energy range, even if there is no sensitivity for this energy range, because of the energyenergy correlation with an energy region where the sensitivity is not negligible. This remark is important for the understanding of the following results, where some nonzero correlation values might be found for incident energies where the considered benchmark has a very small neutron population.
The correlation ρ(E_{ν}, E_{σ}) can be calculated for specific values of the incident neutron energies for ν(E_{ν}) and σ(E_{σ}). In the following, ρ(E_{ν}, E_{σ}), ν(E_{ν}) and σ(E_{σ}) will be noted ρ, ν and σ for simplicity. Considering the vector [ν_{i}, σ_{i}, w_{i}] with n = 10 000, ρ can be calculated as follows. Using the definition of weighted averages:and the definition of the weighted covariance factors:the correlation ρ between ν and σ is given by(3)Such correlation ρ can be obtained for different E_{σ} and E_{ν}, thus defining a full correlation matrix. In the example of Figure 1, the correlations are −0.10 and −0.34 for the two different energies of the fission cross section, 70 keV and 1.45 MeV, respectively (the energy of ν is the same and equal to 900 keV). In the case of χ, two energy grids are considered: the incident neutron energy and the outgoing neutron energy. The proposed set of equations can be applied to χ for a specific incident energy (for instance a incident neutron with an energy of 750 keV) and E represents the outgoing neutron energy.
With these basis, correlation matrices can now be calculated for ^{239}Pu, given a set of specific benchmarks. In the following, the statistical uncertainty of the k_{eff} calculation performed with MCNP6 is in the order of 20 pcm.
A similar approach was applied in reference [31] where correlations between the inelastic cross sections of the Cu isotopes and a shielding benchmarks were obtained. In this reference, the correlation was not based on a weighting approach, but rather on a simple accept/reject method.
As a final remark, the statistical significance of the calculated correlation can be assessed as in reference [32] using the Student's ttest:(4)with the degree of freedom equal to n − 2. The tvalue is evaluated by the usual tdistribution: with n = 10 000, a significance level of 95% is achieved for ρ > 1.6%. In the case of 99.9%, ρ > 3.1%.
3 Application to ^{239}Pu
The isotope ^{239}Pu presents several advantages from the nuclear data point of view. Its fission cross section, ν and χ are relatively well known in the fast neutron range (generally above a few tens of keV) and it is commonly used for criticality benchmarks. Additionally, the production of random ENDF6 and ACE files is relatively fast compared to other actinides, which allows to reach n = 10 000 without a large computer facility. Many criticalitysafety benchmarks are highly sensitive to ^{239}Pu cross sections, namely all the “plutonium metal fast” or pmf benchmarks. The thermal and intermediate benchmarks (plutonium solution thermal and the plutonium metal intermediate) are also very sensitive to ^{239}Pu, but mainly below a few tens of keV.
Using the set of equations presented in the previous section, the correlations between the three quantities ν, σ and χ can be obtained, given a set of benchmarks. This approach is very similar to the nuclear data evaluation with differential nuclear data. When evaluating for instance the ^{239}Pu fission cross section, measurements of the pointwise or groupwise cross section are considered. Due to the experimental correlation of the considered measurements (for instance over an energy range), the evaluated cross section will also be correlated with itself. Additionally, if the constrains from other differential cross section measurements are considered (such as the total or reaction cross sections), crosscorrelations between these different types will be obtained. The method proposed in this paper is therefore the counterpart for integral measurements.
As a consequence, the integral k_{eff} benchmarks are used in the following to obtain correlation matrices for nuclear data which have high sensitivities. The prior correlations used for ^{239}Pu are presented in Figure 2. Three distinct blocks with nonzero correlations can be seen and correspond to energyenergy correlations for the same quantity. As explained, there are no crosscorrelations between these three quantities. Specific patterns are only due to the considered models as no differential experiments are considered (e.g. the strong correlations observed for the χ matrix comes from the use of the MadlandNix model, see Ref. [18] for details).
The correlation matrix between ν and σ for different incident neutron energies are presented in Figure 3 for a specific benchmark.
As presented, such correlation values are depending on the type of considered benchmarks, namely fast or thermal.
In the case of the fast benchmark, two clear zones are visible with a border at the incident energies of the cross section ≈500–600 keV. Anticorrelations between ν and σ are expected from the benchmark itself in the energy region where the neutron population is important: a high neutron multiplicity can be compensated by a low fission cross section to obtain a specific k_{eff}. In the energy region where the neutron population is weak, no noticeable correlation is found, even if the correlation is gradually fading away, due to the mentioned internal correlation. The example of this benchmark is representative of the other fast cases.
In a similar way, the correlations between χ and ν can also be obtained. No explicit examples are provided here, as the final correlation matrix is presented in the next section. The impact of the additional anticorrelations can be seen in Table 1 where the prior and posterior k_{eff} and the calculated uncertainties are presented. The posterior values are calculated using the weights ω_{i} for both the averages and the uncertainties (standard deviations). The calculated over measured k_{eff} are also presented as C/E ratios to indicate the impact of the weights compared to the reference (measured) values. For all the considered benchmarks, the calculated uncertainties are strongly reduced (by more than a factor 2) and the average k_{eff} is in better agreement with the experimental value. To illustrate the change in the pdf, Figure 4 presents the prior and posterior distributions for two benchmarks: pmf1 and pmf13. In the case of the pmf1 benchmark, the central value of the prior distribution is already relatively close to the experimental reference. Therefore the posterior average is not drastically changed. On the contrary, the prior uncertainty is large compared to the experimental value, leading to a posterior distribution much narrower. In the case of the pmf13 benchmark, same observation can be done for the uncertainty, but the average value is also strongly changed, as the prior distribution is shifted compared to the experimental k_{eff}. In both cases, the posterior distributions reflect the constraints from the integral data.
Fig. 2 Prior correlation matrix for ^{239}Pu ν, σ and χ (for the incident neutron energy of 750 keV). The energy axis is for the incident neutrons for ν and σ, and for the outgoing neutron for χ. The X and Yaxis are in linear scale. 
Fig. 3 Correlation matrix between ^{239}Pu ν and σ considering the fast pmf1 benchmark. The X and Yaxis are in log scale. 
Prior and posterior average k_{eff} and uncertainties for selected benchmarks. Uncertainties Δk are given in pcm. C/E values are also indicated.
Fig. 4 Prior and posterior distributions of k_{eff} for two benchmarks: pmf1 and pmf13. The blue line indicate the experimental value with its uncertainty. 
4 Combination
Based on the previous results, the correlation matrix presented in Figure 2 can be updated for the three considered quantities. With a set of fast benchmarks, the weighted average correlation matrices can be calculated, as it can be done during the usual evaluation process with differential data (the weight is driven here by the experimental uncertainty of k_{eff}, but any other definition can be used). Considering 9 fast benchmarks (pmf1, pmf3, pmf9, pmf10, pmf19, pmf20, pmf24, pmf35 and pmf44, only the first configuration for each of these benchmarks), the posterior correlation matrix is calculated and presented in Figure 5.
The number of 9 benchmarks is certainly not enough to achieve a satisfactory convergence of the average correlation factors, but similarly to the evaluation of the cross sections, the number of experimental data is often limited. Nuclear data evaluators frequently consider such a small number of data; in the case of the maximum correlation, the 95% confidence interval is obtained within a ±0.16 correlation band. As mentioned in the previous section, nonzero crosscorrelation terms are observed. The strongest crosscorrelation links the fission cross section and the neutron multiplicity, which globally remain negative for all incident neutron energies. The crosscorrelation for the two other quantities are weaker, although not zero. In the context of a full evaluation for the inclusion in a specific library, a variety of sensitive benchmarks needs to be used. Also, the thermal and resonance energy range needs to be studied, with dedicated thermal and intermediate benchmarks. The present study nevertheless demonstrates the possibility of using integral values to calculate correlations which are not available otherwise. The observed negative correlation between the fission cross section and the neutron multiplicity will certainly impact the calculated uncertainties for various applications, for instance in reactor calculations. Alternatively, such method can be used to obtain correlation matrices between different isotopes. In the case of ^{238}U and ^{235}U, many of their nuclear data affect integral observables. An integral benchmarks with an intermediate ^{235}U enrichment, such as the “Big10 (or imf7 for intermediatemetalfast benchmark)” will certainly lead to nonnegligible crosscorrelations. An other method can be used following the work presented in reference [14]: the effect of nuclear data is observed on the simulation of realistic reactor core cycle operations. Such simulation can be compared to incore measurements and thus replacing the k_{eff} quantity used in this study.
It can be noticed that the update of the uncertainties is not proposed in this work, but only the (cross)correlations. The use of integral experiments for the reduction of uncertainties for differential quantities is a subject of discussion and no consensus is yet achieved. Compensating factors can play an important role and might not be applicable in all cases for a generalpurpose library. It is different for correlation matrices: they are already implicitly included in evaluations by not being present. Missing correlations are equal to zero correlation, and if not entered in a library, they can be considered as shadow correlations. The proposed method merely offers an alternative to turn them into explicit quantities (note that a consequence of modifying correlation factors can be a reevaluation of the uncertainties).
Fig. 5 Posterior correlation matrix for ^{239}Pu ν, σ and χ (for the incident neutron energy of 750 keV). The energy axis is for the incident neutrons for ν and σ, and for the outgoing neutron for χ. The X and Yaxis are in linear scale. 
5 Conclusion
Crosscorrelation matrices between the fission cross section, the neutron multiplicity and their spectra were calculated based on fast integral experiments. The integration of the criticality benchmarks follows a classical Bayesian approach. It is not intended to replace the traditional evaluation methods based on differential measurements, but the use of integral quantities can provide crosscorrelations for nuclear data apparently not connected by the reaction models or by differential experiments. We would like to remind that missing correlations are equivalent to zero correlations. If not explicitly given, these zero correlations are shadow correlations which affect calculated quantities. Therefore, the proposed method turns hidden correlations into actual information. As mentioned, it can be extended to quantify crossisotope correlations, as between ^{235}U and ^{238}U and finally, be part of a general evaluation process for a nuclear data library based on differential and integral constraints.
6 Perspectives
One of the goals of such work is to complement the covariance information provided in nuclear data libraries. The present example does not represent an approach which is general enough for the insertion in such library, but it is one step towards this direction. In this context, the following remarks can be mentioned:
the example of σ, ν and χ can be extended to other nuclear data quantities. Other cross sections, such as capture, will be positively correlated with fission and other ones using integral benchmarks. For the thermal and resonance range, correlations can also be extracted using a multigroup representation for the cross sections since the energy grid is extremely dense;
such correlation factors are complementary to the ones obtained from differential measurements and theoretical calculations of differential quantities. They are not contradictory and should be used together for a final evaluation, for instance using a set of differential and integral weights (as presented in Ref. [31]);
for a complete evaluation using the integral information, a variety of quantities needs to be involved, covering both a large energy region and divers integral quantities (not only k_{eff}, but also reaction rates, spectra indexes, activation, or emitted neutron spectra from pulsed spheres. Additionally, the correlations between benchmarks need to be taken into account, for instance using a generalized chi2 instead of equation (1). In this case, the experimental covariance matrices need to be assessed with care since their impact will be of importance;
as mentioned, integral benchmarks can also bring correlations between isotopes (e.g. Big10 k_{eff} for correlations between ^{235}U and ^{238}U, but also between various reactions using spectra indexes);
the current ENDF6 format used for the nuclear data evaluation does not allow for correlation matrices between quantities such as σ and ν. An update of this format is therefore necessary;
in the context of a full evaluation of covariance matrices, it is also important to provide “statistically converged” correlations. It is therefore necessary to use as many integral experiments as necessary to obtain correlation values which do not changed outside a given limit, when adding new experiments. Again, this problematic is very similar to the case of differential data.
The combination of the integral and differential correlations can be complemented by the adjustment of the variances, taking into account both sources of information. This would lead to fully updated covariance matrices, following a well defined and reproducible scheme. Such work would then be part of the elaboration of a nuclear data library based on models (for differential data), realistic model parameter distributions and integral constraints, as presented in reference [33].
References
 P. Romojaro, F. AlvarezVelarde, I. Kodeli, A. Stankovskiy, C.J. Diez, O. Cabellos, N. GarciaHerranz, J. Heyse, P. Schillebeeckx, G. Van den Eynde, G. Zerovnik, Nuclear data sensitivity and uncertainty analysis of effective neutron multiplication factor in various MYRRHA core configurations, Ann. Nucl. Energy 110, 330 (2017) [CrossRef] (In the text)
 M. Aufiero, M. Martin, M. Fratoni, XGPT: extending Monte Carlo generalized perturbation theory capabilities to continuousenergy sensitivity functions, Ann. Nucl. Energy 96, 295 (2016) [CrossRef] (In the text)
 D. Rochman, A. Vasiliev, H. Ferroukhi, T. Zhu, S.C. van der Marck, A.J. Koning, Nuclear data uncertainty for criticalitysafety: Monte Carlo vs. linear perturbation, Ann. Nucl. Energy 92, 150 (2016) [CrossRef] (In the text)
 M. Pusa, Incorporating sensitivity and uncertainty analysis to a lattice physics code with application to CASMO4, Ann. Nucl. Energy 40, 153 (2012) [CrossRef]
 W. Zwermann, A. Aures, L. Gallner, V. Hannstein, B. KrzykaczHausmann, K. Velkov, J.S. Martinez, Nuclear data uncertainty and sensitivity analysis with XSUSA for fuel assembly depletion calculations, Nucl. Eng. Technol. 46, 343 (2014) [CrossRef] (In the text)
 M.L. Williams, B.T. Rearden, SCALE6 Sensitivity/uncertainty methods and covariance data, Nucl. Data Sheets 109, 2796 (2008) [CrossRef] (In the text)
 A. Hoefer, O. Buss, M. Hennebach, M. Schmid, P. Porsch, MOCABA: a general Monte CarloBayes procedure for improved predictions of integral functions of nuclear data, Ann. Nucl. Energy 77, 514 (2015) [CrossRef] (In the text)
 C.J. Diez, O. Cabellos, D. Rochman, A.J. Koning, J.S. Martinez, Monte Carlo uncertainty propagation approaches in ADS burnup calculations, Ann. Nucl. Energy 54, 27 (2013) [CrossRef]
 N. GarciaHerranz, O. Cabellos, J. Sanz, J. Juan, J.C. Kuijper, Propagation of statistical and nuclear data uncertainties in Monte Carlo burnup calculations, Ann. Nucl. Energy 35, 714 (2008) [CrossRef]
 E. Alhassan, H. Sjostrand, P. Helgesson, M. Osterlund, S. Pomp, A.J. Koning, D. Rochman, On the use of integral experiments for uncertainty reduction of reactor macroscopic parameters within the TMC methodology, Prog. Nucl. Energy 88, 43 (2016) [CrossRef]
 D. Rochman, O. Leray, A. Vasiliev, H. Ferroukhi, A.J. Koning, M. Fleming, J.C. Sublet, A Bayesian Monte Carlo method for fission yield covariance information, Ann. Nucl. Energy 95, 125 (2016) [CrossRef]
 P. Helgesson, D. Rochman, H. Sjostrand, E. Alhassan, A.J. Koning, UO_{2} versus MOX: propagated nuclear data uncertainty for k_{eff}, with burnup, Nucl. Sci. Eng. 177, 321 (2014) [CrossRef]
 D. Rochman, O. Leray, M. Hursin, H. Ferroukhi, A. Vasiliev, A. Aures, F. Bostelmann, W. Zwermann, O. Cabellos, C.J. Diez, J. Dyrda, N. GarciaHerranz, E. Castro, S. van der Marck, H. Sjostrand, A. Hernandez, M. Fleming, J.Ch. Sublet, L. Fiorito, Nuclear data uncertainties for typical LWR fuel assemblies and a simple reactor core, Nucl. Data Sheets 139, 1 (2017) [CrossRef] (In the text)
 O. Leray, H. Ferroukhi, M. Hursin, A. Vasiliev, D. Rochman, Methodology for core analysis with nuclear data uncertainty quantification and application to Swiss PWR operated cycles, Ann. Nucl. Energy (2016) (In the text)
 M.B. Chadwick, P. Oblozinsky, M. Herman, N.M. Greene, R.D. McKnight, D.L. Smith, P.G. Young, R.E. MacFarlane, G.M. Hale, S.C. Frankle, A.C. Kahler, T. Kawano, R.C. Little, D.G. Madland, P. Moller, R.D. Mosteller, P.R. Page, P. Talou, H. Trellue, M.C. White, W.B. Wilson, R. Arcilla, C.L. Dunford, S.F. Mughabghab, B. Pritychenko, D. Rochman, A.A. Sonzogni, C.R. Lubitz, T.H. Trumbull, J.P. Weinman, D.A. Brown, D.E. Cullen, D.P. Heinrichs, D.P. McNabb, H. Derrien, M.E. Dunn, N.M. Larson, L.C. Leal, A.D. Carlson, R.C. Block, J.B. Briggs, E.T. Cheng, H.C. Huria, M.L. Zerkle, K.S. Kozier, A. Courcelle, V. Pronyaev, S.C. van der Marck, ENDF/BVII.0: next generation evaluated nuclear data library for nuclear science and technology, Nucl. Data Sheets 107, 2931 (2006) [NASA ADS] [CrossRef] (In the text)
 M. Chadwick, M. Herman, P. Oblozinsky, M.E. Dunn, Y. Danon, A.C. Kahler, D.L. Smith, B. Pritychenko, G. Arbanas, R. Arcilla, R. Brewer, D.A. Brown, R. Capote, A.D. Carlson, Y.S. Cho, H. Derrien, K. Guber, G.M. Hale, S. Hoblit, S. Holloway, T.D. Johnson, T. Kawano, B.C. Kiedrowski, H. Kim, S. Kunieda, N.M. Larson, L. Leal, J.P. Lestone, R.C. Little, E.A. McCutchan, R.E. MacFarlane, M. MacInnes, C.M. Mattoon, R.D. McKnight, S.F. Mughabghab, G.P.A. Nobre, G. Palmiotti, A. Palumbo, M.T. Pigni, V.G. Pronyaev, R.O. Sayer, A.A. Sonzogni, N.C. Summers, P. Talou, I.J. Thompson, A. Trkov, R.L. Vogt, S.C. van der Marck, A. Wallner, M.C. White, D. Wiarda, P.G. Young, 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] (In the text)
 K. Shibata, O. Iwamoto, T. Nagagawa, N. Iwamoto, A. Ichihara, S. Kunieda, S. Chiba, K. Furutaka, N. Otuka, T. Ohasawa, T. Murata, H. Matsunobu, A. Zukeran, S. Kamada, J. Katakura, JENDL4.0: a new library for nuclear science and engineering, J. Nucl. Sci. Technol. 48, 1 (2011) [CrossRef] (In the text)
 A.J. Koning, D. Rochman, Modern nuclear data evaluation with the TALYS code system, Nucl. Data Sheets 113, 2841 (2012) [CrossRef] (In the text)
 A.J. Koning, D. Rochman, Towards sustainable nuclear energy: putting nuclear physics to work, Ann. Nucl. Energy 35, 2024 (2008) [CrossRef] (In the text)
 E. Bauge, P. DossantosUzarralde, Evaluation of the covariance matrix of ^{239}Pu neutronic cross sections in the continuum using the backwardforward MonteCarlo method, J. Korean Phys. Soc. 59, 1218 (2011) [CrossRef] (In the text)
 H.Y. Lee, R.C. Haight, T.A. Bredeweg, M. Devlin, N. Fotiades, M. Jandel, A. Laptev, R.O. Nelson, J.M. O'Donnell, B.A. Perdue, T.N. Taddeucci, J.L. Ullmann, S.A. Wender, M.C. White, C.Y. Wu, A. Chyzr, R.A. Henderson, E. Kwan, Prompt fission neutron spectrum study at LANSCE: CHINU project, in Proceedings of the Fifth International Conference on ICFN5, Sanibel Island, Florida, USA, 4–10 November 2012 (2012) (In the text)
 H. Leeb, St. Gundacker, D. Neudecker, V. Wildpaner, The GENEUS Project – development of an evaluation tool, J. Korean Phys. Soc. 59, 1230 (2011) [CrossRef] (In the text)
 R. Capote, D.L. Smith, An investigation of the performance of the unified Monte Carlo method of neutron cross section data evaluation, Nucl. Data Sheets 109, 2725 (2008) [CrossRef]
 R. Capote, D. Smith, A. Trkov, M. Meghzifene, A new formulation of the Unified Monte Carlo approach (UMCB), and crosssection evaluation for the dosimetry reaction ^{55}Mn(n,γ)^{56}Mn, J. ASTM Int. 9, JAI104115 (2012) [CrossRef]
 A.J. Koning, Bayesian Monte Carlo method for nuclear data evaluation, Eur. Phys. J. A 51, 184 (2015) [CrossRef] [EDP Sciences] (In the text)
 D. Rochman, A.J. Koning, S.C. van der Marck, Improving neutronics simulations and uncertainties via a selection of nuclear data, Eur. Phys. J. A 51, 182 (2015) [CrossRef] [EDP Sciences] (In the text)
 J.B. Briggs Ed., International Handbook of evaluated Criticality Safety Benchmark Experiments, NEA/NSC/DOC(95)03/I (Organisation for Economic Cooperation and Development, Nuclear Energy Agency, 2004) (In the text)
 T. Goorley, MCNP 6.1.1–Beta Release Notes, Los Alamos National Laboratory, Report LAUR1424680, June 2014 (In the text)
 A.J. Koning, D. Rochman, J. Kopecky, J.Ch. Sublet, E. Bauge, S. Hilaire, P. Romain, B. Morillon, H. Duarte, S. van der Marck, S. Pomp, H. Sjostrand, R. Forrest, H. Henriksson, O. Cabellos, S. Goriely, J. Leppanen, H. Leeb, A. Plompen, R. Mills, TENDL2014: TALYSbased evaluated nuclear data library (2014), ftp://ftp.nrg.eu/pub/www/talys/tendl2014/tendl2014.html (In the text)
 T. Zhu, A. Vasiliev, H. Ferroukhi, D. Rochman, A. Pautz, Testing the samplingbased NUSSRF tool for the nuclear data related global sensitivity analysis with Monte Carlo neutronics calculations, Nucl. Sci. Eng. 69, 184 (2016) (In the text)
 D. Rochman, A.J. Koning, Evaluation and adjustment of the neutroninduced reactions of ^{63,65}Cu, Nucl. Sci. Eng. 170, 265 (2012) [CrossRef] (In the text)
 A. Vasiliev, D. Rochman, M. Pecchia, H. Ferroukhi, Energies 9, 1039 (2016) [CrossRef] (In the text)
 E. Bauge, M. Dupuis, S. Hilaire, S. Péru, A.J. Koning, D. Rochman, S. Goriely, Nucl. Data Sheets 118, 32 (2014) [CrossRef] (In the text)
Cite this article as: Dimitri Rochman, Eric Bauge, Alexander Vasiliev, Hakim Ferroukhi, Correlation ν̅_{p} − σ − χ in the fast neutron range via integral information, EPJ Nuclear Sci. Technol. 3, 14 (2017)
All Tables
Prior and posterior average k_{eff} and uncertainties for selected benchmarks. Uncertainties Δk are given in pcm. C/E values are also indicated.
All Figures
Fig. 1 Example of correlations between ν and σ for two different incident neutron energies for the fission cross section. The incident neutron energy for ν is 900 keV in both cases. The considered criticality benchmark is pmf1. 10 000 samples are presented; the distinct colors are proportional to w_{i}: red for small w_{i} and black for w_{i} ≃ 1. 

In the text 
Fig. 2 Prior correlation matrix for ^{239}Pu ν, σ and χ (for the incident neutron energy of 750 keV). The energy axis is for the incident neutrons for ν and σ, and for the outgoing neutron for χ. The X and Yaxis are in linear scale. 

In the text 
Fig. 3 Correlation matrix between ^{239}Pu ν and σ considering the fast pmf1 benchmark. The X and Yaxis are in log scale. 

In the text 
Fig. 4 Prior and posterior distributions of k_{eff} for two benchmarks: pmf1 and pmf13. The blue line indicate the experimental value with its uncertainty. 

In the text 
Fig. 5 Posterior correlation matrix for ^{239}Pu ν, σ and χ (for the incident neutron energy of 750 keV). The energy axis is for the incident neutrons for ν and σ, and for the outgoing neutron for χ. The X and Yaxis are in linear scale. 

In the text 