Next Article in Journal / Special Issue
Illustrating Important Effects of Second-Order Sensitivities on Response Uncertainties in Reactor Physics
Previous Article in Journal / Special Issue
Choosing Transport Events for Initiating Splitting and Rouletting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

1
Naval Nulcear Laboratory, P.O. Box 1072, Schenectady, NY 12301-1072, USA
2
Naval Nulcear Laboratory, P.O. Box 79, West Mifflin, PA 15122-0079, USA
*
Author to whom correspondence should be addressed.
J. Nucl. Eng. 2021, 2(2), 105-113; https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020011
Submission received: 1 October 2020 / Revised: 10 February 2021 / Accepted: 18 February 2021 / Published: 7 April 2021
(This article belongs to the Special Issue Selected Papers from PHYSOR 2020)

Abstract

:
Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ -ZrHx and ϵ -ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

1. Introduction

Zirconium hydride is used as a moderator material in TRIGA reactors and has been historically used as a moderator in other reactors requiring high proton density for moderation (e.g., SNAP-10A [1]). Additionally, ZrHx is a corrosion product in Zircaloy. ZrHx exists in multiple phases or polymorphs with varying stoichiometry [2]. The most significant of these phases are the δ phase and ϵ phase, the former of which is the primary constituent of TRIGA fuel and Zircaloy corrosion [3,4]. At room temperature, the δ phase is dominant for 1.59 < x < 1.64 and exists in a mixture with α -Zirconium for x < 1.59 below the β Zr α Zr + δ ZrH x reaction eutectoid temperature of approximately 820 K [2]. The ϵ phase is dominant for x > 1.74 and is mixed with the δ phase for 1.64 < x < 1.74 . This stability range for the ϵ phase is approximately from room temperature up to nearly 800 K; however, the stability range for the δ phase grows considerably with increasing temperature, shrinking the mixed-phase region [2]. As a result of crystal binding of hydrogen in ZrHx, energy exchange associated with thermal neutron scattering involves narrow, well-separated energy levels in contrast to the large continuum spectra observed in other materials [5,6]. This behavior is characteristic of metal-hydrides [5,6] and greatly influences neutron thermalization effects, such as the warm neutron principle responsible for the large negative moderator feedback coefficient of TRIGA reactors.
In neutron transport calculations (e.g., Monte Carlo) of reactors, Evaluated Nuclear Data File (ENDF) thermal scattering law (TSL) sub-libraries are used to capture the effects of crystal binding effects on low-energy neutron scattering. The TSL, S ( α , β ) , is a temperature-dependent probability distribution defining the possible scattered states for an incident neutron, and has contributions from distinct effects (d) due to scattered wave interference between different nuclear scattering sites and self (s) non-interference effects [7],
S ( α , β ) = S d ( α , β ) + S s ( α , β ) .
The momentum and energy transfer resulting from scattering are described by the unitless parameters α and β as,
α = E + E μ E E A k B T
β = E E k B T
with μ and A representing the scattering cosine and nuclide to neutron mass ratio, respectively. The incident and scattered neutron energies are E and E , and k B and T are the Boltzmann constant and temperature.
In crystals, the phonon expansion is used to compute the TSL such that [7,8],
S ( α , β ) = p S p ( α , β ) ,
which is a summation of increasing orders of the convolution of the vibrational (phonon) spectra, with p = 0 representing elastic scattering and p > 0 representing inelastic scattering. In current ENDF thermal scattering sub-libraries, the incoherent approximation is applied, such that distinct effects on inelastic scattering are neglected. Within this approximation, the double-differential scattering cross section is related to the TSL using Fermi’s Golden rule as [7],
2 σ E = 1 4 π k B T E E σ c o h S d 0 ( α , β ) + S s 0 ( α , β ) + p > 0 S s p ( α , β ) + σ i n c S s ( α , β )
where σ c o h and σ i n c are the bound coherent and incoherent nuclear cross sections of the nuclide, and T is the temperature of the material.
The current ENDF/B-VIII.0 ZrHx TSLs were generated in the incoherent approximation within the LEAPR module of NJOY using the vibrational spectra calculated with the central force model developed at General Atomics (GA) [8,9,10], which is based upon the ϵ phase (i.e., ZrH2). Nevertheless, these TSLs have historically been used to represent both phases in reactor physics calculations. Elastic scattering was assumed to be incoherent for both hydrogen and zirconium, such that S d 0 = 0 ; an approximation which is reasonable for hydrogen due to its large incoherent nuclear scattering cross section. This approximation, however, neglects the effects of crystal structure important for scattering from zirconium bound in ZrHx.
In this work, the TSLs for hydrogen and zirconium bound in ZrHx are generated from phonon spectra derived from ab initio lattice dynamics (AILD) and ab initio molecular dynamics (AIMD) for both δ -ZrHx and ϵ -ZrH2, thereby expanding the set of available sub-libraries for this material. These TSLs were generated using the Full Law Analysis Scattering System Hub (FLASSH), which includes a generalized coherent elastic scattering capability [11]. At this stage, coherent elastic scattering only incorporates zirconium, enabling the sub-libraries to be extensible, within the limitations of the ENDF format, to the entire stoichiometry domain of each phase and mixed phase.

