Assessment of interatomic potentials for atomistic analysis of static and dynamic properties of screw dislocations in W

Gaussian Process Regression for Materials and Molecules

We provide an introduction to Gaussian process regression (GPR) machine-learning methods in computational materials science and chemistry. The focus of the present review is on the regression of atomistic properties: in particular, on the construction of interatomic potentials, or force fields, in the Gaussian Approximation Potential (GAP) framework; beyond this, we also discuss the fitting of arbitrary scalar, vectorial, and tensorial quantities. Methodological aspects of reference data generation, representation, and regression, as well as the question of how a data-driven model may be validated, are reviewed and critically discussed. A survey of applications to a variety of research questions in chemistry and materials science illustrates the rapid growth in the field. A vision is outlined for the development of the methodology in the years to come.

Atomistic simulations of dislocation mobility in refractory high-entropy alloys and the effect of chemical short-range order

Refractory high-entropy alloys (RHEAs) are designed for high elevated-temperature strength, with both edge and screw dislocations playing an important role for plastic deformation. However, they can also display a significant energetic driving force for chemical short-range ordering (SRO). Here, we investigate mechanisms underlying the mobilities of screw and edge dislocations in the body-centered cubic MoNbTaW RHEA over a wide temperature range using extensive molecular dynamics simulations based on a highly-accurate machine-learning interatomic potential. Further, we specifically evaluate how these mechanisms are affected by the presence of SRO. The mobility of edge dislocations is found to be enhanced by the presence of SRO, whereas the rate of double-kink nucleation in the motion of screw dislocations is reduced, although this influence of SRO appears to be attenuated at increasing temperature. Independent of the presence of SRO, a cross-slip locking mechanism is observed for the motion of screws, which provides for extra strengthening for refractory high-entropy alloy system.

Using first-principles calculations to predict the mechanical properties of transmuting tungsten under first wall fusion power-plant conditions

Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low long-term disposal radioactive footprint. Yet, despite these relevant features, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. Significant neutron-induced transmutation happens in these tungsten components during nuclear fusion reactions, creating transmutant elements including Re, Os and Ta. Density functional theory (DFT) calculations that allow the calculation of defect and solute energetics are critical to better understand the behavior and evolution of tungsten-based materials in a fusion energy environment. In this study, we present a novel computational approach to perform DFT calculations on transmuting materials. In particular, we predict elastic and plastic mechanical properties (such as bulk modulus, shear modulus, ductility parameter, etc) on a variety of W-X compositions that result when pure tungsten is exposed to the EU-DEMO fusion first wall conditions for ten years.

Discrete shear band plasticity through dislocation activities in body-centered cubic tungsten nanowires

Shear band in metallic crystals is localized deformation with high dislocation density, which is often observed in nanopillar deformation experiments. The shear band dynamics coupled with dislocation activities, however, remains unclear. Here, we investigate the dynamic processes of dislocation and shear band in body-centered cubic (BCC) tungsten nanowires via an integrated approach of in situ nanomechanical testing and atomistic simulation. We find a strong effect of surface orientation on dislocation nucleation in tungsten nanowires, in which {111} surfaces act as favorite sites under high strain. While dislocation activities in a localized region give rise to an initially thin shear band, self-catalyzed stress concentration and dislocation nucleation at shear band interfaces cause a discrete thickening of shear band. Our findings not only advance the current understanding of defect activities and deformation morphology of BCC nanowires, but also shed light on the deformation dynamics in other microscopic crystals where jerky motion of deformation band is observed.

Interatomic potentials for modelling radiation defects and dislocations in tungsten

We have developed empirical interatomic potentials for studying radiation defects and dislocations in tungsten. The potentials use the embedded atom method formalism and are fitted to a mixed database, containing various experimentally measured properties of tungsten and ab initio formation energies of defects, as well as ab initio interatomic forces computed for random liquid configurations. The availability of data on atomic force fields proves critical for the development of the new potentials. Several point and extended defect configurations were used to test the transferability of the potentials. The trends predicted for the Peierls barrier of the [Formula: see text] screw dislocation are in qualitative agreement with ab initio calculations, enabling quantitative comparison of the predicted kink-pair formation energies with experimental data.

{110} Slip with {112} slip traces in bcc Tungsten

While propagation of dislocations in body centered cubic metals at low temperature is understood in terms of elementary steps on {110} planes, slip traces correspond often with other crystallographic or non-crystallographic planes. In the past, characterization of slip was limited to post-mortem electron microscopy and slip trace analysis on the sample surface. Here with in-situ Laue diffraction experiments during micro-compression we demonstrate that when two {110} planes containing the same slip direction experience the same resolved shear stress, sharp slip traces are observed on a {112} plane. When however the {110} planes are slightly differently stressed, macroscopic strain is measured on the individual planes and collective cross-slip is used to fulfill mechanical boundary conditions, resulting in a zig-zag or broad slip trace on the sample surface. We anticipate that such dynamics can occur in polycrystalline metals due to local inhomogeneous stress distributions and can cause unusual slip transfer among grains.

Comparison of analytic and numerical bond-order potentials for W and Mo

Bond-order potentials (BOPs) are derived from the tight-binding approximation and provide a linearly-scaling computation of the energy and forces for a system of interacting atoms. While the numerical BOPs involve the numerical integration of the response (Green's) function, the expressions for the energy and interatomic forces are analytical within the formalism of the analytic BOPs. In this paper we present a detailed comparison of numerical and analytic BOPs. We use established parametrizations for the bcc refractory metals W and Mo and test structural energy differences; tetragonal, trigonal, hexagonal and orthorhombic deformation paths; formation energies of point defects as well as phonon dispersion relations. We find that the numerical and analytic BOPs generally are in very good agreement for the calculation of energies. Different from the numerical BOPs, the forces in the analytic BOPs correspond exactly to the negative gradients of the energy. This makes it possible to use the analytic BOPs in dynamical simulations and leads to improved predictions of defect energies and phonons as compared to the numerical BOPs.