Investigation of surface orientation dependent sputtering of Ag.

Sputtering of metal surfaces can be both a beneficial phenomenon, for instance in the coating industry, or an undesired side-effect, for instant materials subjected to irradiation. While the average sputtering yields are well known in common metals, recent studies have shown that the yields can depend on the crystallographic orientation of the surface much stronger than commonly appreciated. In this study, we investigate by computational means, molecular dynamics, the sputtering of single crystalline Ag surfaces under various incoming energies. The results at low and high energy are compared to experimental results for single crystalline Ag nanocubes of different orientations. We observe strong differences between the sputtering yields of different surface directions and ion energies. We analyze the results in terms of the atom cluster size of the sputtered materials, and show that the cluster size distribution is a key factor to understand the correspondence between simulations and experiments. At low energies mainly single atoms are sputtered, whereas at higher energies the sputtered material is mainly in atom clusters.

An empirical potential for simulating hydrogen isotope retention in highly irradiated tungsten.

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.

Revealing hidden defects through stored energy measurements of radiation damage.

With full knowledge of a material's atomistic structure, it is possible to predict any macroscopic property of interest. In practice, this is hindered by limitations of the chosen characterization techniques. For example, electron microscopy is unable to detect the smallest and most numerous defects in irradiated materials. Instead of spatial characterization, we propose to detect and quantify defects through their excess energy. Differential scanning calorimetry of irradiated Ti measures defect densities five times greater than those determined using transmission electron microscopy. Our experiments also reveal two energetically distinct processes where the established annealing model predicts one. Molecular dynamics simulations discover the defects responsible and inform a new mechanism for the recovery of irradiation-induced defects. The combination of annealing experiments and simulations can reveal defects hidden to other characterization techniques and has the potential to uncover new mechanisms behind the evolution of defects in materials.

Effects of precipitates and dislocation loops on the yield stress of irradiated iron.

Plastic deformation of crystalline materials is governed by the features of stress-driven motion of dislocations. In the case of irradiated steels subject to applied stresses, small dislocation loops as well as precipitates are known to interfere with the dislocation motion, leading to an increased yield stress as compared to pure crystals. We study the combined effect of precipitates and interstitial glissile [Formula: see text] dislocation loops on the yield stress of iron, using large-scale three-dimensional discrete dislocation dynamics simulations. Precipitates are included in the simulations using our recent multi-scale implementation [A. Lehtinen et al., Phys. Rev. E 93 (2016) 013309], where the strengths and pinning mechanisms of the precipitates are determined from molecular dynamics simulations. In the simulations we observe dislocations overcoming precipitates with an atypical Orowan mechanism which results from pencil-glide of screw segments in iron. Even if the interaction mechanisms with dislocations are quite different, our results suggest that in relative terms, precipitates and loops of similar sizes contribute equally to the yield stress in multi-slip conditions.

Improving atomic displacement and replacement calculations with physically realistic damage models.

Atomic collision processes are fundamental to numerous advanced materials technologies such as electron microscopy, semiconductor processing and nuclear power generation. Extensive experimental and computer simulation studies over the past several decades provide the physical basis for understanding the atomic-scale processes occurring during primary displacement events. The current international standard for quantifying this energetic particle damage, the Norgett-Robinson-Torrens displacements per atom (NRT-dpa) model, has nowadays several well-known limitations. In particular, the number of radiation defects produced in energetic cascades in metals is only ~1/3 the NRT-dpa prediction, while the number of atoms involved in atomic mixing is about a factor of 30 larger than the dpa value. Here we propose two new complementary displacement production estimators (athermal recombination corrected dpa, arc-dpa) and atomic mixing (replacements per atom, rpa) functions that extend the NRT-dpa by providing more physically realistic descriptions of primary defect creation in materials and may become additional standard measures for radiation damage quantification.

Multiscale modeling of dislocation-precipitate interactions in Fe: From molecular dynamics to discrete dislocations.

The stress-driven motion of dislocations in crystalline solids, and thus the ensuing plastic deformation process, is greatly influenced by the presence or absence of various pointlike defects such as precipitates or solute atoms. These defects act as obstacles for dislocation motion and hence affect the mechanical properties of the material. Here we combine molecular dynamics studies with three-dimensional discrete dislocation dynamics simulations in order to model the interaction between different kinds of precipitates and a 1/2〈111〉{110} edge dislocation in BCC iron. We have implemented immobile spherical precipitates into the ParaDis discrete dislocation dynamics code, with the dislocations interacting with the precipitates via a Gaussian potential, generating a normal force acting on the dislocation segments. The parameters used in the discrete dislocation dynamics simulations for the precipitate potential, the dislocation mobility, shear modulus, and dislocation core energy are obtained from molecular dynamics simulations. We compare the critical stresses needed to unpin the dislocation from the precipitate in molecular dynamics and discrete dislocation dynamics simulations in order to fit the two methods together and discuss the variety of the relevant pinning and depinning mechanisms.