Convergence rate for numerical computation of the lattice Green's function

Abstract

Flexible boundary condition methods couple an isolated defect to bulk through the bulk lattice Green's function. The inversion of the force-constant matrix for the lattice Green's function requires Fourier techniques to project out the singular subspace, corresponding to uniform displacements and forces for the infinite lattice. Three different techniques--relative displacement, elastic Green's function, and discontinuity correction--have different computational complexity for a specified numerical error. We calculate the convergence rates for elastically isotropic and anisotropic cases and compare them to analytic results. Our results confirm that the discontinuity correction is the most computationally efficient method to compute the lattice Green's function.

12 pages, 4 figures

References in corpus (1)

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

Cited by in corpus (4)

• Automatic numerical evaluation of vacancy-mediated transport for arbitrary crystals: Onsager coefficients in the dilute limit using a Green function approach
• Direct calculation of lattice Green function with arbitrary interactions for general crystals
• Lattice Green's function for crystals containing a planar interface
• Stable absorbing boundary conditions for molecular dynamics in general domains