Modelling oxygen self-diffusion in UO2 under pressure

Access to values for oxygen self-diffusion over a range of temperatures and pressures in UO2 is important to nuclear fuel applications. Here, elastic and expansivity data are used in the framework of a thermodynamic model, the cBΩ model, to derive the oxygen self-diffusion coefficient in UO2 over a range of pressures (0–10 GPa) and temperatures (300–1900 K). The significant reduction in oxygen self-diffusion as a function of increasing hydrostatic pressure, and the associated increase in activation energy, is identified

The effect of Nb on the corrosion and hydrogen pick-up of Zr alloys

Abstract Zr-Nb alloys are known to perform better in corrosion and hydrogen pick-up than other Zr alloys but the mechanism by which this happens is not well understood. Atomistic simulations using density functional theory of both tetragonal and monoclinic ZrO2 were performed, with intrinsic defects and Nb dopants. The overall defect populations with respect to oxygen partial pressure were calculated and presented in the form of Brouwer diagrams. Nb is found to favour 5 + in monoclinic ZrO2 at all partial pressures, but can exist in oxidation states ranging from 5 + to 3 + in the tetragonal phase. Nb5+ is charge balanced by Zr vacancies in both phases, suggesting that contrary to previous assumptions, Nb does not act as an n-type dopant in the oxide layer. Clusters containing oxygen vacancies were considered, Nb2+ was shown to exist in the tetragonal phase with a binding energy of 2.4 eV. This supports the proposed mechanism whereby low oxidation state Nb ions (2 + or 3+) charge balance the build-up of positive space-charge in the oxide layer, increasing oxygen vacancy and electron mobility, leading to near-parabolic corrosion kinetics and a reduced hydrogen pick-up. Previous experimental work has shown that tetragonal ZrO2 transforms to the monoclinic phase during transition, and that during transition a sharp drop in the instantaneous hydrogen pick-up fraction occurs. The oxidation of lower charge state Nb defects to Nb5+ during this phase change, and the consequent temporary n-doping of the oxide layer, is proposed as an explanation for the drop in hydrogen pick-up during transition

The effect of Sn-VO defect clustering on Zr alloy corrosion

Density functional theory simulations were used to study Sn defect clusters in the oxide layer of Zr-alloys. Clustering was shown to play a key role in the accommodation of Sn in ZrO2, with the {SnZr:VO}× bound defect cluster dominant at all oxygen partial pressures below 10-20 atm, above which Sn Zr × is preferred. {SnZr:VO}× is predicted to increase the tetragonal phase fraction in the oxide layer, due to the elevated oxygen vacancy concentration. As corrosion progresses, the transition to Sn Zr × , and resultant destabilisation of the tetragonal phase, is proposed as a possible explanation for the early first transition observed in Sn-containing Zr-Nb alloys

The effect of Nb on the corrosion and hydrogen pick-up of Zr alloys

Abstract Zr-Nb alloys are known to perform better in corrosion and hydrogen pick-up than other Zr alloys but the mechanism by which this happens is not well understood. Atomistic simulations using density functional theory of both tetragonal and monoclinic ZrO2 were performed, with intrinsic defects and Nb dopants. The overall defect populations with respect to oxygen partial pressure were calculated and presented in the form of Brouwer diagrams. Nb is found to favour 5 + in monoclinic ZrO2 at all partial pressures, but can exist in oxidation states ranging from 5 + to 3 + in the tetragonal phase. Nb5+ is charge balanced by Zr vacancies in both phases, suggesting that contrary to previous assumptions, Nb does not act as an n-type dopant in the oxide layer. Clusters containing oxygen vacancies were considered, Nb2+ was shown to exist in the tetragonal phase with a binding energy of 2.4 eV. This supports the proposed mechanism whereby low oxidation state Nb ions (2 + or 3+) charge balance the build-up of positive space-charge in the oxide layer, increasing oxygen vacancy and electron mobility, leading to near-parabolic corrosion kinetics and a reduced hydrogen pick-up. Previous experimental work has shown that tetragonal ZrO2 transforms to the monoclinic phase during transition, and that during transition a sharp drop in the instantaneous hydrogen pick-up fraction occurs. The oxidation of lower charge state Nb defects to Nb5+ during this phase change, and the consequent temporary n-doping of the oxide layer, is proposed as an explanation for the drop in hydrogen pick-up during transition

