Direct calculation of the lattice Green function with arbitrary interactions for general crystals

Efficient computation of lattice defect geometries such as point defects, dislocations, disconnections, grain boundaries, interfaces, and free surfaces requires accurate coupling of displacements near the defect to the long-range elastic strain. Flexible boundary condition methods embed a defect in infinite harmonic bulk through the lattice Green function. We demonstrate an efficient and accurate calculation of the lattice Green function from the force-constant matrix for general crystals with an arbitrary basis by extending a method for Bravais lattices. New terms appear due to the presence of optical modes and the possible loss of inversion symmetry. By separately treating poles and discontinuities in reciprocal space, numerical accuracy is controlled at all distances. We compute the lattice Green function for a two-dimensional model with broken symmetry to elucidate the role of different coupling terms. The algorithm is generally applicable in two and three dimensions to crystals with arbitrary number of atoms in the unit cell, symmetry, and interactions.