Atomistic simulations of dynamics of an edge dislocation and its interaction with a void in copper: a comparative study

Abstract

Atomistic simulation methods are appropriate tools for investigating the dynamics of dislocations and their interactions with obstacles in metallic materials. In particular, molecular dynamics (MD) simulations have been widely employed on these two topics in the past several decades. However, even for the same type of simulation, the results can vary. While some of the quantitative differences may be due to the choices of interatomic potential and simulation cell size, they could similarly be attributed to choice of model settings, which have also differed substantially to date. In this paper, we carry out systematic MD simulations to study the effects of a few key model settings on the dynamics of an edge dislocation and its interaction with a void in copper. For a fixed interatomic potential, three modeling parameters, including applied loading mode, boundary conditions, and thermostat, are considered and their influences on the stress–strain response, the dislocation velocity, and the critical stress for a dislocation to bypass a void are compared. For a few select cases, we further examine the influence of temperature, strain rate, and simulation cell size. The results show that (i) compared with flexible boundary conditions, rigid boundary conditions result in greater stress oscillations in simulation cells of certain sizes; (ii) compared with the cases of no thermostat and a full thermostat, a partial thermostat provides better temperature control and lower friction on the dislocation core, respectively; and (iii) for dislocation–void interactions, the critical dislocation bypassing stress in shear loading can be appropriately determined with either a constant applied strain rate or a constant applied stress although the strain rate cannot be controlled in the latter. This analysis reveals that these three settings greatly influence the accuracy and interpretation of the results for the same type of simulation.

Effects of vacancy concentration on the edge dislocation motion in copper by atomic simulations

“…This potential is very suitable to reproduce the energy of various lattice defects in copper and can be used to predict the lattice distortion and local stress field. In this work, the accuracy of the local distortions and stress field effects caused by vacancies is critical [29,30]. In addition, the stacking fault energy (SFE) of copper calculated by this potential function (44.39 mJ m −2 ) is in good agreement with the experimental values (45 mJ m −2 ), which is important for correctly predicting the edge dislocation behaviors and mobility under high stresses [7,28].…”

Confined Layer Slip Process in Nanolaminated Ag and Two Ag/Cu Nanolaminates

“…In these layers, the x , y , and z axes denote coordinates for layer M , while , , and pertain to layer . We have applied PBCs along the x (or ) and y (or ) axes, while traction-free boundary conditions are implemented along the z (or ) axis [ 14 , 47 , 48 , 49 ]. The crystallographic orientations are shown in Table 1 .…”

Moment tensor potential for static and dynamic investigations of screw dislocations in bcc Nb

“…In experiments, tensile or shear deformation is performed at constant strain rate [80]. Moreover, our own MD simulations for the screw dislocation in bcc Nb and published work for an edge dislocation in fcc Cu [81] have shown that the constant-stress loading mode leads initially to shear-stress oscillations using PADs. For these reasons, we used PADs (containing 96 000 atoms) and the constantstrain loading mode under NVE conditions.…”