Atomistic simulations of dynamics of an edge dislocation and its interaction with a void in copper: a comparative study
Authors:
- Wu-Rong JianWJWu-Rong Jian
- Min ZhangMZMin Zhang
- Shuozhi XuSXShuozhi Xu
- Irene J. BeyerleinIBIrene J. Beyerlein
- 27 April 2020
- journal article
- research article
- Published by IOP Publishing in Modelling and Simulation in Materials Science and Engineering
- Vol. 28 (4) , 045004
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.Keywords
Funding Information
- National Science Foundation (ACI-1053575 CNS-1725797 DMR-1720256)
This publication has 85 references indexed in Scilit:
- An analysis on nanovoid growth in body-centered cubic single crystalline vanadiumComputational Materials Science, 2011
- Molecular dynamics simulations of dislocation interaction with voids in nickelComputational Materials Science, 2011
- Dislocation drag at the nanoscaleActa Materialia, 2010
- Interaction of an edge dislocation with voids in α-iron modelled with different interatomic potentialsJournal of Physics: Condensed Matter, 2008
- Stress orientation and relativistic effects on the separation of moving screw dislocationsPhysical Review B, 2008
- Molecular dynamics investigation of dislocation pinning by a nanovoid in copperPhysical Review B, 2005
- Modeling of dislocation–grain boundary interactions in FCC metalsJournal of Nuclear Materials, 2003
- Structural stability and lattice defects in copper:Ab initio, tight-binding, and embedded-atom calculationsPhysical Review B, 2001
- Dynamic Crack Processes via Molecular DynamicsPhysical Review Letters, 1996
- Measurement of the Stress and Strain on Specimens in an Electron MicroscopeJapanese Journal of Applied Physics, 1971