Lattice Green function for extended defect calculations: Computation and error estimation with long-range forces

Dallas R. Trinkle

DOI: https://doi.org/10.1103/PhysRevB.78.014110

2006-07-15

Abstract:Computing the atomic geometry of lattice defects--point defects, dislocations, crack tips, surfaces, or boundaries--requires an accurate coupling of the local strain field to the long-range elastic field. Periodic 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. Recently, flexible boundary conditions have been developed to produce the correct long-range strain field from a defect--effectively "embedding" a defect in a finite cell 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 dynamical matrix, for supercell methods (periodic boundary conditions) it can only be calculated up to a maximum range. We illustrate how to accurately calculate the lattice Green function and estimate the error using a cutoff dynamical 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 quasi-continuum and flexible boundary-condition calculations, and measures the error introduced with periodic boundary conditions.

Materials Science

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is how to accurately couple the local strain field with the long - range elastic field when calculating the atomic geometric structures of lattice defects (such as point defects, dislocations, crack tips, surfaces or boundaries). Specifically, the paper focuses on the inability to accurately reproduce the correct bulk response of isolated defects when using periodic boundary conditions, especially when dealing with dislocations.

To solve the above problems, the author introduced flexible boundary conditions. The core of this is to accurately describe the long - range strain field around the defect by calculating the lattice Green function (LGF). The LGF can effectively "embed" the defect within a finite element and isolate it from the influence of periodic mirror images or free surfaces. However, for the supercell method (i.e., periodic boundary conditions), the LGF can only be calculated up to a maximum range. Therefore, it is necessary to combine the truncated dynamic matrix and the knowledge of the long - range behavior of the LGF to accurately calculate the LGF and estimate the error.

The key contribution of the paper is to provide a general and accurate method for calculating the LGF applicable to various defect geometric structures in density - functional theory (DFT). In addition, the author also proposed and tested an LGF error estimation method based on dynamic matrix calculation and elastic constants. This not only improves the calculation accuracy but also avoids the errors introduced due to the geometric limitations of periodic boundary conditions.

In summary, this paper aims to develop an effective method for calculating the lattice Green function and evaluating its error, thereby providing an important length scale for multi - scale quasi - continuum methods and flexible boundary condition calculations.