Radiation tolerance of nanocrystalline ceramics: insights from Yttria Stabilized Zirconia.

Materials for applications in hostile environments, such as nuclear reactors or radioactive waste immobilization, require extremely high resistance to radiation damage, such as resistance to amorphization or volume swelling. Nanocrystalline materials have been reported to present exceptionally high radiation-tolerance to amorphization. In principle, grain boundaries that are prevalent in nanomaterials could act as sinks for point-defects, enhancing defect recombination. In this paper we present evidence for this mechanism in nanograined Yttria Stabilized Zirconia (YSZ), associated with the observation that the concentration of defects after irradiation using heavy ions (Kr(+), 400 keV) is inversely proportional to the grain size. HAADF images suggest the short migration distances in nanograined YSZ allow radiation induced interstitials to reach the grain boundaries on the irradiation time scale, leaving behind only vacancy clusters distributed within the grain. Because of the relatively low temperature of the irradiations and the fact that interstitials diffuse thermally more slowly than vacancies, this result indicates that the interstitials must reach the boundaries directly in the collision cascade, consistent with previous simulation results. Concomitant radiation-induced grain growth was observed which, as a consequence of the non-uniform implantation, caused cracking of the nano-samples induced by local stresses at the irradiated/non-irradiated interfaces

Stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4-9 $\mu$s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3-4 $\mu$s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept for last-ditch runaway mitigation in ITER, SPARC, and other future burning-plasma tokamaks.Comment: Submitted to: Nuclear Fusion - 16 pages (4 supplementary), 11 figures (5 supplementary), 4 table

Impact of Homogeneous Strain On Uranium Vacancy Diffusion In Uranium Dioxide

We present a detailed mechanism of, and the effect of homogeneous strains on, the migration of uranium vacancies in UO2. Vacancy migration pathways and barriers are identified using density functional theory and the effect of uniform strain fields are accounted for using the dipole tensor approach. We report complex migration pathways and noncubic symmetry associated with the uranium vacancy in UO2 and show that these complexities need to be carefully accounted for to predict the correct diffusion behavior of uranium vacancies. We show that under homogeneous strain fields, only the dipole tensor of the saddle with respect to the minimum is required to correctly predict the change in the energy barrier between the strained and the unstrained case. Diffusivities are computed using kinetic Monte Carlo simulations for both neutral and fully charged state of uranium single and divacancies. We calculate the effect of strain on migration barriers in the temperature range 800–1800 K for both vacancy types. Homogeneous strains as small as 2% have a considerable effect on diffusivity of both single and divacancies of uranium, with the effect of strain being more pronounced for single vacancies than divacancies. In contrast, the response of a given defect to strain is less sensitive to changes in the charge state of the defect. Further, strain leads to anisotropies in the mobility of the vacancy and the degree of anisotropy is very sensitive to the nature of the applied strain field for strain of equal magnitude. Our results suggest that the influence of strain on vacancy diffusivity will be significantly greater when single vacancies dominate the defect structure, such as sintering, while the effects will be much less substantial under irradiation conditions where divacancies dominate

Method for long time scale simulations of solids: Application to crystal growth and dopant clustering

An important challenge in theoretical chemistry is the time scale problem. Atomic motion can be simulated directly by integrating Newton's equations over a time scale of nanoseconds, but most interesting chemical reactions take place on a time scale of seconds. We have developed a methodology to bridge this time scale gap using harmonic transition state theory suitable for solid systems. Possible reactive events and their rates are found with a saddle point finding method called the dimer method. When enough events are found, a kinetic Monte Carlo algorithm is used to choose which event occurs so that the system's position can be advanced in time. This technique has two major advantages over traditional kinetic Monte Carlo -- atoms do not have to map onto lattice sites for classification and kinetic events can be arbitrarily complicated. We have studied the homoepitaxial growth of aluminum and copper using an EAM potential at 80K with experimentally relevant deposition rates of monolayers per minute using a multiple time scale approach. Atomic deposition events are simulated directly with classical dynamics for several picoseconds until the incident energy has dissipated, and the long time between deposition events is simulated with the adaptive kinetic Monte Carlo method. Our simulations indicate that the Al( 100) surface grows much smoother then Cu( 100) at temperature between 0 and 80K due in part to long range multi atom processes which enable aluminum atoms to easily descend from atop islands. The high rate of such processes is due to their low activation energy, which is supported by density functional theory calculations, and the trend that processes involving more atoms tend to have larger prefactors and be favored by entropy. The scheme is efficient enough to model the evolution of systems with ab-initio forces as well, for which I will show an example of the breakup of dopant clusters in silicon