TY - JOUR

T1 - Lattice Green function for extended defect calculations

T2 - Computation and error estimation with long-range forces

AU - Trinkle, D. R.

N1 - Copyright:
Copyright 2008 Elsevier B.V., All rights reserved.

PY - 2008/7/23

Y1 - 2008/7/23

N2 - Computing the atomic geometry of lattice defects-e.g., point defects, dislocations, crack tips, surfaces, or boundaries-requires an accurate coupling of the local deformations to the long-range elastic field. Periodic or fixed boundary conditions used by classical potentials or density-functional theory may not accurately reproduce the correct bulk response to an isolated defect; this is especially true for dislocations. Flexible boundary conditions have been developed to produce the correct long-range strain field from a defect-effectively "embedding" a finite-sized defect with infinite bulk response, isolating it from either periodic images or free surfaces. Flexible boundary conditions require the calculation of the bulk response with the lattice Green function (LGF). While the LGF can be computed from the force-constant matrix, the force-constant matrix is only known to a maximum range. This paper illustrates how to accurately calculate the lattice Green function and estimate the error using a truncated force-constant matrix combined with knowledge of the long-range behavior of the lattice Green function. The effective range of deviation of the lattice Green function from the long-range elastic behavior provides an important length scale in multiscale quasicontinuum and flexible boundary-condition calculations, and measures the error introduced with periodic-boundary conditions.

AB - Computing the atomic geometry of lattice defects-e.g., point defects, dislocations, crack tips, surfaces, or boundaries-requires an accurate coupling of the local deformations to the long-range elastic field. Periodic or fixed boundary conditions used by classical potentials or density-functional theory may not accurately reproduce the correct bulk response to an isolated defect; this is especially true for dislocations. Flexible boundary conditions have been developed to produce the correct long-range strain field from a defect-effectively "embedding" a finite-sized defect with infinite bulk response, isolating it from either periodic images or free surfaces. Flexible boundary conditions require the calculation of the bulk response with the lattice Green function (LGF). While the LGF can be computed from the force-constant matrix, the force-constant matrix is only known to a maximum range. This paper illustrates how to accurately calculate the lattice Green function and estimate the error using a truncated force-constant matrix combined with knowledge of the long-range behavior of the lattice Green function. The effective range of deviation of the lattice Green function from the long-range elastic behavior provides an important length scale in multiscale quasicontinuum and flexible boundary-condition calculations, and measures the error introduced with periodic-boundary conditions.

UR - http://www.scopus.com/inward/record.url?scp=48449097044&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=48449097044&partnerID=8YFLogxK

U2 - 10.1103/PhysRevB.78.014110

DO - 10.1103/PhysRevB.78.014110

M3 - Article

AN - SCOPUS:48449097044

VL - 78

JO - Physical Review B-Condensed Matter

JF - Physical Review B-Condensed Matter

SN - 0163-1829

IS - 1

M1 - 014110

ER -