2. Computational Methods

The δ and ϵ phases of ZrHx are characterized by cubic and tetragonal crystal structures, respectively [12]. For the purposes of ab initio calculations, the stoichiometry is represented as an idealized x = 1.5 for δ phase and x = 2 for ϵ phase. The unit cells of each phase used in simulation are shown in Figure 1. Initially, the geometry of these structures were optimized through electronic structure calculations using the Vienna ab initio Simulation Package (VASP) [13,14] with plane-augmented wave (PAW) potentials and a GGA-PBE exchange correlation functional [15,16]. A plane wave cutoff of 450 eV and k-mesh of less than 1.8 nm 1 were used and found to be sufficient for energy convergence of less than 4 meV/atom. The predicted lattice parameter of the δ phase are a = 0.477 nm, which is in good agreement with the experimental value of a = 0.478 nm [12]. The predicted ϵ phase parameters of a = 0.354 nm and c = 0.440 nm are also in reasonable agreement with experiment, a = 0.352 nm and c = 0.445 nm [12].
The phonon density of states (DOS), which is the fundamental input to the generation of TSL in the incoherent approximation, was generated in the harmonic approximation for δ -ZrH1.5 and ϵ -ZrH2 using AILD methods involving the coupling of the VASP [13,14] and PHONON [17] codes, described in Ref. [18]. These calculations were performed with the MedeA simulation environment [19]. Hellmann–Feynman forces were extracted from VASP calculations of 2 × 2 × 2 (80 atoms) δ -ZrH1.5 and 3 × 3 × 2 (108 atoms) ϵ -ZrH2 supercells with the previously listed convergence parameters. The forces describing the small displacements of zirconium atoms are expected to be well-characterized within the harmonic approximation; however, hydrogen—due to its small mass and large vibrational displacement—may have significant anharmonic contributions to its motion even at 0 K (i.e., zero-point motion). AIMD simulations, which inherently include anharmonicity, were therefore performed using methods developed in previous work [20,21] to generate the hydrogen DOS. Using a reduced 350 eV planewave cutoff and a 1 × 1 × 1 k-mesh, a 3 × 3 × 3 (270 atoms) δ -ZrH1.5 and 4 × 4 × 3 (288 atoms) ϵ -ZrH2 supercell were first allowed to thermodynamically equilibrate for 1 ps at 0.00025 ps time-steps under a 300 K NVT thermostat. Atomic trajectories were then evolved for an additional 1.5 ps. Subsequently, the velocities at each time-step during this second time-period were used to generate the phonon DOS for hydrogen as the Fourier transform of velocity auto-correlation function, which, at 300 K, captures intrinsic anharmonic effects from the atomic forcefield without introducing significant additional spectral broadening.
The resulting hydrogen phonon DOS from AIMD, as seen in Figure 2, are shifted to lower energies compared to AILD calculations, and are more consistent with inelastic neutron spectroscopy (INS) measurements. Additionally, the AIMD DOS are more analogous to experiments than the GA spectrum used in the current ENDF evaluation [8,10], which has been previously noted to be qualitatively different from the experimental spectra [5,22]. Although the difference between the approximate AILD and AIMD average phonon energies is less than 0.01 eV, the deviation has an increasingly large impact on the predicted TSL and corresponding cross section for neutron scattering events with hydrogen involving the emission of multiple phonons. Specifically, the predicted oscillations using an AILD spectra diverge from experimental observation with increasing phonon order. Consequently, the AIMD phonons for hydrogen are expected to be more suitable for use in the generation of TSL.
TSL for zirconium and hydrogen in δ -ZrHx (Zr(ZrHx) and H(ZrHx)), as well as zirconium and hydrogen in ϵ -ZrH2 (Zr(ZrH2) and H(ZrH2)) were generated in the incoherent approximation as File 7 (MF = 7) thermal scattering sub-libraries (NSUB = 12) in ENDF-6 format [23] using the FLASSH code [11]. Each File 7 included evaluation of the TSL for temperatures between 296 K and 1200 K. The bound and free nuclear scattering cross sections (1 b = 10 28 m2) used for each element are listed in Table 1 and were derived from relevant data in the NIST neutron scattering length database and ENDF/B-VIII.0 neutron reaction sub-libraries (at E = 0.0253 eV from File 3 MT = 2 processed at 0 K) [24,25]. The zirconium File 7 corresponds to the natural isotopic mixture, as reported in IUPAC atomic weight and isotopic abundance evaluations [26,27]. The phonon DOS used for zirconium and hydrogen were from AILD and AIMD simulations, respectively. A phonon order of 300 was used in the phonon expansion of all inelastic reactions (MT = 4).
The coherent elastic scattering due to crystalline structure of each phase was computed using the generalized elastic scattering treatment in FLASSH and stored on the corresponding elastic reaction (MT = 2) of the zirconium File 7. The Debye–Waller matrix in the present work was computed in the cubic approximation. Although hydrogen has an appreciable coherent scattering length ( 3.7 fm) when compared to zirconium ( 7.3 fm) [24,25], the contribution of the hydrogen sublattice was neglected so that the TSL is extensible across the entire stoichiometric range of each phase.
The hydrogen TSLs used the incoherent nuclear scattering option for MT = 2. The nuclear scattering cross section stored on this reaction is only the nuclear incoherent scattering cross section, whereas the total free-atom cross section is used for MT = 4 per Equation (5). This treatment, while in contrast to the use of the total cross section for MT = 2 in ENDF/B-VIII.0, is both consistent with the neglection of the hydrogen contribution to coherent elastic scattering and the treatment of nitrogen in uranium mononitride [25], which like hydrogen, has a non-negligible coherent and incoherent nuclear scattering cross section.

