Molecular dynamics simulation of helium–vacancy interaction in plutonium

Abstract

The formation energies of small He n V m clusters ( n and m denote the number of He atoms and vacancy, respectively) in Pu have been calculated with molecular dynamics (MD) simulations using the embedded atom method (EAM) potential, the Mores potential and the Lennard–Jones potential for describing the interactions of Pu–Pu, Pu–He and He–He, respectively. The binding energies of an interstitial He atom, an isolated vacancy and a self-interstitial Pu atom to a He n V m cluster are also obtained from the calculated formation energies of the clusters. All the binding energies mainly depend on the He–vacancy ratio ( n/ m) of clusters rather than the clusters size. With the increase of the n/ m ratio, the binding energies of a He atom and a Pu atom to a He n V m cluster decrease with the ratio, and the binding energy of a vacancy to a He n V m cluster increases. He atoms act as a catalyst for the formation of He n V m clusters.