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

Joseph A. Yasi, Dallas R. Trinkle

Research output: Contribution to journalArticlepeer-review


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.

Original languageEnglish (US)
Article number066706
JournalPhysical Review E - Statistical, Nonlinear, and Soft Matter Physics
Issue number6
StatePublished - Jun 14 2012

ASJC Scopus subject areas

  • Statistical and Nonlinear Physics
  • Statistics and Probability
  • Condensed Matter Physics


Dive into the research topics of 'Direct calculation of the lattice Green function with arbitrary interactions for general crystals'. Together they form a unique fingerprint.

Cite this