3. Thermal Scattering Laws

Inelastic neutron scattering techniques provide a direct measure of the double-differential cross section of a compound, which is also directly calculable from the TSL via Equation (5). The double-differential scattering cross section for ZrH2 generated in the present work at the Naval Nuclear Laboratory is compared to the experiment and ENDF/B-VIII.0 in Figure 3. The anticipated oscillatory behavior due to the localized, higher-energy hydrogen phonon DOS is observed, whereby peaks occur at intervals of approximately 0.14 eV of energy loss. Both the Naval Nuclear Laboratory and current ENDF TSLs demonstrate reasonable agreement with experiment [28]. Despite the similarity in the energy structure of the present and ENDF/B-VIII.0 double-differential cross sections, a slight divergence occurs in each subsequently increasing energy loss peak due to differences in the average hydrogen phonon energy; however, as noted previously, the Naval Nuclear Laboratory phonon DOS more closely resembles the experiment.
The differential scattering cross section of ϵ -ZrH1.92, shown in Figure 3, is computed as an integral of Equation (5) over all scattered neutron energies. The Naval Nuclear Laboratory TSL and ENDF TSL provide similar predictions and trend well with experiment [29]. These calculations exclude the contribution of coherent elastic scattering from the differential cross section; therefore, the locations and relative amplitudes of the coherent elastic reaction (i.e., diffraction pattern) of the Naval Nuclear Laboratory zirconium TSL are also shown. Deviation from the experiment in the calculated differential cross sections at low scattering cosine are due to the use of the incoherent nuclear scattering cross section as opposed to the total elastic reaction, MT = 2, of the Naval Nuclear Laboratory TSLs. At low-incident neutron energies where coherent elastic scattering is significant, diffraction peaks occur in the experimental data that result in a deviation from the continuous differential behavior of the calculations. Overall, the calculated diffraction pattern matches these peaks, thereby demonstrating the importance of this contribution for the prediction of angular distributions.
The inelastic contributions to the Naval Nuclear Laboratory TSLs for both δ -ZrHx (H(ZrHx), Zr(ZrHx)) and ϵ -ZrH2 (H(ZrH2), and Zr(ZrH2)) are illustrated in Figure 3 at temperatures of 296 K and 700 K, where the latter is representative of the maximum fuel temperature of a 2000 MW pulse in a TRIGA reactor [3]. In general, the TSLs for hydrogen and zirconium are similar between the two phases. The zirconium TSLs indicate a quasi-elastic scattering spectrum at low momentum transfers, trending toward free gas behavior at large momentum transfer (e.g., α = 1.979 ). Oscillations in the hydrogen TSLs strongly favor energy loss ( β < 0 ) at 296 K, but have a significant contribution from 0.14 eV 1-phonon up-scattering events at 700 K. This behavior is a contributor to the strong negative moderator feedback coefficient of TRIGA fuel systems. At low momentum transfers (e.g., α = 0.4943 ), 1-phonon events are dominant; however, as the momentum transfer increases the possibility of down-scattering, involving multiple phonons becomes likely, exemplified by α = 5.067 . However, the multi-phonon peaks occurring at integer values of approximately 0.14 eV are similar for the two phases, and oscillations diverge with increasing | β | due to the small increase in the average phonon energy of δ relative to the ϵ phase, that is also observed in the experiment [5].
Integrated total scattering cross sections for δ -ZrHx and ϵ -ZrH2 were generated using the NDEX nuclear data processing package [30] and are shown in Figure 4. NDEX uses the β grid tabulated on the TSL and an adaptive energy mesh to provide high-resolution cross sections [31], which is important for predicting the oscillator behavior of metal hydride moderators that is not fully captured with the coarse, fixed energy meshes used in other codes (e.g., NJOY [9]), as demonstrated in previous work [31,32]. Several compositions were selected for comparison between the two models and experiment, with x = 1.6 for the δ phase and x = 2 for ϵ phase representing TRIGA and SNAP-10A fuel, respectively [1,3]. Each cross section is illustrated at room and elevated temperatures, with 700 K corresponding to pulse characteristics of a TRIGA and 800 K corresponding to the approximate maximum operating temperature of SNAP-10A [1].
At 296 K, the impact of changes in the coherent and incoherent elastic scattering is appreciable for both phases, resulting in a 6–10 b difference in the total scattering cross section below the Bragg edge between the Naval Nuclear Laboratory and current ENDF evaluations. In particular, the Bragg peaks at 0.0027 and 0.0036 eV in the δ -ZrH Naval Nuclear Laboratory evaluation are in reasonable agreement with experimental data [3], but are absent from the current ENDF evaluation. Below the Bragg cut-off (0.0027 eV), the experimental cross section is in between the Naval Nuclear Laboratory and ENDF calculations, consistent with the neglect of coherent elastic effects of hydrogen (e.g., correlated hydrogen vacant sites). Additionally, ZrH1.6 shows a parallel 1 / v trend between the Naval Nuclear Laboratory and ENDF evaluations in this region, whereas ZrH2 diverges with less pronounced 1 / v behavior in the Naval Nuclear Laboratory evaluation. This disparity arises from a decreased contribution of low-energy (<0.03 eV) phonons in the ϵ phase compared to the δ phase, present in both AILD and AIMD calculations.
As temperature increases, inelastic scattering with hydrogen begins to dominate, and the 0.14 eV 1-phonon up-scatter becomes the primary contributor to the 1 / v scattering cross section. Consequently, the total scattering cross sections for the Naval Nuclear Laboratory and ENDF evaluations are nearly the same at both 700 K and 800 K; however, an observable consequence of coherent elastic scattering on the total scattering cross section persists between approximately 0.001–0.004 eV. At incident energies above 0.05 eV, both calculated cross sections trend well with the experiment. In this range, inelastic scattering causes oscillations in the cross section at intervals of approximately 0.14 eV that broaden with increasing temperature, and quantitative differences between Naval Nuclear Laboratory and ENDF evaluations arise primarily from a 0.4 b difference in the free-atom cross section used in the generation of the zirconium TSLs.

