Interplay of defect cluster and the stability of xenon in uranium dioxide by density functional calculations

Self-defect clusters in bulk matrix might affect the thermodynamic behavior of fission gases in nuclear fuel such as uranium dioxide. With first-principles LSDA+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 tri-vacancy sites to incorporate xenon in hyper-stoichiometric regime at high temperatures. The results show that in hypo-stoichiometric regime neutral tri-vacancy sites are the most favored position for diluted xenon gas, whereas in hyper-stoichiometric 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 meta-stable states difficulty that frequently encountered in DFT+U applications, namely, the quasi-annealing procedure, is also discussed.

💡 Research Summary

This paper investigates how intrinsic defect clusters in uranium dioxide (UO₂) influence the thermodynamic behavior of xenon (Xe), a representative fission gas, using first‑principles density‑functional theory (DFT) with the LSDA+U approach. The authors consider both hypo‑stoichiometric (U‑rich) and hyper‑stoichiometric (O‑rich) regimes, constructing realistic defect models that include uranium vacancies (U‑vacancies), tri‑vacancies (a uranium vacancy combined with two neighboring oxygen vacancies), oxygen interstitials, and oxygen clusters of the cuboctahedral (COT) type.

A major methodological contribution is the introduction of a “quasi‑annealing” procedure to overcome the notorious meta‑stable state problem in DFT+U calculations. By initially heating the electronic subsystem to a high‑energy configuration and then slowly cooling it, the electronic density is guided toward the true ground‑state manifold, ensuring consistent convergence for all defect‑gas complexes.

Energetic analysis reveals that the presence of oxygen defect clusters raises the solution energy of Xe by roughly 0.5 eV. This increase accounts for about 43 % of the energy gap between the two lowest‑energy configurations at 700 K, indicating that oxygen clustering has a pronounced destabilizing effect on dissolved xenon. The two competing configurations are: (i) Xe⁺ occupying a uranium vacancy, where the strong Madelung potential ionizes the xenon atom, and (ii) neutral Xe residing in a tri‑vacancy.

Temperature‑dependent free‑energy calculations show a clear crossover. At low temperatures (≈300 K) the entropy contribution is small, and neutral Xe in a tri‑vacancy is thermodynamically favored. As temperature rises, the entropy term grows and the free energy of the ionized Xe⁺ in a U‑vacancy drops, leading to a crossover near 600–700 K. Consequently, in the hyper‑stoichiometric regime, even when oxygen clusters are present, Xe prefers uranium vacancies at low temperatures; at higher temperatures the tri‑vacancy becomes competitive. In the hypo‑stoichiometric regime, tri‑vacancies remain the most stable site across the examined temperature range.

The ionization of Xe at a U‑vacancy is a direct consequence of the strong electrostatic (Madelung) field generated by the surrounding lattice and defect charge distribution. This finding implies that by controlling the fuel temperature one can deliberately switch the preferred xenon site, thereby influencing the nucleation and growth rate of gas bubbles. For instance, operating at temperatures that favor Xe⁺ in U‑vacancies could accelerate bubble formation, whereas maintaining lower temperatures that keep Xe neutral in spacious tri‑vacancies could suppress bubble coalescence.

The authors also discuss the broader implications for fuel performance modeling. Traditional fuel‑behavior codes often treat fission‑gas solubility as a function of bulk composition alone, neglecting the role of specific defect clusters. The present results suggest that accurate predictions of gas swelling, bubble nucleation, and fuel swelling require explicit inclusion of oxygen‑cluster effects and temperature‑dependent site preferences.

In summary, the study provides three key insights: (1) oxygen defect clusters significantly increase xenon solution energy, (2) there is a temperature‑driven competition between uranium vacancies (hosting ionized Xe⁺) and tri‑vacancies (hosting neutral Xe), and (3) the quasi‑annealing technique offers a reliable route to resolve meta‑stable electronic states in DFT+U calculations of complex defect systems. These findings advance our understanding of fission‑gas behavior in realistic operating conditions and open pathways for engineering fuel microstructures and temperature profiles to mitigate gas‑induced degradation.