Bond-order potential for simulations of extended defects in tungsten

We present a bond-order potential (BOP) for the bcc transition metal tungsten. The bond-order potentials are a real-space semiempirical scheme for the description of interatomic interactions based on the tight-binding approximation. In the hierarchy of atomic-scale-modeling methods the BOPs thus provide a direct bridge between electronic-structure and atomistic techniques. Two variants of the BOP were constructed and extensively tested against accurate first-principles methods in order to assess the potentials\u27 reliability and applicability. A comparison of the BOP with a central-force potential is used to demonstrate that a correct description of directional mixed covalent and metallic bonds is crucial for a successful and fully transferable model. The potentials are applied in studies of low-index surfaces, symmetrical tilt grain boundaries, and dislocations

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

Screw dislocations in bcc metals display non-planar cores at zero temperature which result in high lattice friction and thermally activated strain rate behavior. In bcc W, electronic structure molecular statics calculations reveal a compact, non-degenerate core with an associated Peierls stress between 1.7 and 2.8 GPa. However, a full picture of the dynamic behavior of dislocations can only be gained by using more efficient atomistic simulations based on semiempirical interatomic potentials. In this paper we assess the suitability of five different potentials in terms of static properties relevant to screw dislocations in pure W. As well, we perform molecular dynamics simulations of stress-assisted glide using all five potentials to study the dynamic behavior of screw dislocations under shear stress. Dislocations are seen to display thermally-activated motion in most of the applied stress range, with a gradual transition to a viscous damping regime at high stresses. We find that one potential predicts a core transformation from compact to dissociated at finite temperature that affects the energetics of kink-pair production and impacts the mechanism of motion. We conclude that a modified embedded-atom potential achieves the best compromise in terms of static and dynamic screw dislocation properties, although at an expense of about ten-fold compared to central potentials

Bridging the Gap between Quantum Mechanics and Large-Scale Atomistic Simulation

The prospect of modeling across disparate length and time scales to achieve a predictive multiscale description of real materials properties has attracted widespread research interest in the last decade. To be sure, the challenges in such multiscale modeling are many, and in demanding cases, such as mechanical properties or dynamic phase transitions, multiple bridges extending from the atomic level all the way to the continuum level must be built. Although often overlooked in this process, one of the most fundamental and important problems in multiscale modeling is that of bridging the gap between first-principles quantum mechanics, from which true predictive power for real materials emanates, and the large-scale atomistic simulation of thousands or millions of atoms, which is usually essential to describe the complex atomic processes that link to higher length and time scales. For example, to model single-crystal plasticity at micron length scales via dislocation-dynamics simulations that evolve the detailed dislocation microstructure requires accurate large-scale atomistic information on the mobility and interaction of individual dislocations. Similarly, modeling the kinetics of structural phase transitions requires linking accurate large-scale atomistic information on nucleation processes with higher length and time scale growth processes