A thermal modelling of displacement cascades in uranium dioxide

The space and time dependent temperature distribution was studied in uranium dioxide during displacement cascades simulated by classical molecular dynamics (MD). The energy for each simulated radiation event ranged between 0.2 keV and 20 keV in cells at initial temperatures of 700 K or 1400 K. Spheres into which atomic velocities were rescaled (thermal spikes) have also been simulated by MD to simulate the thermal excitation induced by displacement cascades. Equipartition of energy was shown to occur in displacement cascades, half of the kinetic energy of the primary knock-on atom being converted after a few tenths of picoseconds into potential energy. The kinetic and potential parts of the system energy are however subjected to little variations during dedicated thermal spike simulations. This is probably due to the velocity rescaling process, which impacts a large number of atoms in this case and would drive the system away from a dynamical equilibrium. This result makes questionable MD simulations of thermal spikes carried out up to now (early 2014). The thermal history of cascades was compared to the heat equation solution of a punctual thermal excitation in UO2. The maximum volume brought to a temperature above the melting temperature during the simulated cascade events is well reproduced by this simple model. This volume eventually constitutes a relevant estimate of the volume affected by a displacement cascade in UO2. This definition of the cascade volume could also make sense in other materials, like iron. (C) 2014 Elsevier B.V. All rights reserved.