Milestone report: The simulation of radiation driven gas diffusion in UO 2 at low temperature
Below 1000 K it is thought that fission gas diffusion in nuclear fuel during irradiation occurs through atomic mixing due to radiation damage. This is an important process for nuclear reactor performance as it affects fission gas release, particularly from the periphery of the pellet where such temperatures are normal. Here we present a molecular dynamics study of Xe and Kr diffusion due to irradiation. Thermal spikes and cascades have been used to study the electronic stopping and ballistic phases of damage, respectively. Our results predict that O and Kr exhibit the greatest diffusivity and U the least, while Xe lies in between. It is concluded that the ballistic phase does not sufficiently account for the experimentally observed diffusion. Preliminary thermal spike calculations indicate that the electronic stopping phase generates greater fission gas displacement than the ballistic phase, although further calculation must be carried out to confirm this. A good description of the system by the empirical potentials is important over the very wide temperatures induced during thermal spike and damage cascade simulations. This has motivated the development of a parameter set for gas-actinide and gas-oxygen interactions that is complementary for use with a recent many-body potential set. A comprehensive set of density functional theory (DFT) calculations were used to study Xe and Kr incorporation at a number of sites in CeO 2 , ThO 2 , UO 2 and PuO 2 . These structures were used to fit a potential, which was used to generate molecular dynamics (MD) configurations incorporating Xe and Kr at 300 K, 1500 K, 3000 K and 5000 K. Subsequent matching to the forces predicted by DFT for these MD configurations was used to refine the potential set. This fitting approach ensured weighted fitting to configurations that are thermodynamically significant over a broad temperature range, while avoiding computationally expensive DFT-MD calculations. The resultant gas potentials were validated against DFT binding energies and are suitable for simulating combinations of Xe and Kr in solid solutions of CeO 2 , ThO 2 , UO 2 and PuO 2 , providing a powerful tool for the atomistic simulation of conventional nuclear reactor fuel UO 2 as well as advanced MOX fuels.