Atomistic simulations of temperature and direction dependent threshold displacement energies in alpha- and gamma-uranium

EY Chen and C Deo and R Dingreville, COMPUTATIONAL MATERIALS SCIENCE, 157, 75-86 (2019).

DOI: 10.1016/j.commatsci.2018.10.026

We performed a systematic study of the threshold displacement energy (E-d) in metallic uranium as a function of both the recoil direction and temperature using Molecular Dynamics simulations. We developed a novel orientation sampling scheme that utilizes crystallographic symmetrical geodesic grids to select directions from the orientation fundamental zone to study the directional dependency. Additionally, we studied the temperature dependency by considering both the alpha-uranium phase, corresponding to the ground state for temperatures ranging from 0 K to 600 K, and the gamma-uranium phase, corresponding to high-temperature state for temperatures above 900 K. In this study, we compared several definitions of the threshold energy: a direction-specific threshold displacement energy (E-d (theta, phi)) an angle-averaged threshold energy (E-d(ave)), a production probability threshold displacement energy (E-d(pp)), and a defect count threshold displacement energy (E-d(dc)). The directionspecific threshold displacement energies showed large angular anisotropy and variations in accordance with crystallographic considerations. Specifically, preferred defect channeling directions were observed in the 1 2 0, 1 (2) over bar0, 1 (1) over bar1 directions for the alpha-uranium, and 0 0 1, 1 1 1 directions for the gamma-uranium. The production probability threshold displacement energy (E-d(pp)) is calculated as approximately 99.2659 eV at 10 K (alpha-U), 103.4980 eV at 300 K (alpha-U), 76.0915 eV at 600 K (alpha-U), and 42.9929 eV at 900 K (gamma-U). With exception of those calculated at 10 K, threshold displacement energies decrease with increasing temperature. Analyses of the stable defect structures showed that the most commonly observed interstitial configuration in alpha- uranium consists of a 0 1 0 dumbbell-like interstitial; while in y-uranium no preferential defect configuration could be identified due to thermally-induced lattice instabilities at the elevated temperatures.

Return to Publications page