4. Conclusions

Phonon spectra for hydrogen and zirconium in δ and ϵ zirconium hydride were developed using predictive AIMD and AILD simulations. Subsequently, these spectra were used in the FLASSH code to generate new TSL evaluations for zirconium hydride which include, for the first time, the previously neglected contribution of coherent elastic scattering in multiple phases of this material. These evaluations include both δ -ZrHx (H(ZrHx) and Zr(ZrHx)) and ϵ -ZrH2 (H(ZrH2) and Zr(ZrH2)), permitting explicit modeling of each phase in reactor physics calculations or benchmarks (e.g., TRIGA, SNAP-10A) across the entire range of their stoichiometry. Currently, only a single sub-library set (H(ZrH) and Zr(ZrH)) is included in the ENDF evaluation for this material, such that the new Naval Nuclear Laboratory evaluations may be considered not only a re-evaluation, but also an expansion of the ENDF/B library database. Moreover, differences in both the coherent elastic reaction and hydrogen phonon spectra between the evaluations provides the capability to study the importance of separate treatment for the δ and ϵ phases in neutronics calculations.

Author Contributions

Conceptualization, J.W. and M.Z.; methodology, J.W., M.Z. and J.H.; software, J.W.; validation, J.W., M.Z. and J.H.; formal analysis, J.W.; investigation, J.W.; writing—original draft preparation, J.W.; writing—review and editing, M.Z.; visualization, J.W.; funding acquisition, M.Z.; supervision, M.Z.; project administration, M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was performed under contract No. DOE-89233018CNR000004.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The submitted manuscript has been authored by contractors of the US Government under contract No. DOE-89233018CNR000004. Accordingly, the US Government retains a non-exclusive, royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for US Government purposes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lillie, A.F.; McClelland, D.T.; Roberts, W.J.; Walter, J.H. Zirconium Hydride Fuel Element Performance Characteristics; AI-AEC-13084; Atomics International Division: Canoga Park, CA, USA, 1973. [Google Scholar]
  2. Zuzek, E.; Abriata, J.P.; San-Martin, A.; Manchester, F.D. The H-Zr (Hydrogen-Zirconium) System. Bull. Alloy Phase Diagrams 1990, 11, 385–395. [Google Scholar] [CrossRef]
  3. West, G.B.; Whittemore, W.L.; Shoptaugh, J.R.; Dee, J.B.; Coffer, C.O. Kinetic Behavior of TRIGA Reactors; GA-7882; General Atomics: San Diego, CA, USA, 1967. [Google Scholar]
  4. Bair, J.; Zaeem, M.A.; Tonks, M. A Review of Hydride Precipitation in Zirconium Alloys. J. Nucl. Mater. 2015, 446, 12–20. [Google Scholar] [CrossRef]
  5. Couch, J.G.; Harling, O.K.; Clune, L.C. Structure in the Neutron Scattering Spectra of Zirconium Hydride. Phys. Rev. B 1971, 4, 2675–2681. [Google Scholar] [CrossRef]
  6. Whittemore, W.L. Neutron Interactions in Zirconium Hydride; GA-4490; General Atomics: San Diego, CA, USA, 1964. [Google Scholar]
  7. Hawari, A.I. Modern Techniques for Inelastic Thermal Neutron Scattering Analysis. Nucl. Data Sheets 2014, 118, 172–175. [Google Scholar] [CrossRef]
  8. MacFarlane, R.E. New Thermal Neutron Scattering Files for ENDF/B-VI Release 2; LA-12639-MS; Los Alamos National Lab: Los Alamos, NM, USA, 1994. [Google Scholar]
  9. Kahler, A.C.; MarFarlane, R.E. The NJOY Nuclear Data Processing System, Version 2012; LA-UR-12-27079; Los Alamos National Lab: Los Alamos, NM, USA, 2012. [Google Scholar]
  10. Slaggie, E.L. Central Force Lattice Dynamical Model for Zirconium Hydride. J. Phys. Chem. Solids 1968, 29, 923–934. [Google Scholar] [CrossRef]
  11. Zhu, Y.; Hawari, A.I. Full Law Analysis Scattering System Hub (FLASSH). In Proceedings of the PHYSOR 2018, Cancun, Mexico, 22–26 April 2018. [Google Scholar]
  12. Kempter, C.P.; Elliot, R.O.; Gschneider, K.A., Jr. Thermal Expansion of Delta and Epsilon Zirconium Hydrides. J. Chem. Phys. 1960, 33, 837–840. [Google Scholar] [CrossRef]
  13. Kresse, G.; Furthmüller, J. Efficiency of ab initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  14. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total Energy Calc. Using A Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  15. Blöchl, P.E. Projector Augmented-wave Method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [Green Version]
  16. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  17. Parlinski, K.; Li, Z.Q.; Kawazoe, Y. First-principles Determination of the Soft Mode in Cubic ZrO2. Phys. Rev. Lett. 1997, 78, 4063–4066. [Google Scholar] [CrossRef]
  18. Hawari, A.I.; Al-Qasir, I.I.; Gillette, V.H.; Wehring, B.W. Ab initio Generation of Thermal Neutron Scattering cross sections. In Proceedings of the PHYSOR 2004, Chicago, IL, USA, 25–29 April 2004. [Google Scholar]
  19. MedeA Version 3.0; MedeA: San Diego, CA, USA, 2019.
  20. Wormald, J.L.; Hawari, A.I. Generation of Phonon Density of States and Thermal Scattering Law Using Ab Initio Molecular Dynamics. Prog. Nucl. Energy 2017, 101, 461–467. [Google Scholar] [CrossRef]
  21. Wormald, J.L.; Hawari, A.I. Thermal Neutron Scattering Law Calculations using Ab Initio molecular Dynamics. EPJ Web Conf. 2017, 146, 13002. [Google Scholar] [CrossRef] [Green Version]
  22. Malik, S.S.; Rorer, D.C.; Brunhart, G. Optical-phonon Structure and Precision Neutron Total cross section Measurements of Zirconium Hydride. J. Phys. F Met Phys. 1984, 14, 73–81. [Google Scholar] [CrossRef]
  23. Herman, M.; Trkov, A. Data Formats and Procedures for the Evaluated Nuclear Data File ENDF/B-VI and ENDF/B-VII; BNL-90365-2009 (ENDF-102); Brookhaven National Laboratory: Upton, NY, USA, 2009. [Google Scholar]
  24. Sears, V.F. Neutron Scattering Lengths and Cross Sections. Neutron News 1992, 3, 26–37. [Google Scholar] [CrossRef]
  25. Brown, D.A.; Chadwick, M.B.; Capote, R.; Kahler, A.C.; Trkov, A.; Herman, M.W.; Sonzogni, A.A.; Danon, Y.; Carlson, A.D.; Dunn, M.; et al. ENDF/B-VIII.0: The 8th Major Release of the Nuclear Reaction Data Library with CIELO-project Cross Sections, New Standards and Thermal Scattering Data. Nucl. Data Sheets 2018, 148, 1–142. [Google Scholar] [CrossRef]
  26. Coursey, J.S.; Schwab, D.J.; Tsai, J.J.; Dragoset, R.A. Atomic Weights and Isotopic Compositions (Version 4.1); National Institute of Standards and Technology: Gaithersburg, MD, USA, 2015. Available online: http://physics.nist.gov/Comp (accessed on 14 April 2019).
  27. Berglund, M.; Wieser, M.E. Isotopic Compositions of the Elements 2009 (IUPAC Technical Report). Pure Appl. Chem. 2011, 83, 397–410. [Google Scholar] [CrossRef] [Green Version]
  28. Purohit, S.N.; Pan, S.S.; Bischoff, F.; Bryant, W.A.; Lajeunesse, C.; Yeater, M.L.; Moore, W.E.; Kirouac, G.J.; Esch, L.J.; Francis, N.C. Inelastic Neutron Scattering in Metal Hydrides, UC and UO2, and Applications of the Scattering Law. In Proceedings of the Symposium on Neutron Thermalization and Reactor Spectra, Ann Arbor, MI, USA, 17–21 July 1967; Volume I, pp. 407–433. [Google Scholar]
  29. Kornbichler, S. Bestimmung der Diffusionskonstanten D(E0,T) und <Dv(T)> Thermischer Neutronen in H2O, Phenylen, ZrH1.92 und D2O durch Messung der Streuwinkelverteilungen. II. Zirkonhydrid und Schweres Wasser. Nukleonik 1965, 78, 281. [Google Scholar]
  30. Griesheimer, D.P.; Gill, D.F.; Nease, B.R.; Sutton, T.M.; Stedry, M.H.; Dobreff, P.S.; Carpenter, D.C.; Trumbull, T.H.; Caro, E.; Joo, H.; et al. MC21 v6.0—A continuous-energy Monte Carlo Particle Transport Code with Integrated Reactor Feedback Capabilities. Ann. Nucl. Energy 2015, 82, 29. [Google Scholar] [CrossRef]
  31. Wormald, J.L.; Thopmson, J.T.; Trumbull, T.H. Implementation of an Adaptive Energy Grid Routine in NDEX for the Processing of Thermal Neutron Scattering Cross Sections. Ann. Nucl. Energy 2020, 149, 107773. [Google Scholar] [CrossRef]
  32. Zerkle, M.; Holmes, J. A Thermal Neutron Scattering Law for Yttrium Hydride. EPJ Web Conf. 2017, 146, 13005. [Google Scholar] [CrossRef] [Green Version]
  33. Schmidt, U. Untersuchung der Wasserstoff- und Deuteriumschwingungen in Metallhydriden, -Hydrodeuteriden und -Deuteriden Mittles Totaler Neutronen- Wirkunsquerschnitte. ATKE 1967, 12, 385. [Google Scholar]
