Recognition: 2 theorem links
· Lean TheoremGeneration of fission yield covariance matrices and its application in uncertainty analysis of decay heat
Pith reviewed 2026-05-10 20:15 UTC · model grok-4.3
The pith
Fission yield covariance matrices from generalized least squares reduce decay heat uncertainties and highlight decay energy contributions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Applying the generalized least squares updating approach to fission yield data, constrained by conservation equations and chain yields, produces covariance matrices that, when used in decay heat summation calculations, strongly reduce the propagated uncertainty from fission yields compared to uncorrelated data, with decay energy data becoming the major contributor; the uncorrelated case gives a constant approximately 4 percent uncertainty that dominates beyond 100 seconds cooling time.
What carries the argument
The generalized least squares updating procedure that enforces basic physical conservation equations and chain yield data to derive complete covariance matrices for fission yields from evaluated nuclear data files.
If this is right
- Decay heat uncertainty from yields drops below that from decay data across the cooling time range.
- Relative uncertainties reach approximately 10 percent for ENDF/B-VIII.0 and JEFF-3.3 and 5 percent for JENDL-5 at 0.1 seconds cooling time.
- Uncertainties decrease to about 1 percent at 10^5 seconds for all three libraries.
- The correlations alter the contributions and sensitivity coefficients of important fission products to the total decay heat.
Where Pith is reading between the lines
- Similar covariance generation could be applied to fission yields for other isotopes to improve predictions in different reactor fuel cycles.
- Direct experimental measurements of decay heat variability would provide a test for the accuracy of these generated covariances.
- Further reduction in decay heat uncertainty would require improved precision in decay energy measurements for key fission products.
Load-bearing premise
The generalized least squares procedure constrained solely by conservation equations and chain yield data produces covariance matrices that accurately capture all relevant physical correlations in fission yield uncertainties.
What would settle it
Comparison of the calculated decay heat uncertainties against direct experimental measurements of decay heat for thermal fission of U-235 at multiple cooling times, checking whether the reduced uncertainty with covariances matches the observed spread better than the uncorrelated case.
Figures
read the original abstract
The uncertainties and covariance matrices of fission yield are important in the uncertainty analysis of decay heat. At present, there are no covariance matrixes of fission yield given in the evaluated nuclear data library, although they have provided the uncertainties with good estimates. In this work, the generalized least squares (GLS) updating approach was adopted to evaluate the fission yield covariances with the constraints from basic physical conservation equation and chain yield data, using the nuclear data files from ENDF/B-VIII.0, JENDL-5 and JEFF-3.3. Based on these original and updated data, summation calculation was performed for fission pulse decay heat of thermal neutron-induced fission of $^{235}$U. The uncertainties of decay heat were obtained through generalized perturbation theory, including the uncertainties propagated from fission yield, decay energy, decay constant and branching ratio. The original uncorrelated yield data contributes a $\sim 4 \%$ uncertainty at all times and dominates the decay heat uncertainty at cooling times longer than \SI{100}{s}. With the generated covariance matrixes, the uncertainty of calculated decay heat is strongly reduced and decay energy data makes a major contribution in general. The relative uncertainties at cooling time \SI{0.1}{\second} are $\sim$10$\%$ for ENDF/V-VIII.0 and JEFF-3.3 and $\sim$5$\%$ for JENDL-5 and those at cooling time 10$^{5}$ s are about 1$\%$ for three libraries. The influence of the GLS updating procedure on the contributions of important fission products to decay heat and their sensitive coefficients was also discussed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper describes the generation of fission yield covariance matrices for thermal neutron-induced fission of 235U using a generalized least squares (GLS) updating procedure applied to data from ENDF/B-VIII.0, JENDL-5, and JEFF-3.3. The GLS step incorporates constraints from basic physical conservation equations and chain yield values. These updated covariances (along with the original uncorrelated yields) are then used in summation calculations of fission pulse decay heat, with uncertainties propagated via generalized perturbation theory including contributions from fission yields, decay energies, decay constants, and branching ratios. The central results are that the original uncorrelated yields produce a roughly constant ~4% uncertainty that dominates decay heat uncertainty for cooling times >100 s, while the GLS-updated covariances strongly reduce the yield-driven uncertainty component, shifting dominance to decay energy data and yielding relative decay heat uncertainties of ~5-10% at 0.1 s and ~1% at 10^5 s.
Significance. If the generated covariances are physically complete, the work demonstrates a practical route to incorporating yield correlations into decay heat uncertainty quantification, which is relevant for nuclear safety analyses and spent-fuel heat-load assessments. The multi-library comparison and explicit separation of uncertainty sources provide a clear illustration of how proper covariance treatment can alter the dominant contributors. The approach is reproducible in principle since it relies on publicly available libraries and standard GLS updating.
major comments (2)
- [§2] §2 (GLS updating procedure): The method enforces only conservation laws and chain-yield constraints; no additional physical correlations (e.g., from pairing, shell effects, or fission dynamics) are imposed. This is load-bearing for the central claim because the reported strong reduction in decay-heat uncertainty (from the original ~4% level) relies on the completeness of the resulting negative correlations; without external validation against model-generated or measured covariances, it is unclear whether the reduction reflects physics or the limited constraint set.
- [§4] §4 (decay heat uncertainty results): The dominance shift to decay energy and the numerical uncertainty values (e.g., ~5-10% at 0.1 s) are presented without sensitivity tests to the choice of GLS constraints or to the inclusion of additional correlation sources; this directly affects the robustness of the claim that yield covariances no longer dominate after 100 s.
minor comments (2)
- [Abstract and §1] The abstract and §1 would benefit from a brief statement of the specific fission products or mass chains that drive the largest sensitivity coefficients after GLS updating.
- [§2] Notation for the GLS-updated covariance matrix (e.g., how the posterior covariance is denoted relative to the prior) should be defined explicitly in the methods section to aid reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. The points raised highlight important aspects of the GLS procedure and the robustness of the uncertainty results. We address each major comment below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [§2] §2 (GLS updating procedure): The method enforces only conservation laws and chain-yield constraints; no additional physical correlations (e.g., from pairing, shell effects, or fission dynamics) are imposed. This is load-bearing for the central claim because the reported strong reduction in decay-heat uncertainty (from the original ~4% level) relies on the completeness of the resulting negative correlations; without external validation against model-generated or measured covariances, it is unclear whether the reduction reflects physics or the limited constraint set.
Authors: We agree that the GLS updating relies on the specified constraints and does not incorporate additional nuclear-structure correlations. The procedure is designed to generate covariances consistent with the independent-yield uncertainties and the conservation/chain-yield data already present in the evaluated libraries. The observed reduction in decay-heat uncertainty therefore reflects the effect of these minimal constraints rather than a claim of complete physical realism. In the revised manuscript we will add an explicit discussion of this limitation, stating that the generated covariances represent those implied by the chosen constraints and that further correlations (e.g., from pairing or fission dynamics) could be introduced in future extensions. We will also note that the current results provide a baseline for the impact of yield correlations derived from basic physical laws. revision: partial
-
Referee: [§4] §4 (decay heat uncertainty results): The dominance shift to decay energy and the numerical uncertainty values (e.g., ~5-10% at 0.1 s) are presented without sensitivity tests to the choice of GLS constraints or to the inclusion of additional correlation sources; this directly affects the robustness of the claim that yield covariances no longer dominate after 100 s.
Authors: We acknowledge that the manuscript does not include explicit sensitivity tests on the GLS constraint set or on hypothetical additional correlations. To strengthen the robustness discussion, the revised version will incorporate a short sensitivity study: (i) repeating the GLS update with relaxed chain-yield tolerances and (ii) a qualitative assessment of how stronger negative correlations would further suppress the yield contribution. These additions will confirm that the shift in dominance to decay-energy uncertainties remains under moderate variations of the constraint strength, while clearly stating the assumptions involved. revision: yes
Circularity Check
No significant circularity; GLS covariance generation and GFT uncertainty propagation are independent of inputs by construction.
full rationale
The paper applies the standard generalized least squares procedure to produce fission-yield covariance matrices from external evaluated libraries (ENDF/B-VIII.0, JENDL-5, JEFF-3.3) subject only to conservation laws and chain-yield constraints. Decay-heat uncertainties are then obtained by generalized perturbation theory summation calculations that propagate the resulting matrices together with independent decay-data uncertainties. Neither the covariance construction nor the reported uncertainty reduction is equivalent to the original uncorrelated yields or to any self-citation; the numerical outcome is a direct, falsifiable consequence of the external data and the chosen constraints. No load-bearing step reduces to a definition, a fitted parameter renamed as prediction, or an author-specific uniqueness theorem.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Fission product yields must obey mass, charge, and other conservation laws.
- domain assumption Chain yield measurements provide reliable independent constraints on the fission yield distributions.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
generalized least squares (GLS) updating approach was adopted to evaluate the fission yield covariances with the constraints from basic physical conservation equation and chain yield data
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The original uncorrelated yield data contributes a ∼4% uncertainty at all times
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Z. Zhang, Q. Zou, N. Gui et al., Prediction and anal- ysis of decay heat transfer in the core of the peb- ble bed reactor, Prog. Nucl. Energy 173 (2024) 105253. doi:10.1016/j.pnucene.2024.105253
-
[2]
Y . Chang, M. Wang, J. Zhang et al., Best estimate plus uncer- tainty analysis of the China advanced large-scale PWR during LBLOCA scenarios, Int. J. Adv. Nucl. Reactor Des. Technol. 2 (2020) 34–42. doi:10.1016/j.jandt.2020.07.002
-
[3]
J. Li, R. Li, D. She et al., V alidation of NUIT for De- cay Heat Predictions of PWR Spent Fuel Assemblies, V ol. V olume 2: Nuclear Fuel and Material, Reactor Physics and Transport Theory, and Fuel Cycle Technology of International Conference on Nuclear Engineering, 2022, p. V002T02A008. doi:10.1115/ICONE29-89241
-
[4]
J.-i. Katakura, Uncertainty Analyses of Decay Heat Sum- mation Calculations Using JENDL, JEFF, and ENDF Files, J. Nucl. Sci. Technol. 50 (8) (2013) 799–807. doi:10.1080/00223131.2013.808004
-
[5]
G. Chiba, Consistent Adjustment of Radioactive Decay and Fission Yields Data with Measurement Data of Decay Heat and β -Delayed Neutron Activities, Ann. Nucl. Energy 101 (2017) 23–30. doi:10.1016/j.anucene.2016.09.054
-
[7]
F. Wang, X. Huang, Study on Decay Heat Calculation Method for Fast Neutron Fission of 235U, Atomic En- ergy Science and Technology 55 (zengkan) (2021) 10-16 doi:10.7538/yzk.2020.youxian.0859
-
[8]
J. Ma, H. Guo, H. Huang, Decay Heat Uncertainty Quan- tificationBased on Stochastic Sampling Method, Atomic En- ergy Science and Technology 55 (6) (2024) 1280-1286 doi:10.7538/yzk.2023.youxian.0766
-
[9]
D. Brown, M. Chadwick, R. Capote et al., ENDF/B-VIII.0: The 8 th Major Release of the Nuclear Reaction Data Library with CIELO-project Cross Sections, New Standards and Ther- mal Scattering Data, Nucl. Data Sheets 148 (2018) 1–142. doi:10.1016/j.nds.2018.02.001
-
[10]
O. Iwamoto, N. Iwamoto, S. Kunieda et al., Japanese evaluated 11 nuclear data library version 5: JENDL-5, J. Nucl. Sci. Technol. 60 (1) (2023) 1–60. doi:10.1080/00223131.2022.2141903
-
[11]
A. J. M. Plompen, O. Cabellos, C. De Saint Jean et al., The joint evaluated fission and fusion nuclear data library, JEFF- 3.3, Eur. Phys. J. A 56 (2020) 181. doi:10.1140/epja/s10050- 020-00141-9
-
[12]
Katakura, JENDL FP decay data file 2011 and fission yields data file 2011, Tech
J.-i. Katakura, JENDL FP decay data file 2011 and fission yields data file 2011, Tech. Rep. JAEA-Data/Code 2011-025, Japan Atomic Energy Agency (2012)
work page 2011
-
[13]
L. Fiorito, D. Piedra, O. Cabellos et al., Inventory calculation and nuclear data uncertainty propagation on light water reactor fuel using ALEPH-2 and SCALE 6.2, Ann. Nucl. Energy 83 (2015) 137–146. doi:10.1016/j.anucene.2015.03.046
-
[14]
Y . Wang, M. Cui, J. Guo et al., Lognormal-based sam- pling for fission product yields uncertainty propagation in pebble-bed htgr, Sci. Technol. Nucl. Install. 2020 (2020) 21. doi:10.1155/2020/8014521
-
[15]
Z. Lu, T. Zu, L. Cao et al., Generation and anal- ysis of independent fission yield covariances based on GEF model code, EPJ Web Conf. 281 (2023) 00015. doi:10.1051/epjconf/202328100015
-
[16]
L. Fiorito, A. Stankovskiy, G. V an Den Eynde et al., Gener- ation of fission yield covariances to correct discrepancies in the nuclear data libraries, Ann. Nucl. Energy 88 (2016) 12–23. doi:10.1016/j.anucene.2015.10.027
-
[17]
K. Tsubakihara, S. Okumura, C. Ishizuka et al., Evalu- ation of fission product yields and associated covariance matrices, J. Nucl. Sci. Technol. 58 (2) (2021) 151–165. doi:10.1080/00223131.2020.1813643
-
[18]
E. F. Matthews, L. A. Bernstein, W. Y ounes, Stochastically estimated covariance matrices for independent and cumula- tive fission yields in the ENDF/B-VIII.0 and JEFF-3.3 eval- uations, At. Data Nucl. Data Tables 140 (2021) 101441. doi:10.1016/j.adt.2021.101441
-
[19]
A. E. Lovell, T. Kawano, P . Talou, Calculated covariance ma- trices for fission product yields using BeoH, EPJ Web Conf. 281 (2023) 00018. doi:10.1051/epjconf/202328100018
-
[20]
N. Shu, S. Zhu, T. Zu et al., Fission Yield Evaluation Method Based upon Zp Model and Uncertainty Analysis in Burnup Calculation, Atomic Energy Science and Technology 56 (2022) 927–936. doi:10.7538/yzk.2022.youxian.0188
-
[21]
Z.-A. Wang, J. Pei, Y . Liu et al., Bayesian Evaluation of In- complete Fission Yields, Phys. Rev. Lett. 123 (2019) 122501. doi:10.1103/PhysRevLett.123.122501
-
[22]
Z.-A. Wang, J. Pei, Optimizing multilayer Bayesian neural net- works for evaluation of fission yields, Phys. Rev. C 104 (2021) 064608. doi:10.1103/PhysRevC.104.064608
-
[23]
C. Y . Qiao, J. C. Pei, Z. A. Wang, et al., Bayesian evaluation of charge yields of fission fragments of 239U, Phys. Rev. C 103 (2021) 034621. doi:10.1103/PhysRevC.103.034621
-
[24]
M.-X. Xiao, X.-J. Bao, Z. Wei et al., Bayesian evaluation of energy dependent neutron induced fission yields, Chin. Phys. C 47 (12) (2023) 124102. doi:10.1088/1674-1137/acf7b5
-
[25]
Q.-F. Song, L. Zhu, H. Guo et al., V erification of neutron- induced fission product yields evaluated by a tensor decomp- sition model in transport-burnup simulations, Nucl. Sci. Tech. 34 (2) (2023) 32. doi:10.1007/s41365-023-01176-5
-
[26]
D. Y . Huo, Z. Wei, K. Wu et al., Evaluation of pre-neutron- emission mass distributions in induced fission of typical ac- tinides based on Monte Carlo dropout neural network, Eur. Phys. J. A 59 (11) (2023) 265. doi:10.1140/epja/s10050-023- 01189-z
- [27]
-
[28]
D. L. Smith, Probability, Statistics, and Data Uncertainties in Nuclear Science and Technology, American Nuclear Society, 1991
work page 1991
-
[29]
Lewins, The Time-Dependent Importance of Neutrons and Precursors, Nucl
J. Lewins, The Time-Dependent Importance of Neutrons and Precursors, Nucl. Sci. Eng. 7 (3) (1960) 268–274. doi:10.13182/NSE60-A25713
-
[30]
G. Chiba, T. Narabayashi, Uncertainty quantification of total delayed neutron yields and time-dependent de- layed neutron emission rates in frame of summation calculations, Ann. Nucl. Energy 85 (2015) 846–855. doi:10.1016/j.anucene.2015.07.001
-
[31]
N. Schunck, D. Regnier, Theory of nuclear fis- sion, Prog. Part. Nucl. Phys. 125 (2022) 103963. doi:10.1016/j.ppnp.2022.103963
-
[32]
R. W. Mills, Fission product yield evaluation, Ph.D. thesis, University of Birmingham (1995)
work page 1995
-
[33]
T. England, B. Rider, Evaluation and compilation of fission product yields 1993, Tech. Rep. LA-UR-94-3106, Los Alamos National Lab. (1995)
work page 1993
-
[34]
A. Nichols, D. Aldama, M. V erpelli, Handbook of Nuclear Data for Safeguards: Database Extensions, August 2008., Tech. Rep. INDC(NDS)-0534, International Atomic Energy Agency (2008)
work page 2008
-
[35]
R. W. Mills, A new UK fission yield evaluation UKFY3.7, EPJ Web of Conferences 146 (2017) 04008. doi:10.1051/epjconf/201714604008
-
[36]
H. Bateman, Solution of a system of differential equations oc- curring in the theory of radioactive transformations, in: Proc. Cambridge Philos. Soc., V ol. 15, 1910, pp. 423–427
work page 1910
-
[37]
A. Ladshaw, A. I. Wiechert, Y .-h. Kim et al., Algorithms and algebraic solutions of decay chain differential equations for sta- ble and unstable nuclide fractionation, Comput. Phys. Com- mun. 246 (2020) 106907. doi:10.1016/j.cpc.2019.106907
-
[38]
IAEA, CoNDERC: Compilation of Nuclear Data Experi- ments for Radiation Characterisation, https://www-nds. iaea.org/conderc/ (2019)
work page 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.