-
Simulated TEM imaging of a heavily irradiated metal
Authors:
D. R. Mason,
M. Boleininger,
J. Haley,
E. Prestat,
G. He,
F. Hofmann,
S. L. Dudarev
Abstract:
We recast the Howie-Whelan equations for generating simulated transmission electron microscope (TEM) images, replacing the dependence on local atomic displacements with atomic positions only. This allows very rapid computation of simulated TEM images for arbitrarily complex atomistic configurations of lattice defects and dislocations in the dynamical two beam approximation. Large scale massively-o…
▽ More
We recast the Howie-Whelan equations for generating simulated transmission electron microscope (TEM) images, replacing the dependence on local atomic displacements with atomic positions only. This allows very rapid computation of simulated TEM images for arbitrarily complex atomistic configurations of lattice defects and dislocations in the dynamical two beam approximation. Large scale massively-overlapping cascade simulations performed with molecular dynamics, are used to generate representative high-dose nanoscale irradiation damage in tungsten at room temperature, and we compare the simulated TEM images to experimental TEM images with similar irradiation and imaging conditions. The simulated TEM shows 'white-dot' damage in weak-beam dark-field imaging conditions, in line with our experimental observations and as expected from previous studies, and in bright-field conditions a dislocation network is observed. In this work we can also compare the images to the nanoscale lattice defects in the original atomic structures, and find that at high dose the white spots are not only created by small dislocation loops, but rather arise from nanoscale fluctuations in strains around curved sections of dislocation lines.
△ Less
Submitted 26 January, 2024;
originally announced January 2024.
-
Fast low-temperature irradiation creep driven by athermal defect dynamics
Authors:
Alexander Feichtmayer,
Max Boleininger,
Johann Riesch,
Daniel R. Mason,
Luca Reali,
Till Höschen,
Maximilian Fuhr,
Thomas Schwarz-Selinger,
Rudolf Neu,
Sergei L. Dudarev
Abstract:
The occurrence of high stress concentrations in reactor components is a still intractable phenomenon encountered in fusion reactor design. We observe and quantitatively model a non-linear high-dose radiation mediated microstructure evolution effect that facilitates fast stress relaxation in the most challenging low-temperature limit. In situ observations of a tensioned tungsten wire exposed to a h…
▽ More
The occurrence of high stress concentrations in reactor components is a still intractable phenomenon encountered in fusion reactor design. We observe and quantitatively model a non-linear high-dose radiation mediated microstructure evolution effect that facilitates fast stress relaxation in the most challenging low-temperature limit. In situ observations of a tensioned tungsten wire exposed to a high-energy ion beam show that internal stress of up to 2 GPa relaxes within minutes, with the extent and time-scale of relaxation accurately predicted by a parameter-free multiscale model informed by atomistic simulations. As opposed to conventional notions of radiation creep, the effect arises from the self-organisation of nanoscale crystal defects, athermally coalescing into extended polarized dislocation networks that compensate and alleviate the external stress.
△ Less
Submitted 24 January, 2024;
originally announced January 2024.
-
Atomistic Study of Irradiation-Induced Plastic and Lattice Strain in Tungsten
Authors:
Jintong Wu,
Daniel R. Mason,
Fredric Granberg
Abstract:
We demonstrate a practical way to perform decomposition of the elasto-plastic deformation directly from atomistic simulation snapshots. Through molecular dynamics simulations on a large single crystal, we elucidate the intricate process of converting plastic strain, atomic strain, and rigid rotation during irradiation. Our study highlights how prismatic dislocation loops act as initiators of plast…
▽ More
We demonstrate a practical way to perform decomposition of the elasto-plastic deformation directly from atomistic simulation snapshots. Through molecular dynamics simulations on a large single crystal, we elucidate the intricate process of converting plastic strain, atomic strain, and rigid rotation during irradiation. Our study highlights how prismatic dislocation loops act as initiators of plastic strain effects in heavily irradiated metals, resulting in experimentally measurable alterations in lattice strain. We show the onset of plastic strain starts to emerge at high dose, leading to the spontaneous emergence of dislocation creep and irradiation-induced lattice swelling. This phenomenon arises from the agglomeration of dislocation loops into a dislocation network. Furthermore, our numerical framework enables us to categorize the plastic transformation into two distinct types: pure slip events and slip events accompanied by lattice swelling. The latter type is particularly responsible for the observed divergence in interstitial and vacancy counts, and also impacts the behavior of dislocations, potentially activating non-conventional slip systems.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
An empirical potential for simulating hydrogen isotope retention in highly irradiated tungsten
Authors:
Daniel R. Mason,
Duc Nguyen-Manh,
Victor W. Lindblad,
Fredric G. Granberg,
Mikhail Yu. Lavrentiev
Abstract:
We describe the parameterization of a tungsten-hydrogen empirical potential designed for use with large-scale molecular dynamics simulations of highly irradiated tungsten containing hydrogen isotope atoms, and report test results. Particular attention has been paid to getting good elastic properties, including the relaxation volumes of small defect clusters, and to the interaction energy between h…
▽ More
We describe the parameterization of a tungsten-hydrogen empirical potential designed for use with large-scale molecular dynamics simulations of highly irradiated tungsten containing hydrogen isotope atoms, and report test results. Particular attention has been paid to getting good elastic properties, including the relaxation volumes of small defect clusters, and to the interaction energy between hydrogen isotopes and typical irradiation-induced defects in tungsten. We conclude that the energy ordering of defects changes with the ratio of H atoms to point defects, indicating that this potential is suitable for exploring mechanisms of trap mutation, including vacancy loop to plate-like void transformations.
△ Less
Submitted 18 May, 2023;
originally announced May 2023.
-
Microstructure of a heavily irradiated metal exposed to a spectrum of atomic recoils
Authors:
Max Boleininger,
Daniel R. Mason,
Andrea E. Sand,
Sergei L. Dudarev
Abstract:
At temperatures below the onset of vacancy migration, metals exposed to energetic ions develop dynamically fluctuating steady-state microstructures. Statistical properties of these microstructures in the asymptotic high exposure limit are not universal and vary depending on the energy and mass of the incident ions. We develop a model for the microstructure of an ion-irradiated metal under athermal…
▽ More
At temperatures below the onset of vacancy migration, metals exposed to energetic ions develop dynamically fluctuating steady-state microstructures. Statistical properties of these microstructures in the asymptotic high exposure limit are not universal and vary depending on the energy and mass of the incident ions. We develop a model for the microstructure of an ion-irradiated metal under athermal conditions, where internal stress fluctuations dominate the kinetics of structural evolution. The balance between defect production and recombination depends sensitively not only on the total exposure to irradiation, defined by the dose, but also on the energy of the incident particles. The model predicts the defect content in the high dose limit as an integral of the spectrum of primary knock-on atom energies, with the finding that low energy ions produce a significantly higher amount of damage than high energy ions.
△ Less
Submitted 12 October, 2022; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Volume of a dislocation network
Authors:
Max Boleininger,
Sergei L. Dudarev,
Daniel R. Mason,
Enrique Martínez
Abstract:
We derive a simple analytical line integral expression for the relaxation volume tensor of an arbitrary interconnected dislocation network. This quantity determines the magnitude of dislocation contribution to the dimensional changes and volumetric swelling of a material, and highlights the fundamental dual role of dislocations as sources of internal strain as well as carriers of plastic deformati…
▽ More
We derive a simple analytical line integral expression for the relaxation volume tensor of an arbitrary interconnected dislocation network. This quantity determines the magnitude of dislocation contribution to the dimensional changes and volumetric swelling of a material, and highlights the fundamental dual role of dislocations as sources of internal strain as well as carriers of plastic deformation. To illustrate applications of the method, we compute the relaxation volume of a stacking fault tetrahedron, a defect commonly occurring in fcc metals; the volume of an unusual tetrahedral configuration formed by the (a/2)<111> and a<001> dislocations in a bcc metal; and estimate the relative contribution of extended dislocations to the volume relaxation of heavily irradiated tungsten.
△ Less
Submitted 11 October, 2022; v1 submitted 10 March, 2022;
originally announced March 2022.
-
An estimate for thermal diffusivity in highly irradiated tungsten using Molecular Dynamics simulation
Authors:
Daniel R Mason,
Abdallah Reza,
Fredric Granberg,
Felix Hofmann
Abstract:
The changing thermal conductivity of an irradiated material is among the principal design considerations for any nuclear reactor, but at present few models are capable of predicting these changes starting from an arbitrary atomistic model. Here we present a simple model for computing the thermal diffusivity of tungsten, based on the conductivity of the perfect crystal and resistivity per Frenkel p…
▽ More
The changing thermal conductivity of an irradiated material is among the principal design considerations for any nuclear reactor, but at present few models are capable of predicting these changes starting from an arbitrary atomistic model. Here we present a simple model for computing the thermal diffusivity of tungsten, based on the conductivity of the perfect crystal and resistivity per Frenkel pair, and dividing a simulation into perfect and athermal regions statistically. This is applied to highly irradiated microstructures simulated with Molecular Dynamics. A comparison to experiment shows that simulations closely track observed thermal diffusivity over a range of doses from the dilute limit of a few Frenkel pairs to the high dose saturation limit at 3 displacements per atom (dpa).
△ Less
Submitted 24 November, 2021; v1 submitted 25 June, 2021;
originally announced June 2021.
-
Parameter-free quantitative simulation of high dose microstructure and hydrogen retention in ion-irradiated tungsten
Authors:
Daniel R. Mason,
Fredric Granberg,
Max Boleininger,
Thomas Schwarz-Selinger,
Kai Nordlund,
Sergei L. Dudarev
Abstract:
Hydrogen isotopes are retained in materials for fusion power applications, changing both hydrogen embrittlement and tritium inventory as the microstructure undergoes irradiation damage. But modelling of highly damaged materials - exposed to over 0.1 displacements per atom (dpa) - where asymptotic saturation is observed, for example tungsten facing the plasma in a fusion tokamak reactor, is difficu…
▽ More
Hydrogen isotopes are retained in materials for fusion power applications, changing both hydrogen embrittlement and tritium inventory as the microstructure undergoes irradiation damage. But modelling of highly damaged materials - exposed to over 0.1 displacements per atom (dpa) - where asymptotic saturation is observed, for example tungsten facing the plasma in a fusion tokamak reactor, is difficult because a highly damaged microstructure cannot be treated as weakly interacting isolated defect traps. In this paper we develop computational techniques to find the defect content in highly irradiated materials without adjustable parameters. First we show how to generate converged high dose (>1 dpa) microstructures using a combination of the creation-relaxation algorithm and molecular dynamics simulations of collision cascades. Then we make robust estimates of point defects and void regions with simple developments of the Wigner-Seitz decomposition of lattice sites. We use our estimates of the void surface area to predict the deuterium retention capacity of tungsten as a function of dose. This is then compared to 3He nuclear reaction analysis (NRA) measurements of tungsten samples self-irradiated at 290 K to different damage doses and exposed to deuterium plasma at low energy at 370 K. We show that our simulated microstructures give an excellent match to the experimental data, with both model and experiment showing 1.5-2.0 at.% deuterium retained in tungsten in the limit of high dose.
△ Less
Submitted 24 June, 2021;
originally announced June 2021.
-
Observation of transient and asymptotic driven structural states of tungsten exposed to irradiation
Authors:
Daniel R. Mason,
Suchandrima Das,
Peter M. Derlet,
Sergei L. Dudarev,
Andrew London,
Hongbing Yu,
Nicholas W. Phillips,
David Yang,
Kenichiro Mizohata,
Ruqing Xu,
Felix Hofmann
Abstract:
Combining spatially resolved X-ray Laue diffraction with atomic-scale simulations, we observe how ion-irradiated tungsten undergoes a series of non-linear structural transformations with increasing irradiation exposure. Nanoscale defect-induced deformations accumulating above 0.02 displacements per atom (dpa) lead to highly fluctuating strains at ~0.1 dpa, collapsing into a driven quasi-steady str…
▽ More
Combining spatially resolved X-ray Laue diffraction with atomic-scale simulations, we observe how ion-irradiated tungsten undergoes a series of non-linear structural transformations with increasing irradiation exposure. Nanoscale defect-induced deformations accumulating above 0.02 displacements per atom (dpa) lead to highly fluctuating strains at ~0.1 dpa, collapsing into a driven quasi-steady structural state above ~1 dpa. The driven asymptotic state is characterized by finely dispersed vacancy defects coexisting with an extended dislocation network, and exhibits positive volumetric swelling due to the creation of new crystallographic planes through self-interstitial coalescence, but negative lattice strain.
△ Less
Submitted 5 October, 2020; v1 submitted 30 July, 2020;
originally announced July 2020.
-
Statistical mechanics of kinks on a gliding screw dislocation
Authors:
M. Boleininger,
M. Gallauer,
S. L. Dudarev,
T. D. Swinburne,
D. R. Mason,
D. Perez
Abstract:
The ability of a body-centered cubic metal to deform plastically is limited by the thermally activated glide motion of screw dislocations, which are line defects with a mobility exhibiting complex dependence on temperature, stress, and dislocation segment length. We derive an analytical expression for the velocity of dislocation glide, based on a statistical mechanics argument, and identify an app…
▽ More
The ability of a body-centered cubic metal to deform plastically is limited by the thermally activated glide motion of screw dislocations, which are line defects with a mobility exhibiting complex dependence on temperature, stress, and dislocation segment length. We derive an analytical expression for the velocity of dislocation glide, based on a statistical mechanics argument, and identify an apparent phase transition marked by a critical temperature above which the activation energy for glide effectively halves, changing from the formation energy of a double kink to that of a single kink. The analysis is in quantitative agreement with direct kinetic Monte Carlo simulations.
△ Less
Submitted 27 May, 2020; v1 submitted 27 May, 2020;
originally announced May 2020.
-
Morphological analysis of 3d atom probe data using Minkowski functionals
Authors:
Daniel R Mason,
Andrew J London
Abstract:
We present a morphological analysis of atom probe data of nanoscale microstructural features, using methods developed by the astrophysics community to describe the shape of superclusters of galaxies. We describe second-phase regions using Minkowski functionals, representing the regions' volume, surface area, mean curvature and Euler characteristic. The alloy data in this work show microstructures…
▽ More
We present a morphological analysis of atom probe data of nanoscale microstructural features, using methods developed by the astrophysics community to describe the shape of superclusters of galaxies. We describe second-phase regions using Minkowski functionals, representing the regions' volume, surface area, mean curvature and Euler characteristic. The alloy data in this work show microstructures that can be described as sponge-like, filament-like, plate-like, and sphere-like at different concentration levels, and we find quantitative measurements of these features. To reduce user decision-making in constructing isosurfaces and to enhance the accuracy of the analysis a maximum likelihood based denoising filter was developed. We show that this filter performs significantly better than a simple Gaussian smoothing filter. We also interpolate the data using natural cubic splines, to refine voxel sizes and to refine the surface. We demonstrate that it is possible to find a mathematically well-defined, quantitative description of microstructure from atomistic datasets, to sub-voxel resolution, without user-tuneable parameters.
△ Less
Submitted 29 January, 2020;
originally announced January 2020.
-
Atomistic-Object Kinetic Monte Carlo simulations of irradiation damage in tungsten
Authors:
Daniel R. Mason,
Andrea E. Sand,
Sergei L. Dudarev
Abstract:
We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-t…
▽ More
We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-the-fly. Elastic interactions are handled within objects with empirical potentials at short distances, and between spatially distinct objects using the dipole tensor formalism. The model is shown to evolve mobile interstitial clusters in tungsten faster than an equivalent molecular dynamics simulation, even at elevated temperatures. We apply the model to the evolution of complex defects generated using molecular dynamics simulations of primary radiation damage in tungsten. We show that we can evolve defect structures formed in cascade simulations to experimentally observable timescales of seconds while retaining atomistic detail. We conclude that the first few nanoseconds of simulation following cascade initiation would be better performed using molecular dynamics, as this will capture some of the near-temperature-independent evolution of small highly-mobile interstitial clusters. We also conclude that, for the 20keV PKA cascades annealing simulations considered here, internal relaxations of sessile objects difficult to capture using conventional object KMC with idealised object geometries establish the conditions for long timescale evolution.
△ Less
Submitted 8 April, 2019;
originally announced April 2019.
-
Relaxation volumes of microscopic and mesoscopic irradiation-induced defects in tungsten
Authors:
Daniel R Mason,
Duc Nguyen-Manh,
Mihai-Cosmin Marinica,
Rebecca Alexander,
Sergei L Dudarev
Abstract:
The low energy structures of irradiation-induced defects have been studied in detail, as these determine the available modes by which a defect can diffuse or relax. As a result, there are many studies concerning the relative energies of possible defect structures, and empirical potentials are commonly fitted to or evaluated with respect to these energies. But recently [Dudarev et al Nuclear Fusion…
▽ More
The low energy structures of irradiation-induced defects have been studied in detail, as these determine the available modes by which a defect can diffuse or relax. As a result, there are many studies concerning the relative energies of possible defect structures, and empirical potentials are commonly fitted to or evaluated with respect to these energies. But recently [Dudarev et al Nuclear Fusion 2018], we have shown how to determine the stresses, strains and swelling of reactor components under irradiation from the elastic properties of ensembles of irradiation-induced defects. These elastic properties have received comparatively little attention. Here we evaluate relaxation volumes of irradiation-induced defects in tungsten computed with empirical potentials, and compare to density functional theory results where available. Different empirical potentials give different results, but some potential-independent trends in relaxation volumes can be identified. We show that the relaxation volume of small defects can be predicted to within 10% from their point-defect count. For larger defects we provide empirical fits for the relaxation volume of as a function of size. We demonstrate that the relaxation volume associated with a single primary-damage cascade can be estimated from the primary knock-on atom (PKA) energy. We conclude that while annihilation of vacancy- and interstitial- character defects will invariably reduce the total relaxation volume of the cascade debris, empirical potentials disagree whether coalescence of defects will reduce or increase the total relaxation volume.
△ Less
Submitted 17 December, 2018;
originally announced December 2018.
-
The electron elevator: excitations across the band gap via a dynamical gap state?
Authors:
A. Lim,
W. M. C. Foulkes,
A. P. Horsfild,
D. R. Mason,
A. Schleife,
E. W. Draeger,
A. A. Correa
Abstract:
We have used time-dependent density functional theory to study self-irradiated Si. We calculate the electronic stopping power of Si in Si by evaluating the energy transferred to the electrons per unit path length by an ion of kinetic energy from $1$ eV to $100 $ keV moving through the host. Electronic stopping is found to be significant below the threshold velocity normally identified with direct…
▽ More
We have used time-dependent density functional theory to study self-irradiated Si. We calculate the electronic stopping power of Si in Si by evaluating the energy transferred to the electrons per unit path length by an ion of kinetic energy from $1$ eV to $100 $ keV moving through the host. Electronic stopping is found to be significant below the threshold velocity normally identified with direct transitions across the band gap. A structured crossover at low velocity exists in place of a hard threshold. An analysis of the time dependence of the transition rates using coupled linear rate equations enables one of the excitation mechanisms to be clearly identified: a defect state induced in the gap by the moving ion acts like an elevator and carries electrons across the band gap.
△ Less
Submitted 29 September, 2015;
originally announced September 2015.
-
Non-Contact Measurement of Thermal Diffusivity in Ion-Implanted Nuclear Materials
Authors:
Felix Hofmann,
Daniel R. Mason,
Jeffrey K. Eliason,
Alexei A. Maznev,
Keith A. Nelson,
Sergei L. Dudarev
Abstract:
Knowledge of mechanical and physical property evolution due to irradiation damage is essential for the development of future fission and fusion reactors. Ion-irradiation provides an excellent proxy for studying irradiation damage, allowing high damage doses without sample activation. Limited ion-penetration-depth means that only few-micron-thick damaged layers are produced. Substantial effort has…
▽ More
Knowledge of mechanical and physical property evolution due to irradiation damage is essential for the development of future fission and fusion reactors. Ion-irradiation provides an excellent proxy for studying irradiation damage, allowing high damage doses without sample activation. Limited ion-penetration-depth means that only few-micron-thick damaged layers are produced. Substantial effort has been devoted to probing the mechanical properties of these thin implanted layers. Yet, whilst key to reactor design, their thermal transport properties remain largely unexplored due to a lack of suitable measurement techniques. Here we demonstrate non-contact thermal diffusivity measurements in ion-implanted tungsten for nuclear fusion armour. Alloying with transmutation elements and the interaction of retained gas with implantation-induced defects both lead to dramatic reductions in thermal diffusivity. These changes are well captured by our modelling approaches. Our observations have important implications for the design of future fusion power plants.
△ Less
Submitted 21 May, 2015;
originally announced May 2015.
-
Direct observation of size scaling and elastic interaction between nano-scale defects in collision cascades
Authors:
X. Yi,
A. E. Sand,
D. R. Mason,
M. A. Kirk,
S. G. Roberts,
K. Nordlund,
S. L. Dudarev
Abstract:
Using in-situ transmission electron microscopy, we have directly observed nano-scale defects formed in ultra-high purity tungsten by low-dose high energy self-ion irradiation at 30K. At cryogenic temperature lattice defects have reduced mobility, so these microscope observations offer a window on the initial, primary damage caused by individual collision cascade events. Electron microscope images…
▽ More
Using in-situ transmission electron microscopy, we have directly observed nano-scale defects formed in ultra-high purity tungsten by low-dose high energy self-ion irradiation at 30K. At cryogenic temperature lattice defects have reduced mobility, so these microscope observations offer a window on the initial, primary damage caused by individual collision cascade events. Electron microscope images provide direct evidence for a power-law size distribution of nano-scale defects formed in high-energy cascades, with an upper size limit independent of the incident ion energy, as predicted by Sand et al. [Eur. Phys. Lett., 103:46003, (2013)]. Furthermore, the analysis of pair distribution functions of defects observed in the micrographs shows significant intra-cascade spatial correlations consistent with strong elastic interaction between the defects.
△ Less
Submitted 10 March, 2015;
originally announced March 2015.
-
Elastic trapping of dislocation loops in cascades in ion-irradiated tungsten foils
Authors:
D R Mason,
X Yi,
M A Kirk,
S L Dudarev
Abstract:
Using _in situ_ transmission electron microscopy (TEM), we have observed nanometre scale dislocation loops formed when an ultra-high-purity tungsten foil is irradiated with a very low fluence of self-ions. Analysis of the TEM images has revealed the largest loops to be predominantly of prismatic 1/2<111> type and of vacancy character. The formation of such dislocation loops is surprising since iso…
▽ More
Using _in situ_ transmission electron microscopy (TEM), we have observed nanometre scale dislocation loops formed when an ultra-high-purity tungsten foil is irradiated with a very low fluence of self-ions. Analysis of the TEM images has revealed the largest loops to be predominantly of prismatic 1/2<111> type and of vacancy character. The formation of such dislocation loops is surprising since isolated loops are expected to be highly mobile, and should escape from the foil. In this work we show that the observed size and number density of loops can be explained by the fact that the loops are _not_ isolated - the loops formed in close proximity in the cascades interact with each other and with vacancy clusters, also formed in cascades, through long-range elastic fields, which prevent the escape of loops from the foil. We find that experimental observations are well reproduced by object Kinetic Monte Carlo simulations of evolution of cascades _only_ if elastic interaction between the loops is taken into account. Our analysis highlights the profound effect of elastic interaction between defects on the microstructural evolution of irradiated materials.
△ Less
Submitted 4 February, 2014;
originally announced February 2014.