Figure 1. Crystal unit cell for ϵ -ZrH2 (left) and δ -ZrH1.5 (right). Zirconium atoms are colored green (large) and hydrogen atoms are rose (small).
Figure 1. Crystal unit cell for ϵ -ZrH2 (left) and δ -ZrH1.5 (right). Zirconium atoms are colored green (large) and hydrogen atoms are rose (small).
Jne 02 00011 g001
Figure 2. Phonon DOS for hydrogen in ϵ -ZrH2 (left) and δ -ZrH1.5 (right) compared to INS [5]. The experimental, AILD, and AIMD phonon spectra differ between phases. The AIMD spectra is used for subsequent evaluation of the hydrogen TSLs in this work.
Figure 2. Phonon DOS for hydrogen in ϵ -ZrH2 (left) and δ -ZrH1.5 (right) compared to INS [5]. The experimental, AILD, and AIMD phonon spectra differ between phases. The AIMD spectra is used for subsequent evaluation of the hydrogen TSLs in this work.
Jne 02 00011 g002
Figure 3. Computed double-differential cross section of ϵ -ZrH2 for E = 0.572 eV and μ = 0.766 (top left) and differential cross section of ϵ -ZrH1.92 for E = 0.0225 eV (top right) compared to the experiment, Purohit [28] and Korbichler [29]. TSL for H (bottom left), and Zr (bottom right) in δ -ZrHx and ϵ -ZrH2 are shown at 296 K and 700 K. Both α and β are referenced to 293.6 K (i.e., k B T = 0.0253 eV). The double-differential scattering cross sections for the present Naval Nuclear Laboratory (NNL) and ENDF/B-VIII.0 TSL evaluations were generated with the FLASSH code [11]. The diffraction peaks in the differential scattering cross section indicate the location of contributions from coherent elastic scattering in the NNL evaluation that are not included in the NNL plot. Diffraction presents as a deviation from continuous behavior in the experimental differential cross section. Experimental double-differential cross sections were obtained from the time of flight measurements.
Figure 3. Computed double-differential cross section of ϵ -ZrH2 for E = 0.572 eV and μ = 0.766 (top left) and differential cross section of ϵ -ZrH1.92 for E = 0.0225 eV (top right) compared to the experiment, Purohit [28] and Korbichler [29]. TSL for H (bottom left), and Zr (bottom right) in δ -ZrHx and ϵ -ZrH2 are shown at 296 K and 700 K. Both α and β are referenced to 293.6 K (i.e., k B T = 0.0253 eV). The double-differential scattering cross sections for the present Naval Nuclear Laboratory (NNL) and ENDF/B-VIII.0 TSL evaluations were generated with the FLASSH code [11]. The diffraction peaks in the differential scattering cross section indicate the location of contributions from coherent elastic scattering in the NNL evaluation that are not included in the NNL plot. Diffraction presents as a deviation from continuous behavior in the experimental differential cross section. Experimental double-differential cross sections were obtained from the time of flight measurements.
Jne 02 00011 g003
Figure 4. Scattering cross section of δ -ZrH compared to experiments at room temperature, Whittemore [6] and Schmidt [33] (left). Scattering cross sections of δ -ZrH1.6 and ϵ -ZrH2 compared to current ENDF/B-VIII.0 evaluations (right). The integrated cross section for both the NNL TSL evaluation and ENDF/B-VIII.0 were generated with the NDEX nuclear data processing code [30,31].
Figure 4. Scattering cross section of δ -ZrH compared to experiments at room temperature, Whittemore [6] and Schmidt [33] (left). Scattering cross sections of δ -ZrH1.6 and ϵ -ZrH2 compared to current ENDF/B-VIII.0 evaluations (right). The integrated cross section for both the NNL TSL evaluation and ENDF/B-VIII.0 were generated with the NDEX nuclear data processing code [30,31].
Jne 02 00011 g004
Table 1. Elemental scattering cross sections used in the present TSL evaluations for ZrHx. The coherent and incoherent cross sections are bound. Each free atom cross section represents an isotopic average of the ENDF/B-VIII.0 neutron reaction sub-libraries at 0.0253 eV [25]. For hydrogen, the incoherent cross section corresponds to the Sears database [24], whereas the coherent scattering cross section is estimated as the difference between the bound total and bound incoherent scattering cross sections.
Table 1. Elemental scattering cross sections used in the present TSL evaluations for ZrHx. The coherent and incoherent cross sections are bound. Each free atom cross section represents an isotopic average of the ENDF/B-VIII.0 neutron reaction sub-libraries at 0.0253 eV [25]. For hydrogen, the incoherent cross section corresponds to the Sears database [24], whereas the coherent scattering cross section is estimated as the difference between the bound total and bound incoherent scattering cross sections.
ElementA σ coh (b) σ inc (b) σ free (b)
Zr90.43996.7809906.633485
H0.9991671.758380.0541620.43608
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wormald, J.; Zerkle, M.; Holmes, J. Generation of the TSL for Zirconium Hydrides from Ab Initio Methods. J. Nucl. Eng. 2021, 2, 105-113. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020011

AMA Style

Wormald J, Zerkle M, Holmes J. Generation of the TSL for Zirconium Hydrides from Ab Initio Methods. Journal of Nuclear Engineering. 2021; 2(2):105-113. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020011

Chicago/Turabian Style

Wormald, Jonathan, Michael Zerkle, and Jesse Holmes. 2021. "Generation of the TSL for Zirconium Hydrides from Ab Initio Methods" Journal of Nuclear Engineering 2, no. 2: 105-113. https://0-doi-org.brum.beds.ac.uk/10.3390/jne2020011

Article Metrics

Back to TopTop