Computer simulation of structural modifications induced by highly energetic ions in uranium dioxide

The structural modification caused by the high-energy-ion irradiation of single-crystalline uranium dioxide was simulated by the molecular dynamics method. As the initial condition, high kinetic energy was supplied to the individual atoms within a cylindrical region of nanometer-order radius located in the center of the specimen. The potential proposed by Basak et al. [C.B. Basak, A.K. Sengupta, H.S. Kamath, J. Alloys Compd. 360 (2003) 210–216] was utilized to calculate interaction between atoms. The supplied kinetic energy was first spent to change the crystal structure into an amorphous one within a short period of about 0.3 ps, then it dissipated in the specimen. The amorphous track radius Ra was determined as a function of the effective stopping power gSe, i.e., the kinetic energy of atoms per unit length created by ion irradiation (Se: electronic stopping power, g: energy transfer ratio from stopping power to lattice vibration energy). It was found that the relationship between Ra and gSe follows the relation Ra2=aln(gSe)+b. Compared to the case of Si and β-cristobalite single crystals, it was harder to produce amorphous track because of the long range interaction between U atoms