A theoretical study of intrinsic point defects and defect clusters in magnesium aluminate spinel

Point and small cluster defects in magnesium aluminate spinel have been studied from a first principles viewpoint. Typical point defects that occur during collision cascade simulations are cation anti-site defects, which have a small formation energy and are very stable, O and Mg split interstitials and vacancies. Isolated Al interstitials were found to be energetically unfavourable but could occur as part of a split Mg-Al pair or as a three atom-three vacancy Al ‘ring’ defect, previously observed in collision cascades using empirical potentials. The structure and energetics of the defects were investigated using density functional theory (DFT) and the results compared to simulations using empirical fixed-charge potentials. Each point defect was studied in a variety of supercell sizes in order to ensure convergence. It was found that empirical potential simulations significantly overestimate formation energies, but that the type and relative stability of the defects are well-predicted by the empirical potentials both for point defects and small defect clusters

Molecular dynamics modelling of radiation damage in normal, partly inverse and inverse spinels

The radiation response of perfect crystals of MgAl2O4, partially inverted MgGa2O4 and fully inverse MgIn2O4 were investigated using molecular dynamics. Dynamical cascades were initiated in these spinels over a range of trajectories with energies of 400 eV and 2 keV for the primary knock-on event. Collision cascades were set up on each of the cation and anion sublattices and were monitored up to 10 ps. Simulations in the normal MgAl2O4 spinel for the 2 keV energy regime resulted in similar defect structures as obtained at the post-threshold 400 eV energies, with little clustering occurring. The predominant defect configurations were split interstitials and cation antisites. For the inverse spinels, a much wider variety of lattice imperfections was observed. More defects were also produced due to the formation of interstitialvacancy cation chains and oxygen crowdions

Crystallization of a classical two-dimensional electron system: Positional and orientational orders

Crystallization of a classical two-dimensional one-component plasma (electrons interacting with the Coulomb repulsion in a uniform neutralizing positive background) is investigated with a molecular dynamics simulation. The positional and the orientational correlation functions are calculated for the first time. We have found an indication that the solid phase has a quasi-long-range (power-law) positional order along with a long-range orientational order. This indicates that, although the long-range Coulomb interaction is outside the scope of Mermin's theorem, the absence of ordinary crystalline order at finite temperatures applies to the electron system as well. The `hexatic' phase, which is predicted between the liquid and the solid phases by the Kosterlitz-Thouless-Halperin-Nelson-Young theory, is also discussed.Comment: 3 pages, 4 figures; Corrected typos; Double columne

Hydrogen accommodation in Zr second phase particles: Implications for H pick-up and hydriding of Zircaloy-2 and Zircaloy-4

Ab-initio computer simulations have been used to predict the energies associated with the accommodation of H atoms at interstitial sites in {\alpha}, {\beta}-Zr and Zr.M intermetallics formed with common alloying additions (M = Cr, Fe, Ni). Intermetallics that relate to the Zr2(Ni,Fe) second phase particles (SPPs) found in Zircaloy-2 exhibit favourable solution enthalpies for H. The intermetallic phases that relate to the Zr(Cr,Fe)2 SPPs, found predominantly in Zircaloy-4, do not offer favourable sites for interstitial H. It is proposed that Zr(Cr,Fe)2 particles may act as bridges for the migration of H through the oxide layer, whilst the Zr2(Ni,Fe)-type particles will trap the migrating H until these are dissolved or fully oxidised