Paul C. M. Fossati

Thermophysical and anion diffusion properties of (Ux,Th1−x)O2

Using molecular dynamics, the thermophysical properties of the (Ux,Th1−x)O2 system have been investigated between 300 and 3600 K. The thermal dependence of lattice parameter, linear thermal expansion coefficient, enthalpy and specific heat at constant pressure is explained in terms of defect formation and diffusivity on the oxygen sublattice. Vegards law is approximately observed for solid solution thermal expansion below 2000 K. Different deviations from Vegards law above this temperature occur owing to the different temperatures at which the solid solutions undergo the superionic transition (2500–3300 K). Similarly, a spike in the specific heat, associated with the superionic transition, occurs at lower temperatures in solid solutions that have a high U content. Correspondingly, oxygen diffusivity is higher in pure UO2 than in pure ThO2. Furthermore, at temperatures below the superionic transition, oxygen mobility is notably higher in solid solutions than in the end members. Enhanced diffusivity is promoted by lower oxygen-defect enthalpies in (Ux,Th1−x)O2 solid solutions. Unlike in UO2 and ThO2, there is considerable variety of oxygen vacancy and oxygen interstitial sites in solid solutions generating a wide range of property values. Trends in the defect enthalpies are discussed in terms of composition and the lattice parameter of (Ux,Th1−x)O2.

Prediction and characterisation of radiation damage in fluorapatite

Molecular dynamics simulations, used in conjunction with a set of classical pair potentials, have been employed to examine simulated radiation damage cascades in the fluorapatite structure. Regions of damage have subsequently been assessed for their ability to recover and the effect that damage has on the important structural units defining the crystal structure, namely phosphate tetrahedra and calcium meta-prisms. Damage was considered by identifying how the phosphorous coordination environment changed during a collision cascade. This showed that PO4 units are substantially retained, with only a very small number of under or over coordinated phosphate units being observed, even at peak radiation damage. By comparison the damaged region of the material showed a marked change in the topology of the phosphate polyhedra, which polymerised to form chains up to seven units in length. Significantly, the fluorine channels characteristic of the fluorapatite structure and defined by the structures calcium meta-prisms stayed almost entirely intact throughout. This meant that the damaged region could be characterised as amorphous phosphate chains interlaced with regular features of the original undamaged apatite structure.

Pipe and grain boundary diffusion of He in UO2

Molecular dynamics simulations have been conducted to study the effects of dislocations and grain boundaries on He diffusion in [Formula: see text]. Calculations were carried out for the {1 0 0}, {1 1 0} and {1 1 1} [Formula: see text] edge dislocations, the screw [Formula: see text] dislocation and Σ5, Σ13, Σ19 and Σ25 tilt grain boundaries. He diffusivity as a function of distance from the dislocation core and grain boundaries was investigated for the temperature range 2300-3000 K. An enhancement in diffusivity was predicted within 20 Å of the dislocations or grain boundaries. Further investigation showed that He diffusion in the edge dislocations follows anisotropic behaviour along the dislocation core, suggesting that pipe diffusion occurs. An Arrhenius plot of He diffusivity against the inverse of temperature was also presented and the activation energy calculated for each structure, as a function of distance from the dislocation or grain boundary.

Cation ordering and oxygen transport behaviour in Sr1−3x/2LaxTiO3 perovskites

Molecular dynamics simulations and genetic algorithms are used to identify the mechanisms by which oxygen is transported through the Sr1−3x/2LaxTiO3 family of perovskites as a function of x. Across this compositional range the relative stability of ordered structures and random arrangements of cations and vacancies on the A sublattice is established. O and Ti Frenkel pair formation is then predicted for the considered compositions. These results show that Ti defects are more favourable for all La-containing compositions, but have a larger defect volume than O defects for low x. Oxygen diffusion at high temperature is determined using Molecular dynamics simulations. Two types of oxygen transport mechanisms are identified, that each contribute to oxygen diffusion in ordered structures. These mechanism involve transient O and Ti defects respectively, and provide means for O ions to swap positions. Depending on cation ordering, they can be correlated in space, which causes qualitatively different diffusion patterns.

Simulations of threshold displacement in beryllium

Atomic scale molecular dynamics simulations of radiation damage have been performed on beryllium. Direct threshold displacement simulations along a geodesic projection of directions were used to investigate the directional dependence with a high spatial resolution. It was found that the directionally averaged probability of displacement increases from 0 at 35 eV, with the energy at which there is a 50% chance of a displacement occurring is 70 eV and asymptotically approaching 1 for higher energies. This is, however, strongly directionally dependent with a 50% probability of displacement varying from 35 to 120 eV, with low energy directions corresponding to the nearest neighbour directions. A new kinetic energy dependent expression for the average maximum displacement of an atom as a function of energy is derived which closely matches the simulated data.