An Improved Version of the Green's Function Molecular Dynamics Method.
Ling Ti Kong,Colin Denniston,Martin H. Mueser
DOI: https://doi.org/10.1016/j.cpc.2010.10.006
IF: 4.717
2010-01-01
Computer Physics Communications
Abstract:This work presents an improved version of the Green's function molecular dynamics method (Kong et al., 2009; Campana and Maser, 2004 [1,2]), which enables one to study the elastic response of a three-dimensional solid to an external stress field by taking into consideration only atoms near the surface. In the previous implementation, the effective elastic coefficients measured at the Gamma-point were altered to reduce finite size effects: their eigenvalues corresponding to the acoustic modes were set to zero. This scheme was found to work well for simple Bravais lattices as long as only atoms within the last layer were treated as Green's function atoms. However, it failed to function as expected in all other cases. It turns out that a violation of the acoustic sum rule for the effective elastic coefficients at Gamma (Kong, 2010 [3]) was responsible for this behavior. In the new version, the acoustic sum rule is enforced by adopting an iterative procedure, which is found to be physically more meaningful than the previous one. In addition, the new algorithm allows one to treat lattices with bases and the Green's function slab is no longer confined to one layer.New version program summaryProgram title: FixGFC/FixGFMD v1.12Catalogue identifier: AECW_v1_1Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AECW_v1_1.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 206 436No. of bytes in distributed program, including test data, etc.: 4 314 850Distribution format: tar.gzProgramming language: C++Computer: AllOperating system: LinuxHas the code been vectorized or parallelized?: Yes. Code has been parallelized using MPI directives. RAM: Depends on the problemClassification: 7.7External routines: LAMMPS (http://lammps.sandia.gov/), MPI (http://www.mcs.anl.gov/research/projects/mpi/), FFT (http://www.fftw.org/)Catalogue identifier of previous version: AECW_v1_0Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1004Does the new version supersede the previous version?: YesNature of problem: Green's function molecular dynamics (GFMD) is a coarse-graining method that enables one to investigate the full elastic response of an interface between a semi-infinite solid and a contact while taking only the surface atoms in the solid into consideration. The effect of long-range elastic deformations on the surface atoms from the semi-infinite solid is replaced by effective elastic interactions, thus reducing the problem from three dimensions to two dimensions without compromising the physical essence of the problem.