Plastic anisotropy and dislocation trajectory in BCC metals.

Plasticity in body-centred cubic (BCC) metals at low temperatures is atypical, marked in particular by an anisotropic elastic limit in clear violation of the famous Schmid law applicable to most other metals. This effect is known to originate from the behaviour of the screw dislocations; however, the underlying physics has so far remained insufficiently understood to predict plastic anisotropy without adjustable parameters. Here we show that deviations from the Schmid law can be quantified from the deviations of the screw dislocation trajectory away from a straight path between equilibrium configurations, a consequence of the asymmetrical and metal-dependent potential energy landscape of the dislocation. We propose a modified parameter-free Schmid law, based on a projection of the applied stress on the curved trajectory, which compares well with experimental variations and first-principles calculations of the dislocation Peierls stress as a function of crystal orientation.

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.

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. Moreover, 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.

Dislocation core energies and core fields from first principles.

Ab initio calculations in bcc iron show that a 111 screw dislocation induces a short-range dilatation field in addition to the Volterra elastic field. This core field is modeled in anisotropic elastic theory using force dipoles. The elastic modeling thus better reproduces the atom displacements observed in ab initio calculations. Including this core field in the computation of the elastic energy allows deriving a core energy which converges faster with the cell size, thus leading to a result which does not depend on the geometry of the dislocation array used for the simulation.