A comprehensive investigation on point defects and their clustering behavior in nonstoichiometric uranium dioxide UO2±x is carried out using LSDA+U method based on density functional theory. Accurate energetic information and charge transfers available so far are obtained. With these energies that have improved more than 50% over that of pure GGA and LDA, we show the density functional theory predicts the predominance of oxygen defects over uranium ones at any compositions, which is possible only after treated the localized 5f electrons properly. Calculations also suggest an upper bound of x ∼ 0.03 for oxygen clusters to start off. The volume change induced by point uranium defects is monotonic but nonlinear, whereas for oxygen defects, increase x always reduces the system volume linearly, except dimers that require extra space for accommodation, which has been identified as meta-stable ionic molecule. Though oxygen dimers usually occupy Willis O ′′ sites and mimic a single oxygen in energetics and charge state, they are rare at ambient conditions. Its decomposition process and vibrational properties have been studied carefully. To obtain a general clustering mechanism in anion-excess fluorites systematically, we also analyze the local stabilities of possible basic clustering modes of oxygen defects. The result shows an unified way to understand the structure of Willis type and cuboctahedral clusters in UO2+x and β-U4O9. Finally we generalize the point defect model to the independent clusters approximation to include clustering effects, the impact on defect populations is discussed.

A first-principles calculation for uranium dioxide (UO 2 ) in an antiferromagnetic structure with four types of point defects, uranium vacancy, oxygen vacancy, uranium interstitial, and oxygen interstitial, has been performed by the projector-augmented-wave method with generalized gradient approximation combined with the Hubbard U correction. Defect formation energies are estimated under lattice relaxation for supercells containing 1, 2, and 8 unit cells of UO 2 . The electronic structure, the atomic displacement and the stability of defected systems are obtained, and the effects of cell sizes on these properties are discussed. The results form a self-consistent dataset of formation energies and atomic distance variations of various point defects in UO 2 with relatively high precision. We show that a supercell with 8 UO 2 unit cells or larger is necessary to investigate the defect behavior with reliable precision, since point defects have a wide-ranging effect, not only on the first nearest neighbor atoms of the defect, but on the second neighbors and on more distant atoms.

Structural behavior of UO2 under high-pressure up to 300 GPa has been studied by first-principles calculations with LSDA+U approximation. The results show that a pressure induced structural transition to the cotunnite-type (orthorhombic Pnma) phase occurs at 38 GPa. It agrees well with experiment observed ∼42 GPa. A new iso-structural transition following that is also predicted taking place from 80 to 130 GPa, which has not yet been observed in experiments. Further high compression beyond 226 GPa will result in a metallic and paramagnetic transition. It corresponds to a volume of 90Å 3 per cell, in a good agreement with previous theoretical analysis in the reduction of volume required to delocalize 5f states.

Self-defect clusters in bulk matrix might affect the thermodynamic behavior of fission gases in nuclear fuel such as uranium dioxide. With first-principles local spin-density approximation plus U calculations and taking xenon as a prototype, we find that the influence of oxygen defect clusters on the thermodynamics of gas atoms is prominent, which increases the solution energy of xenon by a magnitude of 0.5 eV, about 43% of the energy difference between the two lowest lying states at 700 K. Calculation also reveals a thermodynamic competition between the uranium vacancy and trivacancy sites to incorporate xenon in hyperstoichiometric regime at high temperatures. The results show that in hypostoichiometric regime neutral trivacancy sites are the most favored position for diluted xenon gas, whereas in hyperstoichiometric condition they prefer to uranium vacancies even after taking oxygen self-defect clusters into account at low temperatures, which not only confirms previous studies but also extends the conclusion to more realistic fuel operating conditions. The observation that gas atoms are ionized to a charge state of Xe + when at a uranium vacancy site due to strong Madelung potential implies that one can control temperature to tune the preferred site of gas atoms and then the bubble growth rate. A solution to the notorious metastable states difficulty that frequently encountered in density functional theory plus U applications, namely, the quasiannealing procedure, is also discussed.