Molecular Dynamics Simulation of Thermodynamic Properties in Uranium Dioxide

In the present study, we investigated the thermodynamic properties of uranium dioxide (UO₂) by molecular dynamics (MD) simulations. As for solid UO₂, the lattice parameter, density, and enthalpy obtained by MD simulations were in good agreement with existing experimental data and previous theoretical predictions. The calculated thermal conductivities matched the experiment results at the midtemperature range but were underestimated at very low and very high temperatures. The calculation results of mean square displacement represented the stability of uranium at all temperatures and the high mobility of oxygen toward 3000 K. By fitting the diffusivity constant of oxygen with the Vogel-Fulcher-Tamman law, we noticed a secondary phase transition near 2006.4 K, which can be identified as a “strong” to “fragile” supercooled liquid or glass phase transition in UO₂. By fitting the oxygen diffusion constant with the Arrhenius equation, activation energies of 2.0 and 2.7 eV that we obtained were fairly close to the recommended values of 2.3 to 2.6 eV.