@article{mahbuba_beeler_jokisaari_2025, title={Grain boundary self-diffusion and point defect interactions in α-U via molecular dynamics}, url={https://doi.org/10.1016/j.jnucmat.2024.155521}, DOI={10.1016/j.jnucmat.2024.155521}, journal={Journal of Nuclear Materials}, author={Mahbuba, Khadija and Beeler, Benjamin and Jokisaari, Andrea}, year={2025}, month={Jan} } @article{mao_li_park_beeler_hu_2025, title={A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity}, volume={246}, ISSN={["1879-0801"]}, url={https://doi.org/10.1016/j.commatsci.2024.113371}, DOI={10.1016/j.commatsci.2024.113371}, journal={COMPUTATIONAL MATERIALS SCIENCE}, author={Mao, Zirui and Li, Yulan and Park, Gyuchul and Beeler, Benjamin and Hu, Shenyang}, year={2025}, month={Jan} } @article{kadambi_aagesen_zhang_beeler_2024, title={Assessment of effective elastic constants of U-10Mo fuel: A multiscale modeling and homogenization study}, volume={599}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2024.155225}, abstractNote={The significant microstructural changes that U-Mo fuel undergoes during operation degrades its mechanical properties and structural integrity. Microstructural evolution entails the formation, evolution, and redistribution of porosity in conjunction with grain refinement. In the present paper, we employ numerical approaches to assess the impact of the various microstructural features—grains, nanoscale intragranular fission gas bubbles, and mesoscale intergranular voids—on the degradation of elastic constants. Phase-field microstructure models are combined with the asymptotic expansion homogenization technique in order to derive the effective elastic constants as a function of porosity and fission density. The results are verified and compared against theoretical bounds. Using this approach, elastic degradation in operating nuclear fuels can be quantified when the distributions of microstructural features are known.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Kadambi, Sourabh B. and Aagesen, Larry K. and Zhang, Yongfeng and Beeler, Benjamin}, year={2024}, month={Oct} } @article{abdulhameed_beeler_galvin_cooper_2024, title={Assessment of uranium nitride interatomic potentials}, volume={60}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2024.155247}, DOI={10.1016/j.jnucmat.2024.155247}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={AbdulHameed, Mohamed and Beeler, Benjamin and Galvin, Conor O. T. and Cooper, Michael W. D.}, year={2024}, month={Nov} } @article{hasan_beeler_2024, title={Calculation of grain boundary diffusion coefficients in γ U-Mo using atomistic simulations}, volume={598}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2024.155190}, DOI={10.1016/j.jnucmat.2024.155190}, abstractNote={The γU-Mo alloy has been selected for the conversion of U.S. High-Performance Research Reactor (HPRR) fuel from highly enriched uranium to low enriched uranium as a part of the effort to reduce nuclear proliferation risks. Although γU-Mo alloys have the advantage of high uranium density and good overall irradiation performance, the irradiation-induced swelling and creep in the metal remain important design parameters since they influence the mechanical and thermal integrity of the fuel. To account for these design criteria, engineering scale models need fundamental properties, such as diffusion coefficients, as input. In this study, the diffusion of related species along grain boundaries in γU-Mo fuel is quantified considering that grain boundaries act as sinks for point defects, nucleation sites for gas bubbles, avenues for Coble creep, etc. The diffusivities of U and Mo in selected grain boundaries of γU-Mo alloys (γU-7Mo, γU-10Mo, and γU-12Mo) are obtained utilizing molecular dynamics simulations for a temperature range of 600 K - 1200 K with an interval of 100 K. The structures analyzed include symmetric tilt, asymmetric tilt, and twist grain boundaries. The grain boundary diffusion coefficients of U and Mo in the examined γU-Mo alloys are on the order of 10−14 to 10−11 m2 s−1. It is observed that the U diffusivity in the grain boundary is higher than the Mo diffusivity in all cases and that the increase in Mo content of the alloy correlates to a decrease in the grain boundary diffusion. Xe diffusion along γU-10Mo grain boundaries is also calculated in this work, and the diffusivity of Xe in the γU-10Mo grain boundaries is found to be 8 to 15 orders of magnitude higher than the intrinsic Xe diffusivity in γU-10Mo depending on the temperature. The information gathered in this work can inform fuel swelling and creep models and help understand various other phenomena related to γU-Mo fuel performance.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Hasan, A. T. M. Jahid and Beeler, Benjamin}, year={2024}, month={Sep} } @article{shousha_beeler_aagesen_beausoleil ii_okuniewski_2024, title={First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport}, volume={601}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2024.155310}, DOI={10.1016/j.jnucmat.2024.155310}, abstractNote={The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. The calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Shousha, Shehab and Beeler, Benjamin and Aagesen, Larry K. and Beausoleil II, Geoffrey L. and Okuniewski, Maria A.}, year={2024}, month={Dec} } @article{duemmler_andersson_beeler_2024, title={First-principles investigation of the thermophysical properties of NaCl, PuCl3, and NaCl-PuCl3 Molten salts}, volume={591}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2024.154902}, DOI={10.1016/j.jnucmat.2024.154902}, abstractNote={Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. The heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl3.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Duemmler, Kai and Andersson, David and Beeler, Benjamin}, year={2024}, month={Apr} } @article{andersson_wang_yang_beeler_2024, title={KCl-UCl3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models}, volume={599}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2024.155226}, abstractNote={Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl3 to form network structures.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Andersson, D. A. and Wang, G. and Yang, P. and Beeler, B. W.}, year={2024}, month={Oct} } @article{shousha_beeler_2024, title={Magnetism and finite-temperature effects in UZr 2: A density functional theory analysis}, volume={595}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2024.155037}, DOI={10.1016/j.jnucmat.2024.155037}, abstractNote={The structure and the thermophysical properties of δ-UZr2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Shousha, Shehab and Beeler, Benjamin}, year={2024}, month={Jul} } @article{tranchida_nicaud_beeler_bourasseau_2024, title={Thermophysical properties and unexpected viscosity of liquid (U, Zr): An atomistic investigation}, volume={160}, ISSN={["1089-7690"]}, url={https://doi.org/10.1063/5.0203177}, DOI={10.1063/5.0203177}, abstractNote={In this study, we performed a numerical investigation of the thermophysical properties of liquid (U, Zr) mixtures, which are particularly relevant in the context of hypothetical nuclear accidents and the formation of in-vessel coriums. To do so, atomistic simulations leveraging classical molecular dynamics and an interatomic potential developed for solid (U, Zr) structures are performed. Our methodology is first validated by comparing the predictions of our model for the melting temperature and the structure factors to experimental, phase diagram, and ab initio data. We then use the approach to evaluate the temperature and composition dependence of four fundamental properties in the context of coriums: density, heat capacity, compressibility, and viscosity. Systematic comparisons to the existing experimental data are performed and discussed. In particular, the viscosity of liquid (U, Zr) mixtures is investigated by comparing diffusion calculations and the Stokes–Einstein formula as well as the results obtained with the Green–Kubo methodology, empirical predictions, and experimental data. Notably, the viscosity of the mixtures is predicted to be significantly higher than that of the single-element liquids, which is unexpected and could have crucial consequences on the early stages of the formation and flow of in-vessel corium.}, number={21}, journal={JOURNAL OF CHEMICAL PHYSICS}, author={Tranchida, J. and Nicaud, F. and Beeler, B. W. and Bourasseau, E.}, year={2024}, month={Jun} } @article{shousha_kadambi_beeler_kombaiah_2024, title={Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study}, volume={8}, ISSN={["2475-9953"]}, url={https://doi.org/10.1103/PhysRevMaterials.8.083605}, DOI={10.1103/PhysRevMaterials.8.083605}, abstractNote={The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. The transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.}, number={8}, journal={PHYSICAL REVIEW MATERIALS}, author={Shousha, Shehab and Kadambi, Sourabh Bhagwan and Beeler, Benjamin and Kombaiah, Boopathy}, year={2024}, month={Aug} } @article{heyl_beeler_2023, title={An ab initio molecular dynamics study of varied compositions of the LiF-NaF-KF molten salt}, volume={585}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2023.154641}, DOI={10.1016/j.jnucmat.2023.154641}, abstractNote={With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. This work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Heyl, Veronica and Beeler, Benjamin}, year={2023}, month={Nov} } @article{wang_beeler_jokisaari_2023, title={An atomistic study of fundamental bulk and defect properties in alpha-uranium}, volume={576}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2023.154289}, abstractNote={Alpha-uranium plays an important role in the performance and structure evolution of metallic fuels under irradiation. It is a highly anisotropic material that demonstrates a complex microstructural response to irradiation dependent on the irradiation temperature, and the mechanisms that control this behavior are not well understood. In this work, fundamental bulk properties and energetic and thermodynamic properties of single point defects and di-defects in α-uranium are determined using molecular dynamics over the range of 400 K - 750 K. The constant pressure heat capacity, total volumetric thermal expansion, anisotropy of the thermal expansion, and single vacancy and single interstitial formation energies all agree well with previous experimental and ab initio results, lending support to the di-defect formation and binding energies and the diffusivity results. We find that the diffusion of vacancies and di-vacancies is strongly anisotropic, while interstitial and di-interstitial diffusion is slightly anisotropic, and the most rapidly migrating species changes depending on the temperature. The formation energy of vacancies and di-vacancies increases slightly with temperature, while the formation energy of interstitials and di-interstitials increases strongly with temperature. In addition, di-vacancies are very loosely bound over the entire temperature range, with an average binding energy of 0.017 eV. Conversely, di-interstitials are strongly bound (0.64 eV) until approximately 600 K, above which the binding energy decreases significantly, with an average binding energy of 0.31 eV. These results provide valuable new insights into the possible mechanisms of the complex irradiation damage behavior of α-uranium.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Wang, Yuhao and Beeler, Benjamin and Jokisaari, Andrea}, year={2023}, month={Apr} } @article{takasugi_aly_holler_abarca_beeler_avramova_ivanov_2023, title={Development of an efficient and improved core thermal-hydraulics predictive capability for fast reactors: Summary of research and development activities at the North Carolina state University}, volume={412}, ISSN={["1872-759X"]}, DOI={10.1016/j.nucengdes.2023.112474}, abstractNote={The improved understanding of the safety, technical gaps, and major uncertainties of advanced fast reactors will result in designing their safe and economical operation. This paper focuses on the development of efficient and improved core thermal-hydraulics predictive capabilities for fast reactor modeling and simulation at the North Carolina State University. The described research and development activities include applying results of high-fidelity thermal-hydraulic simulations to inform the improved use of lower-order models within fast-running design and safety analysis tools to predict improved estimates of local safety parameters for efficient evaluation of realistic safety margins for fast reactors. The above-described high-to-low model information improvements are being verified and validated using benchmarks such as the OECD/NRC Liquid Metal Fast Reactor Core Thermal-Hydraulic Benchmark and code-to-code comparisons.}, journal={NUCLEAR ENGINEERING AND DESIGN}, author={Takasugi, C. and Aly, A. and Holler, D. and Abarca, A. and Beeler, B. and Avramova, M. and Ivanov, K.}, year={2023}, month={Oct} } @article{duemmler_woods_karlsson_gakhar_beeler_2023, title={First-principles-derived transport properties of Molten chloride salts}, volume={585}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2023.154601}, DOI={10.1016/j.jnucmat.2023.154601}, abstractNote={Molten salts have many applications ranging from a heat transfer medium in both generation IV nuclear reactor designs and the solar industry to thermal storage systems. While molten salts show promising properties for these applications, there still exists a knowledge gap for the transport properties of molten salts at elevated temperatures. This work uses ab initio Molecular Dynamics to investigate the transport properties of KCl, LiCl, KCl-LiCl eutectic, NaCl, MgCl2, and NaCl-MgCl2 eutectic molten salt systems. The properties presented here are the diffusion coefficient, viscosity, and isochoric heat capacity. These properties are compared to experimental data where available and other computational work in cases where no experimental data is available. This is the first work to explore timescales over 100 ps via AIMD for the determination of transport properties in molten salts.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Duemmler, Kai and Woods, Michael and Karlsson, Toni and Gakhar, Ruchi and Beeler, Benjamin}, year={2023}, month={Nov} } @article{ye_oaks_hu_beeler_rest_mei_yacout_2023, title={Integrated simulation of U-10Mo monolithic fuel swelling behavior}, volume={583}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2023.154542}, DOI={10.1016/j.jnucmat.2023.154542}, abstractNote={A separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to consider various processes occurring in RTR fuels during irradiation. To accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions versus fission density) predicted using the phase-field method. The remainder of the fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration, detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Ye, Bei and Oaks, Aaron and Hu, Shenyang and Beeler, Benjamin and Rest, Jeff and Mei, Zhi-Gang and Yacout, Abdellatif}, year={2023}, month={Sep} } @article{andersson_beeler_2022, title={Ab initio molecular dynamics (AIMD) simulations of NaCl, UCl3 and NaCl-UCl3 molten salts}, volume={568}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2022.153836}, DOI={10.1016/j.jnucmat.2022.153836}, abstractNote={Ab initio molecular dynamics (AIMD) simulations are used to calculate select thermophysical and thermodynamic properties of NaCl, UCl3 and NaCl-UCl3 molten salts. Following established approaches, the AIMD simulations include a model for Van der Waals interactions. The Langreth & Lundqvist (vDW-DF2), DFT-D3, and density-dependent energy correction (DFT-dDsC) dispersion models are first tested for molten NaCl in order to assess their accuracy for density and heat capacity predictions across a range of temperatures. Based on the NaCl results, the vdW-DF2 and DFT-dDsC methods are extended to the UCl3 system and compared to available experimental data. Next, mixtures of NaCl-UCl3 are investigated and analyzed with respect to the deviation from ideal solution behavior. For the DFT-dDsC simulations, density deviates by up to 2% from an ideal mixture, with the maximum occurring close to the eutectic composition. The mixing energy also deviates from an ideal solution and exhibits a minimum of about -0.074 to eV per formula unit, again close to the eutectic composition. Finally, the compressibility and species diffusivity of the pure and mixed salt systems are calculated. The diffusivities are slightly reduced in the mixed compared to the pure systems and the compressibilities loosely follow a linear correlation as a function of the UCl3 composition. The trends observed for mixing properties are rationalized by correlating them to the coordination chemistry, which emphasizes the importance of maintaining the extended network formed by U and Cl ions as the NaCl concentration increases. The concentration at which breakdown of the extended network occurs, roughly coincides with the minimum of both the mixing energy and the deviation from ideal solution behavior for density.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Andersson, D. A. and Beeler, B. W.}, year={2022}, month={Sep} } @article{aly_beeler_avramova_2022, title={Ab initio molecular dynamics investigation of gamma-(U,Zr) structural and thermal properties as a function of temperature and composition}, volume={561}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2022.153523}, abstractNote={Uranium in its metallic form is considered as a fuel for sodium fast reactors due to its higher thermal conductivity and high fissile material density relative to UO 2 fuel. The metal is alloyed with zirconium to increase its stability at high temperatures and increase its solidus temperature . This work uses ab initio molecular dynamics to perform an evaluation of the mechanical and thermophysical properties of the γ -(U,Zr) system at temperatures between 1000 K and 1400 K. Among these properties are the equilibrium volume, bulk modulus , molar heat capacity , heat of formation , and the surface energy. The obtained results are compared to experimental data and previous computational work available in the literature. This is the first study of γ -(U,Zr) utilizing ab initio molecular dynamics, and reduces thermophysical property knowledge gaps that are currently present in the literature.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Aly, Ahmed and Beeler, Benjamin and Avramova, Maria}, year={2022}, month={Apr} } @article{duemmler_woods_karlsson_gakhar_beeler_2022, title={An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl2}, volume={570}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2022.153916}, DOI={10.1016/j.jnucmat.2022.153916}, abstractNote={Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. A wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl2. Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Duemmler, Kai and Woods, Michael and Karlsson, Toni and Gakhar, Ruchi and Beeler, Benjamin}, year={2022}, month={Nov} } @article{beeler_zhang_hasan_park_hu_mei_2022, title={Analyzing the effect of pressure on the properties of point defects in ?U-Mo through atomistic simulations (vol 15, pg 874. 2022)}, volume={12}, ISSN={["2059-8521"]}, url={https://doi.org/10.1557/s43580-022-00436-7}, DOI={10.1557/s43580-022-00436-7}, journal={MRS ADVANCES}, author={Beeler, Benjamin and Zhang, Yongfeng and Hasan, A. T. M. Jahid and Park, Gyuchul and Hu, Shenyang and Mei, Zhi-Gang}, year={2022}, month={Dec} } @article{beeler_zhang_hasan_park_hu_mei_2022, title={Analyzing the effect of pressure on the properties of point defects in gamma U-Mo through atomistic simulations}, volume={10}, ISSN={["2059-8521"]}, url={https://doi.org/10.1557/s43580-022-00350-y}, DOI={10.1557/s43580-022-00350-y}, abstractNote={Abstract Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. This work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels. Graphical abstract}, journal={MRS ADVANCES}, author={Beeler, Benjamin and Zhang, Yongfeng and Hasan, A. T. M. Jahid and Park, Gyuchul and Hu, Shenyang and Mei, Zhi-Gang}, year={2022}, month={Oct} } @article{park_beeler_okuniewski_2023, title={Computational determination of a primary diffusion mode in gamma U-10Mo under irradiation}, volume={574}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2022.154137}, abstractNote={Low enriched uranium ( < 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high-performance research and test reactors and is in the process of being qualified to replace highly-enriched uranium ( ≥ 20 % 235 U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in γ U-10wt. % Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 × 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime where radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. The total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Park, Gyuchul and Beeler, Benjamin and Okuniewski, Maria A.}, year={2023}, month={Feb} } @article{aly_beeler_avramova_2022, title={Investigation of ?-(U, Zr) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics}, volume={567}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2022.153835}, DOI={10.1016/j.jnucmat.2022.153835}, abstractNote={In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ−(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at.% and 100 at.%. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at.% which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ−(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ−(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. This work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.}, journal={JOURNAL OF NUCLEAR MATERIALS}, publisher={Elsevier BV}, author={Aly, Ahmed and Beeler, Benjamin and Avramova, Maria}, year={2022}, month={Aug} } @article{jokisaari_mahbuba_wang_beeler_2022, title={The impact of anisotropic thermal expansion on the isothermal annealing of polycrystalline alpha-uranium}, volume={205}, ISSN={["1879-0801"]}, url={https://doi.org/10.1016/j.commatsci.2022.111217}, DOI={10.1016/j.commatsci.2022.111217}, abstractNote={Although grain growth impacts microstructural evolution in a wide variety of materials systems, the effect of anisotropic thermal expansion on grain boundary mobility and texture evolution has not been widely studied. Anisotropic thermal expansion occurs in multiple non-cubic metals, and the thermomechanical processing behavior of these materials can be better understood with further study into the impact of thermal expansion on grain boundary mobility and texture evolution. In this work, we develop a mesoscale phase field model of grain growth that includes the effect of anisotropic thermal expansion, which is applied to study polycrystalline α-uranium, a highly anisotropic metal. Three-dimensional simulations on polycrystalline α-uranium with and without thermal expansion eigenstrains are performed to study the grain boundary mobility and texture evolution as a function of temperature. A strain-free temperature of 933 K is selected, and the system is studied within the range of 873–933 K at intervals of ten degrees, resulting in increasing thermal eigenstrain with decreasing temperature. We also estimate a grain boundary mobility prefactor and activation energy based on existing experimental data of isothermal annealing of α-uranium. The grain boundary mobility is found to display significant deviation from Arrhenius behavior with the inclusion of thermal expansion eigenstrain as the amount of thermal eigenstrain (and thus elastic strain energy within the system) increases. This result explains an experimentally observed grain boundary mobility deviation from Arrhenius behavior. Furthermore, the texture evolution is affected, such that the grain orientations become less random with increasing thermal eigenstrain, which could explain experimentally observed texture behavior. These results indicate that the effect of thermal expansion should be considered when predicting the thermomechanical processing behavior of α-uranium and other materials with anisotropic thermal expansion.}, journal={COMPUTATIONAL MATERIALS SCIENCE}, publisher={Elsevier BV}, author={Jokisaari, Andrea M. and Mahbuba, Khadija and Wang, Yuhao and Beeler, Benjamin}, year={2022}, month={Apr} } @article{duemmler_lin_woods_karlsson_gakhar_beeler_2022, title={

Evaluation of thermophysical properties of the LiCl-KCl system via ab initio and experimental methods

}, volume={559}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2021.153414}, abstractNote={Molten Salt Reactors (MSRs) are envisioned as a potential pathway to safer, more economical nuclear electricity generation and supply of industrial heat. MSRs under consideration today are either solid-fueled salt-cooled designs or liquid-salt-fueled designs with chloride or fluoride based salts. A significant knowledge gap exists in the data for the fundamental properties relevant to fuels and coolants for MSRs that needs to be addressed in order to expedite the technical readiness level of the MSR design concepts. With the rapid development and improvement of computational materials science, computational methods such as Density Functional Theory (DFT) calculations and ab initio Molecular Dynamics (AIMD) simulations are widely used as an effective and reliable tool to investigate the atomic interaction in materials. In this article, the density of the LiCl-KCl system was determined via AIMD calculations and verified using new experimental analyses. AIMD was further utilized to calculate the compressibility, heat capacity, enthalpy of mixing, and Gibbs free energy of mixing. This work spans a wider range of compositions and temperatures than have previously been explored computationally for this pseudo-binary system and provides the basis for further advanced thermophysical property evaluation utilizing AIMD methods.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Duemmler, Kai and Lin, Yuxiao and Woods, Michael and Karlsson, Toni and Gakhar, Ruchi and Beeler, Benjamin}, year={2022}, month={Feb} } @article{park_beeler_okuniewski_2021, title={An atomistic study of defect energetics and diffusion with respect to composition and temperature in gamma U and gamma U-Mo alloys}, volume={552}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2021.152970}, DOI={10.1016/j.jnucmat.2021.152970}, abstractNote={Uranium-molybdenum (U-Mo) alloys are promising candidates for high-performance research and test reactors, as well as fast reactors. The metastable γ phase, which shows acceptable irradiation performance, is retained by alloying U with Mo with specific quenching conditions. Point defects contribute to the atomic diffusion process, defect clustering, creep, irradiation hardening, and swelling of nuclear fuels, all of which play a role in fuel performance. In this work, properties of point defects in γU and γU-xMo (x = 7, 10, 12 wt.%) were investigated. Vacancy and self-interstitial formation energies in γU and γU-xMo were calculated with molecular dynamics (MD) simulations using an embedded atom method interatomic potential for the U-Mo system. Formation energies of point defects were calculated in the temperature range between 400 K and 1200 K. The vacancy formation energy was higher than the self-interstitial formation energy in both γU and γU-xMo in the evaluated temperature range, which supports the previous results obtained via first-principles calculations and MD simulations. In γU-xMo, the vacancy formation energy decreased with increasing Mo content, whereas the self-interstitial formation energy increased with increasing Mo content in the temperature range of 400 K to 1200 K. The self-diffusion and interdiffusion coefficients were also determined in γU-xMo as a function of temperature. Diffusion of U and Mo atoms in γU-xMo were negligible below 800 K. The self-diffusion and interdiffusion coefficients decreased with increasing Mo concentration, which qualitatively agreed with the previous experimental observations. Point defect formation energies, self-diffusion coefficients, and interdiffusion coefficients in γU-xMo calculated in the present work can be used as input parameters in mesoscale and engineering scale fuel performance modeling.}, journal={JOURNAL OF NUCLEAR MATERIALS}, publisher={Elsevier BV}, author={Park, Gyuchul and Beeler, Benjamin and Okuniewski, Maria A.}, year={2021}, month={Aug} } @article{andersson_matthews_zhang_beeler_2021, title={Density functional theory calculations of the thermodynamic and kinetic properties of point defects in beta-U}, volume={557}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2021.153238}, abstractNote={Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2×2×2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. The diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Andersson, D. A. and Matthews, C. and Zhang, Y. and Beeler, B.}, year={2021}, month={Dec} } @article{beeler_mahbuba_wang_jokisaari_2021, title={Determination of Thermal Expansion, Defect Formation Energy, and Defect-Induced Strain of alpha-U Via ab Initio Molecular Dynamics}, volume={8}, ISSN={["2296-8016"]}, DOI={10.3389/fmats.2021.661387}, abstractNote={Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize its high-temperature body-centered cubic phase for use in nuclear reactors. However, in all metallic fuel forms, the α phase of U remains in some fraction. This phase decomposition due to temperature or compositional variance can play an outsized role on fuel performance and microstructural evolution. Relatively little is known about fundamental point defect properties in α-U at non-zero temperatures, from either computational or experimental studies. This work performs the first thorough evaluation of the α phase of U via ab initio molecular dynamics (AIMD). A number of thermophysical properties are calculated as a function of temperature, including equilibrium lattice parameters, thermal expansion, and heat capacity. These results indicate a two-region behavior, with the transition at 400 K. The thermal expansion/contraction in the a/b direction occurs rapidly from 100 up to 400 K, after which a more linear and gradual change in the lattice constant takes place. The volumetric expansion matches experiments quantitatively, but the individual lattice constant expansion only matches experiments qualitatively. Point defect formation energies and induced lattice strains are also determined as a function of temperature, providing insight on defect populations and the fundamentals of irradiation growth in α-U. Interstitials induce significantly more strain than vacancies, and the nature of that strain is highly dependent on the individual lattice directions. The direction of point defect-induced lattice strain is contrary to the irradiation growth behavior of α-U. This work shows that isolated point defects cannot be the primary driving force responsible for the significant irradiation-induced growth of α-U observed experimentally.}, journal={FRONTIERS IN MATERIALS}, author={Beeler, Benjamin and Mahbuba, Khadija and Wang, Yuhao and Jokisaari, Andrea}, year={2021}, month={Jun} } @article{mahbuba_beeler_jokisaari_2021, title={Evaluation of the anisotropic grain boundaries and surfaces of alpha-U via molecular dynamics}, volume={554}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2021.153072}, DOI={10.1016/j.jnucmat.2021.153072}, abstractNote={Alloys based on uranium-zirconium are gaining renewed interest as fuels for the Versatile Test Reactor and a number of microreactor designs. Implementing metallic fuel in reactors creates the need for robust descriptive and predictive fuel performance modeling. The current state of metallic fuel performance modeling relies on empirical equations derived from historical experiments, which may be unreliable when applied outside of their temperature, power, and composition phase space. One area where such data is lacking is the irradiation behavior of α-U, specifically tearing and porosity formation at the early stages of irradiation. While grain boundaries likely play a key role in this fuel behavior, relatively little is known about grain boundaries in α-U. Thus, we evaluate the grain boundary, surface energy, and work of adhesion of α-U utilizing molecular dynamics. Symmetric tilt grain boundaries (STGBs) are analyzed with the tilt plane oriented along each major crystallographic axis, for a total of eighty unique grain boundaries. The effect of temperature, tilt plane, and misorientation angle on interfacial energies are analyzed. The interfacial energies typically increase with temperature and there is significant variance as a function of misorientation angle, irrespective of the tilt plane. At 500 K, the average surface energy (1.23 J/m2) is approximately 1.5 times the grain boundary energy (0.79 J/m2), and the work of adhesion is approximately twice the grain boundary energy (1.68 J/m2). Orientations for the likely formation of twins and likely failure planes are identified.}, journal={JOURNAL OF NUCLEAR MATERIALS}, publisher={Elsevier BV}, author={Mahbuba, Khadija and Beeler, Benjamin and Jokisaari, Andrea}, year={2021}, month={Oct} } @article{hu_beeler_2021, title={Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity}, volume={8}, ISSN={["2296-8016"]}, DOI={10.3389/fmats.2021.682667}, abstractNote={In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.}, journal={FRONTIERS IN MATERIALS}, author={Hu, Shenyang and Beeler, Benjamin}, year={2021}, month={Jun} } @article{gamble_pastore_cooper_andersson_matthews_beeler_aagesen_barani_pizzocri_2021, title={Improvement of the BISON U3Si2 modeling capabilities based on multiscale developments to modeling fission gas behavior}, volume={555}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2021.153097}, abstractNote={Uranium silicide (U3Si2) is a concept explored as a potential alternative to UO2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U3Si2 has been previously used in research reactors in the form of dispersion fuel which operates at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U3Si2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U3Si2 in the BISON fuel performance code. To support the implementation of a U3Si2 fission gas model in BISON, cluster dynamics simulations of irradiation enhanced Xe diffusion have been carried out. Similarly, MD simulations were used to predict the athermal contribution due to atomic mixing during ballistic damage cascades. By combining our results with literature DFT data for thermal equilibrium diffusion, Xe diffusivity has been described over a wide range of temperatures for in-reactor conditions. These lower length scale informed models are then utilized in the assessment of BISON U3Si2 modeling capabilities by simulating the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the assessment process to identify where further experiments and lower length scale modeling would be beneficial. The multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Gamble, K. A. and Pastore, G. and Cooper, M. W. D. and Andersson, D. A. and Matthews, C. and Beeler, B. and Aagesen, L. K. and Barani, T. and Pizzocri, D.}, year={2021}, month={Nov} } @article{cooper_gamble_capolungo_matthews_andersson_beeler_stanek_metzger_2021, title={Irradiation-enhanced diffusion and diffusion-limited creep in U3Si2}, volume={555}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2021.153129}, abstractNote={U3Si2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U3Si2. Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U3Si2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. A demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U3Si2 in-reactor, a consequence of its high thermal conductivity.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Cooper, M. W. D. and Gamble, K. A. and Capolungo, L. and Matthews, C. and Andersson, D. A. and Beeler, B. and Stanek, C. R. and Metzger, K.}, year={2021}, month={Nov} } @article{beeler_andersson_jiang_zhang_2021, title={Ab initio molecular dynamics investigation of point defects in gamma-U}, volume={545}, ISSN={["1873-4820"]}, url={https://doi.org/10.1016/j.jnucmat.2020.152714}, DOI={10.1016/j.jnucmat.2020.152714}, abstractNote={Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize the high-temperature body-centered cubic γ phase of uranium for use in nuclear reactors. However, relatively little experimental or computational investigation has centered on γ-U, largely due to the mechanical instability of this phase at room temperature. This is particularly problematic for density functional theory calculations that typically investigate 0 K properties. However, ab initio molecular dynamics (AIMD) allows for quantum mechanical-based calculations to be performed at non-zero temperatures. In this work, AIMD simulations are performed to calculate the equilibrium volume for the γ phase of U from 900 K to 1400 K. Utilizing the volume at each temperature, the bulk modulus, the radial distribution function, the interstitial and vacancy formation energies, and the diffusion coefficients are determined.}, journal={JOURNAL OF NUCLEAR MATERIALS}, publisher={Elsevier BV}, author={Beeler, Benjamin and Andersson, David and Jiang, Chao and Zhang, Yongfeng}, year={2021}, month={Mar} } @article{beeler_casagranda_aagesen_zhang_novascone_2020, title={Atomistic calculations of the surface energy as a function of composition and temperature in gamma U-Zr to inform fuel performance modeling}, volume={540}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2020.152271}, abstractNote={Uranium-zirconium alloy fuels are candidates for advanced sodium cooled fast reactors due to their high uranium density, high thermal conductivity, inherent safety and ability to incorporate minor actinides into the fuel. Unlike traditional ceramic UO2 fuel, U–Zr alloys swell rapidly and substantially, but the actual mechanistic process of swelling, including the rate of swelling, is not well understood. Fuel performance models are being developed to describe the swelling process, but these models currently lack the requisite underlying physics and fundamental property data to be truly predictive. In this work, molecular dynamics simulations are utilized to investigate a number of bulk thermophysical properties in γU-Zr, the void surface energy as a function of temperature and composition, and the void free energy. Finally, the effect of surface energy on fuel swelling behavior is demonstrated via finite element based fuel performance simulations, emphasizing the importance of the inclusion of accurate fundamental material properties.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Beeler, Benjamin and Casagranda, Albert and Aagesen, Larry and Zhang, Yongfeng and Novascone, Stephen}, year={2020}, month={Nov} } @article{hu_setyawan_beeler_gan_burkes_2020, title={Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – A cluster dynamics and phase field model}, volume={542}, url={https://doi.org/10.1016/j.jnucmat.2020.152441}, DOI={10.1016/j.jnucmat.2020.152441}, abstractNote={Irradiation examination shows that gas bubble swelling kinetics are much faster after irradiation-induced recrystallization than that prior to recrystallization in U-10 wt% Mo alloy (UMo) fuels. These kinetics imply that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. For the first time, researchers developed a phase-field model of gas bubble evolution integrating microstructure-dependent cluster dynamics to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission, and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of nonequilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission, and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission, and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms that results in fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second-phase precipitates often observed in irradiated materials.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Hu, Shenyang and Setyawan, Wahyu and Beeler, Benjamin W. and Gan, Jian and Burkes, Douglas E}, year={2020}, month={Dec}, pages={152441} } @article{cheniour_tonks_gong_yao_he_harp_beeler_zhang_lian_2020, title={Development of a grain growth model for U3Si2 using experimental data, phase field simulation and molecular dynamics}, volume={532}, ISSN={0022-3115}, url={http://dx.doi.org/10.1016/j.jnucmat.2020.152069}, DOI={10.1016/j.jnucmat.2020.152069}, abstractNote={The purpose of this work is to develop a model for normal grain growth in U3Si2. The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. The mobility was determined by comparing phase field grain growth simulations to the experimental data. From these various methods, we found that the average grain size D in U3Si2 can be estimated over time t using the equation D2−D02=2αMγt, where D0 is the initial average grain size, the geometry factor α=0.96, the average grain boundary mobility M=6.30×10−18e−0.33[eV]kbTm4/(Js) with the Boltzmann constant kb and temperature T, and the average grain boundary energy has been found as a function of temperature, e.g. γ¯=0.83 J/m2 at 673 K.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Cheniour, Amani and Tonks, Michael R. and Gong, Bowen and Yao, Tiankai and He, Lingfeng and Harp, Jason M. and Beeler, Benjamin and Zhang, Yongfeng and Lian, Jie}, year={2020}, month={Apr}, pages={152069} } @article{mei_ye_yacout_beeler_gao_2020, title={First-principles study of the surface properties of uranium carbides}, volume={542}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2020.152257}, abstractNote={Uranium carbides have attracted renewed interest as advanced nuclear fuels for Generation IV reactors. As an important property required for gas bubble modeling in nuclear fuels, the surface energy of uranium carbides is scarce in literature. In this work, we study the surface properties of uranium carbides by first-principles density functional theory calculations. Surface orientations with maximum Miller index up to 3, 2 and 2 are investigated for UC, U2C3 and α-UC2, respectively. By studying the effects of surface termination and chemical potential on surface energy, we identify the factors that determines the surface stability. From the calculated surface energies, the surface properties of uranium carbide single crystals are obtained from Wulff construction, including equilibrium morphology, dominant surface orientation and area weighted surface energy.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Mei, Zhi-Gang and Ye, Bei and Yacout, Abdellatif M. and Beeler, Benjamin and Gao, Yipeng}, year={2020}, month={Dec} } @article{aagesen_andersson_beeler_cooper_gamble_miao_pastore_tonks_2020, title={Phase-field simulations of intergranular fission gas bubble behavior in U3Si2 nuclear fuel}, volume={541}, url={https://doi.org/10.1016/j.jnucmat.2020.152415}, DOI={10.1016/j.jnucmat.2020.152415}, abstractNote={U3Si2 is a potential accident-tolerant fuel that shows promise due to its high thermal conductivity and higher uranium density relative to UO2. However, its swelling and fission gas release behavior in light water reactor (LWR) conditions is relatively unknown. To provide mechanistic insight and determine parameters for engineering-scale fuel performance modeling of pellet-form U3Si2, phase-field simulations of the growth, interconnection, and venting of intergranular fission gas bubbles were performed. The fractional coverage of the grain boundary and the fraction of bubble area that is vented were calculated as a function of time. From the simulation data, the fractional grain boundary coverage at saturation, an important parameter needed in engineering-scale modeling of swelling and fission gas release, was determined. Multiple simulations were run to determine the uncertainty in the calculated value. The effect of model assumptions and input parameters that are not well known was evaluated. Simulation results are compared to related theoretical and computational work. Based on the simulation results, a value of 0.60 for the fractional grain boundary coverage at saturation is recommended for U3Si2 fuel.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Aagesen, Larry K. and Andersson, David and Beeler, Benjamin W. and Cooper, Michael W.D. and Gamble, Kyle A. and Miao, Yinbin and Pastore, Giovanni and Tonks, Michael R.}, year={2020}, month={Dec}, pages={152415} } @article{beeler_cooper_mei_schwen_zhang_2021, title={Radiation driven diffusion in γU-Mo}, volume={543}, url={https://doi.org/10.1016/j.jnucmat.2020.152568}, DOI={10.1016/j.jnucmat.2020.152568}, abstractNote={A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the U. S. High-Performance Research Reactors. A critical phenomenon of interest regarding U-Mo monolithic fuel is the large amount of swelling that takes place during operation, particularly at high fission densities. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic and kinetic properties into mesoscale and engineering fuel performance modeling codes. One such property where there exists incomplete data is the diffusion of relevant species under irradiation. Fuel performance swelling predictions rely on an accurate representation of diffusion in order to determine the rate of fission gas swelling and the local microstructural evolution. In this work, we present molecular dynamics simulations of the radiation driven diffusion of U, Mo, and Xe in U-Mo nuclear fuels. Diffusion coefficients for each species are determined over a range of temperatures and compositions. Updated diffusion coefficients are presented that are applicable under irradiation and incorporate both intrinsic and radiation driven diffusion.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Beeler, Benjamin and Cooper, Michael W.D. and Mei, Zhi-Gang and Schwen, Daniel and Zhang, Yongfeng}, year={2021}, month={Jan}, pages={152568} } @article{beeler_hu_zhang_gao_2020, title={A improved equation of state for Xe gas bubbles in gamma U-Mo fuels}, volume={530}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2019.151961}, abstractNote={A monolithic fuel design based on a U–Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U–Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in γU-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Beeler, Benjamin and Hu, Shenyang and Zhang, Yongfeng and Gao, Yipeng}, year={2020}, month={Mar} } @article{dacus_beeler_schwen_2019, title={Calculation of threshold displacement energies in UO2}, volume={520}, ISSN={["1873-4820"]}, DOI={10.1016/j.jnucmat.2019.04.002}, abstractNote={Despite the extensive utilization of uranium dioxide (UO2) as a fuel in commercial nuclear reactors, there is only minimal information regarding the fundamental nature of radiation damage at high temperatures, such as those experienced by the fuel under operation. In this work, molecular dynamics simulations have been performed to determine the threshold displacement energy (Ed) for oxygen and uranium in UO2 at 1500 K. Three definitions of displacement energy were employed to fully study the nature of low energy radiation damage: 1) the probability of having the primary knock-on atom (PKA) leave its original lattice site, 2) the probability that the PKA will permanently displace atoms from their original lattice site, and 3) the probability of forming a stable Frenkel pair. Additionally, four unique interatomic potentials were utilized to investigate uncertainties associated with potential choice in high temperature radiation damage studies in UO2. This work provides critical insight into the high temperature behavior of radiation damage in UO2, as well as the variation in behavior between oxygen and uranium PKAs.}, journal={JOURNAL OF NUCLEAR MATERIALS}, author={Dacus, Benjamin and Beeler, Benjamin and Schwen, Daniel}, year={2019}, month={Jul}, pages={152–164} } @article{beeler_baskes_andersson_cooper_zhang_2017, title={A modified Embedded-Atom Method interatomic potential for uranium-silicide}, volume={495}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-85027967820&partnerID=MN8TOARS}, DOI={10.1016/j.jnucmat.2017.08.025}, abstractNote={Uranium-silicide (U-Si) fuels are being pursued as a possible accident tolerant fuel (ATF). This uranium alloy fuel benefits from higher thermal conductivity and higher fissile density compared to uranium dioxide (UO2). In order to perform engineering scale nuclear fuel performance simulations, the material properties of the fuel must be known. Currently, the experimental data available for U-Si fuels is rather limited. Thus, multiscale modeling efforts are underway to address this gap in knowledge. In this study, a semi-empirical modified Embedded-Atom Method (MEAM) potential is presented for the description of the U-Si system. The potential is fitted to the formation energy, defect energies and structural properties of U3Si2. The primary phase of interest (U3Si2) is accurately described over a wide temperature range and displays good behavior under irradiation and with free surfaces. The potential can also describe a variety of U-Si phases across the composition spectrum.}, journal={Journal of Nuclear Materials}, author={Beeler, B. and Baskes, M. and Andersson, D. and Cooper, M.W.D. and Zhang, Y.}, year={2017}, pages={267–276} } @inproceedings{beeler_baskes_andersson_zhang_2017, title={A modified embedded-atom method interatomic potential for U-Si}, volume={116}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-85033481319&partnerID=MN8TOARS}, booktitle={Transactions of the American Nuclear Society}, author={Beeler, B. and Baskes, M. and Andersson, D. and Zhang, Y.}, year={2017}, pages={407–409} } @book{aagesen_ahmed_beeler_schwen_zhang_andersson_2017, title={Multi-Scale Modeling of Swelling in Accident-Tolerant U3Si2 Fuel}, url={https://doi.org/10.2172/1472119}, DOI={10.2172/1472119}, abstractNote={U3Si2 is a promising candidate for use as an accident-tolerant fuel for light water reactors. A multi-scale computational approach was used to calculate swelling in pellet-form U3Si2 fuel. Swelling was assumed to be equal to the volume fraction of fission gas bubbles in the fuel, and the evolution of bubble volume fraction was determined from phase-field simulations. To parameterize the phase-field model, density-functional theory and molecular dynamics simulations were performed. To enable molecular dynamics simulations, a new interatomic potential for the U-Si system was developed based on the modified embedded atom method. A new phase-field model based on a grand-potential functional was also developed. The model is applicable to both intergranular and intragranular bubbles. To calculate the volume fraction of bubbles, the microstructure was decomposed into regions consisting of only intragranular or intergranular bubbles based on a truncated octahedral grain structure, and growth of the bubbles in the two regions was simulated separately. The total swelling was then calculated based on a a weighted average of bubble volume fraction in the two regions. Total swelling was of the same order of magnitude, but larger than that predicted by the existing empirical swelling model used in BISON and a rate-theory based model. Limitations of the present approach and suggestions for improvement are presented.}, institution={Office of Scientific and Technical Information (OSTI)}, author={Aagesen, Larry and Ahmed, Karim and Beeler, Benjamin and Schwen, Daniel and Zhang, Yongfeng and Andersson, David}, year={2017}, month={Sep} } @book{zhang_beeler_aagesen_jiang_ahmed_yu_schwen_andersson_baskes_cooper_et al._2017, title={Progress update on lower length scale research and development on U3Si2 fuel and FeCrAl cladding}, url={https://doi.org/10.2172/1472101}, DOI={10.2172/1472101}, abstractNote={This report summarizes the lower length scale model development and studies in regards to two accident tolerant (ATF) concepts, U3Si2 and FeCrAl. U3Si2 is proposed to replace UO2 fuel for its high thermal conductivity and higher uranium density. FeCrAl alloys are proposed to replace Zr based cladding for their high mechanical strength and good corrosion resistance. Before applying them in reactors, sufficient data are needed to establish their predictable in-pile performance, but largely missing, and they are costly and time consuming to obtain by experimental approaches. Under the NEAMS ATF high-impact-problem (HIP), multiscale modeling and simulations are used to develop atomistic and mesoscale tools, obtain fundamental material properties, and assess the behaviors of these two ATF concepts. The accomplishments in FY17 are briefly introduced in this report including the development of interatomic potentials, assessment of swelling in U3Si2, and Cr precipitation in FeCrAl alloys.}, institution={Office of Scientific and Technical Information (OSTI)}, author={Zhang, Yongfeng and Beeler, Benjamin and Aagesen, Larry and Jiang, Chao and Ahmed, K. and Yu, J. and Schwen, D. and Andersson, D. and Baskes, M. and Cooper, M. and et al.}, year={2017}, month={Sep} } @article{beeler_asta_hosemann_gr?nbech-jensen_2016, title={Effect of strain and temperature on the threshold displacement energy in body-centered cubic iron}, volume={474}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84963585213&partnerID=MN8TOARS}, DOI={10.1016/j.jnucmat.2016.03.017}, abstractNote={The threshold displacement energy (TDE) is the minimum amount of kinetic energy required to displace an atom from its lattice site. The magnitude of the TDE displays significant variance as a function of the crystallographic direction, system temperature and applied strain, among a variety of other factors. It is critically important to determine an accurate value of the TDE in order to calculate the total number of displacements due to a given irradiation condition, and thus to understand the materials response to irradiation. In this study, molecular dynamics simulations have been performed to calculate the threshold displacement energy in body-centered cubic iron as a function of strain and temperature. With applied strain, a decrease of the TDE of up to approximately 14 eV was observed. A temperature increase from 300 K to 500 K can result in an increase of the TDE of up to approximately 9 eV.}, journal={Journal of Nuclear Materials}, author={Beeler, B. and Asta, M. and Hosemann, P. and Gr?nbech-Jensen, N.}, year={2016}, pages={113–119} } @article{moore_beeler_deo_baskes_okuniewski_2015, title={Atomistic modeling of high temperature uranium-zirconium alloy structure and thermodynamics}, volume={467}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84946780407&partnerID=MN8TOARS}, DOI={10.1016/j.jnucmat.2015.10.016}, abstractNote={A semi-empirical Modified Embedded Atom Method (MEAM) potential is developed for application to the high temperature body-centered-cubic uranium–zirconium alloy (γ-U–Zr) phase and employed with molecular dynamics (MD) simulations to investigate the high temperature thermo-physical properties of U–Zr alloys. Uranium-rich U–Zr alloys (e.g. U–10Zr) have been tested and qualified for use as metallic nuclear fuel in U.S. fast reactors such as the Integral Fast Reactor and the Experimental Breeder Reactors, and are a common sub-system of ternary metallic alloys like U–Pu–Zr and U–Zr–Nb. The potential was constructed to ensure that basic properties (e.g., elastic constants, bulk modulus, and formation energies) were in agreement with first principles calculations and experimental results. After which, slight adjustments were made to the potential to fit the known thermal properties and thermodynamics of the system. The potentials successfully reproduce the experimental melting point, enthalpy of fusion, volume change upon melting, thermal expansion, and the heat capacity of pure U and Zr. Simulations of the U–Zr system are found to be in good agreement with experimental thermal expansion values, Vegard's law for the lattice constants, and the experimental enthalpy of mixing. This is the first simulation to reproduce the experimental thermodynamics of the high temperature γ-U–Zr metallic alloy system. The MEAM potential is then used to explore thermodynamics properties of the high temperature U–Zr system including the constant volume heat capacity, isothermal compressibility, adiabatic index, and the Grüneisen parameters.}, journal={Journal of Nuclear Materials}, author={Moore, A.P. and Beeler, B. and Deo, C. and Baskes, M.I. and Okuniewski, M.A.}, year={2015}, pages={802–819} } @inproceedings{beeler_asta_hosemann_gr?nbech-jensen_2015, title={Effect of interfaces on radiation damage accumulation in FeNiAl maraging steels}, volume={112}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84988891995&partnerID=MN8TOARS}, booktitle={Transactions of the American Nuclear Society}, author={Beeler, B. and Asta, M. and Hosemann, P. and Gr?nbech-Jensen, N.}, year={2015}, pages={328–329} } @article{beeler_asta_hosemann_grønbech-jensen_2015, title={Effects of applied strain on radiation damage generation in body-centered cubic iron}, volume={459}, DOI={10.1016/j.jnucmat.2014.12.111}, abstractNote={Radiation damage in body-centered cubic (BCC) Fe has been extensively studied by computer simulations to quantify effects of temperature, impinging particle energy, and the presence of extrinsic particles. However, limited investigation has been conducted into the effects of mechanical stresses and strain. In a reactor environment, structural materials are often mechanically strained, and an expanded understanding of how this strain affects the generation of defects may be important for predicting microstructural evolution and damage accumulation under such conditions. In this study, we have performed molecular dynamics simulations in which various types of homogeneous strains are applied to BCC Fe and the effect on defect generation is examined. It is found that volume-conserving shear strains yield no statistically significant variations in the stable number of defects created via cascades in BCC Fe. However, strains that result in volume changes are found to produce significant effects on defect generation.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Beeler, Benjamin and Asta, Mark and Hosemann, Peter and Grønbech-Jensen, Niels}, year={2015}, month={Apr}, pages={159–165} } @article{miao_beeler_deo_baskes_okuniewski_stubbins_2015, title={Defect structures induced by high-energy displacement cascades in γ uranium}, volume={456}, DOI={10.1016/j.jnucmat.2014.09.016}, abstractNote={Abstract Displacement cascade simulations were conducted for the γ uranium system based on molecular dynamics. A recently developed modified embedded atom method (MEAM) potential was employed to replicate the atomic interactions while an embedded atom method (EAM) potential was adopted to help characterize the defect structures induced by the displacement cascades. The atomic displacement process was studied by providing primary knock-on atoms (PKAs) with kinetic energies from 1 keV to 50 keV. The influence of the PKA incident direction was examined. The defect structures were analyzed after the systems were fully relaxed. The states of the self-interstitial atoms (SIAs) were categorized into various types of dumbbells, the crowdion, and the octahedral interstitial. The voids were determined to have a polyhedral shape with {1 1 0} facets. The size distribution of the voids was also obtained. The results of this study not only expand the knowledge of the microstructural evolution in irradiated γ uranium, but also provide valuable references for the radiation-induced defects in uranium alloy fuels.}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Miao, Yinbin and Beeler, Benjamin and Deo, Chaitanya and Baskes, Michael I. and Okuniewski, Maria A. and Stubbins, James F.}, year={2015}, month={Jan}, pages={1–6} } @inproceedings{beeler_deo_baskes_okuniewski_2013, title={Atomistic investigations of intrinsic and extrinsic point defects in bcc uranium}, volume={1547 STP}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84875784961&partnerID=MN8TOARS}, DOI={10.1520/STP104141}, abstractNote={Metallic alloys of uranium show great potential as transmutation fuels that could be used to burn long-lived and high-heat-producing minor actinides and fission products in nuclear reactors. In fuels, fission and radiation damage result in the production of a large number of intrinsic point defects and extrinsic fission atoms. Radiation damage and diffusion (processes heavily dependent on point defects), as well as fission product behavior, are important to the understanding of the behavior of these metallic fuel alloys. Of the fission products, fission gases Xe, Kr, and the decay product He are of special importance, as they migrate and form bubbles detrimentally affecting fuel properties. In this work, several systems of body-centered-cubic (gamma) U are examined through a semi-empirical interatomic potential based on the modified embedded-atom method. The vacancy formation energy is analyzed as a function of pressure and is used to determine the stable vacancy formation energy at ambient conditions. The vacancy formation energy as a function of temperature is analyzed for high-temperature systems. Interatomic potentials are developed and implemented in the investigation of He, Xe, and Kr point defects in the gamma phase of uranium. For all fission gases studied, the most energetically favorable location is the substitutional position, with helium having the lowest formation energies of the species investigated.}, booktitle={ASTM Special Technical Publication}, publisher={ASTM International}, author={Beeler, Benjamin and Deo, Chaitanya and Baskes, Michael and Okuniewski, Maria}, year={2013}, pages={231–247} } @inproceedings{moore_beeler_baskes_okuniewski_deo_2013, title={Atomistic ordering in body centered cubic Uranium-Zirconium alloy}, volume={1514}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84888074338&partnerID=MN8TOARS}, DOI={10.1557/opl.2013.517}, abstractNote={ABSTRACTThe metallic binary-alloy fuel Uranium-Zirconium is important for the use of the new generation of advanced fast reactors. Uranium-Zirconium goes through a phase transition at higher temperatures to a (gamma) Body Centered Cubic (BCC) phase. The BCC high temperature phase is particularly important, since the BCC phase corresponds to the temperature range in which the fast reactors will operate. A semi-empirical MEAM (Modified Embedded Atom Method) potential is presented for Uranium-Zirconium. The physical properties of the Uranium-Zirconium binary alloy were reproduced using Molecular Dynamics (MD) simulations and Monte Carlo (MC) simulations with the MEAM potential. This is a large step in making a computationally acceptable fuel performance code.}, booktitle={Materials Research Society Symposium Proceedings}, author={Moore, A.P. and Beeler, B. and Baskes, M. and Okuniewski, M. and Deo, C.S.}, year={2013}, pages={27–35} } @article{beeler_deo_baskes_okuniewski_2012, title={Atomistic properties of γ uranium}, volume={24}, DOI={10.1088/0953-8984/24/7/075401}, abstractNote={The properties of the body-centered cubic γ phase of uranium (U) are calculated using atomistic simulations. First, a modified embedded-atom method interatomic potential is developed for the high temperature body-centered cubic (γ) phase of U. This phase is stable only at high temperatures and is thus relatively inaccessible to first principles calculations and room temperature experiments. Using this potential, equilibrium volume and elastic constants are calculated at 0 K and found to be in close agreement with previous first principles calculations. Further, the melting point, heat capacity, enthalpy of fusion, thermal expansion and volume change upon melting are calculated and found to be in reasonable agreement with experiment. The low temperature mechanical instability of γ U is correctly predicted and investigated as a function of pressure. The mechanical instability is suppressed at pressures greater than 17.2 GPa. The vacancy formation energy is analyzed as a function of pressure and shows a linear trend, allowing for the calculation of the extrapolated zero pressure vacancy formation energy. Finally, the self-defect formation energy is analyzed as a function of temperature. This is the first atomistic calculation of γ U properties above 0 K with interatomic potentials.}, number={7}, journal={J. Phys.: Condens. Matter}, publisher={IOP Publishing}, author={Beeler, Benjamin and Deo, Chaitanya and Baskes, Michael and Okuniewski, Maria}, year={2012}, month={Feb}, pages={075401} } @inproceedings{beeler_deo_baskes_okuniewski_2012, title={Calculation of the displacement energy in b.c.c. U at 800 K}, volume={106}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84876400924&partnerID=MN8TOARS}, booktitle={Transactions of the American Nuclear Society}, author={Beeler, B. and Deo, C. and Baskes, M. and Okuniewski, M.}, year={2012}, pages={1214–1215} } @article{beeler_deo_baskes_okuniewski_2013, title={First principles calculations of the structure and elastic constants of α, β and γ uranium}, volume={433}, DOI={10.1016/j.jnucmat.2012.09.019}, abstractNote={This study analyzes structural and elastic properties of five uranium crystal structures: the face centered orthorhombic A20 (α phase), the tetragonal D8b (β phase), body centered tetragonal (bct), body centered cubic (γ phase) and face centered cubic structures. Calculations are performed within the density functional theory framework employing the Projector Augmented Wave method and the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE–GGA) of the exchange correlation. The elastic constants are used to compute polycrystalline elastic moduli, Poisson’s ratio and the Debye temperature for all five structures. The α and γ phase properties are compared with theoretical and experimental results. The complex tetragonal 30 atom beta phase is examined in detail. Representation of the β phase by a bct structure is examined; we find that the structure of the β phase is significantly different from the bct phase but exhibits similar elastic properties. This is the first comprehensive investigation into the elastic constants of uranium utilizing the PBE–GGA.}, number={1-3}, journal={Journal of Nuclear Materials}, publisher={Elsevier BV}, author={Beeler, Benjamin and Deo, Chaitanya and Baskes, Michael and Okuniewski, Maria}, year={2013}, month={Feb}, pages={143–151} } @article{hayward_beeler_deo_2012, title={Multiple hydrogen trapping at monovacancies}, volume={92}, DOI={10.1080/09500839.2012.657702}, abstractNote={Novel structures for multiple hydrogen atoms trapped at a monovacancy are discussed. Using atomistic simulations based on semiempirical interatomic potentials and density functional theory, we find low-energy configurations for four, five, and six hydrogen atoms around a monovacancy different than those that have been previously studied in the literature. The energetics of hydrogen binding are compared to results, both theoretical and experimental, previously published in the literature. We argue that up to four hydrogen atoms may be exothermically bound to monovacancy.}, number={5}, journal={Philosophical Magazine Letters}, publisher={Informa UK Limited}, author={Hayward, Erin and Beeler, Benjamin and Deo, Chaitanya}, year={2012}, month={May}, pages={217–225} } @article{beeler_good_rashkeev_deo_baskes_okuniewski_2012, title={First-principles calculations of the stability and incorporation of helium, xenon and krypton in uranium}, volume={425}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-84860435889&partnerID=MN8TOARS}, DOI={10.1016/j.jnucmat.2011.08.014}, abstractNote={While metallic fuels have a long history of reactor use, their fundamental physical and thermodynamic properties are not well understood. Many metallic nuclear fuels are body-centered cubic alloys of uranium that swell under fission conditions, creating fission product gases such as helium, xenon and krypton. In this paper, helium, xenon, and krypton point defects are investigated in the α and γ phases of metallic uranium using first principles calculations. A density functional theory (DFT) framework is utilized with projector augmented-wave (PAW) pseudopotentials. Formation and incorporation energies of He, Xe, and Kr are calculated at various defect positions for the prediction of fission gas behavior in uranium. In most cases, defect energies follow a size effect, with helium incorporation and formation energies being the smallest. The most likely position for the larger Xe and Kr atoms in uranium is the substitutional site. Helium atoms are likely to be found in a wide variety of defect positions due to the comparable formation energies of all defect configurations analyzed. This is the first detailed study of the stability and incorporation of fission gases in uranium.}, number={1-3}, journal={Journal of Nuclear Materials}, author={Beeler, B. and Good, B. and Rashkeev, S. and Deo, C. and Baskes, M. and Okuniewski, M.}, year={2012}, pages={2–7} } @article{beeler_good_rashkeev_deo_baskes_okuniewski_2010, title={First principles calculations for defects in U}, volume={22}, url={http://www.scopus.com/inward/record.url?eid=2-s2.0-78649825146&partnerID=MN8TOARS}, DOI={10.1088/0953-8984/22/50/505703}, abstractNote={Uranium (U) exhibits a high temperature body-centered cubic (bcc) allotrope that is often stabilized by alloying with transition metals such as Zr, Mo, and Nb for technological applications. One such application involves U–Zr as nuclear fuel, where radiation damage and diffusion (processes heavily dependent on point defects) are of vital importance. Several systems of U are examined within a density functional theory framework utilizing projector augmented wave pseudopotentials. Two separate generalized gradient approximations of the exchange-correlation are used to calculate defect properties and are compared. The bulk modulus, the lattice constant, and the Birch–Murnaghan equation of state for the defect free bcc uranium allotrope are calculated. Defect parameters calculated include energies of formation of vacancies in the α and γ allotropes, as well as self-interstitials, Zr interstitials, and Zr substitutional defects for the γ allotrope. The results for vacancies agree very well with experimental and previous computational studies. The most probable self-interstitial site in γ-U is the ⟨110⟩ dumbbell, and the most probable defect location for dilute Zr in γ-U is the substitutional site. This is the first detailed study of self-defects in the bcc allotrope of U and also the first comprehensive study of dilute Zr defects in γ-U.}, number={50}, journal={Journal of Physics Condensed Matter}, author={Beeler, B. and Good, B. and Rashkeev, S. and Deo, C. and Baskes, M. and Okuniewski, M.}, year={2010} }