Geometries of edge and mixed dislocations in bcc Fe from first principles calculations

Abstract

We use density functional theory (DFT) to compute the core structures of a_0[100](010) edge, a_0[100](011) edge, a_0/2[-1-11](1-10) edge, and a_0/2[111](1-10) 71 degree mixed dislocations in body-centered cubic (bcc) Fe. The calculations are performed using flexible boundary conditions (FBC), which effectively allow the dislocations to relax as isolated defects by coupling the DFT core to an infinite harmonic lattice through the lattice Green function (LGF). We use the LGFs of the dislocated geometries in contrast to most previous FBC-based dislocation calculations that use the LGF of the bulk crystal. The dislocation LGFs account for changes in the topology of the crystal in the core as well as local strain throughout the crystal lattice. The standard deviations of the dislocation Nye tensor distributions quantify the widths of the dislocation cores. The relaxed cores are compact, and the local magnetic moments on the Fe atoms closely follow the volumetric strain distributions in the cores. We also compute the core structures of these dislocations using eight different classical interatomic potentials, and quantify symmetry differences between the cores using the Fourier coefficients of their Nye tensor distributions. Most of the core structures computed using the classical potentials agree well with the DFT results. The DFT core geometries provide substitutional and interstitial sites for computing solute-dislocation interactions and can serve as inputs for mesoscale models of solute diffusion near dislocations.