Harmonic interlayer potential with no higher-order contributions.

Description

Perl script for creating a table-style inter-atomic potential for using in lammps scripts. Contains the potential file itself as well. While the inter-atomic potential is not exactly harmonic, it results in a perfectly harmonic interaction between crystallographic layers in FCC crystal.

Steps to reproduce

Reading of a section in the following is a must: Zbigniew Kozioł, Stretched-exponential stress dynamics in chain of springs and masses model of crystals: analytical results and MD simulations (2025). Otherwise, the provided Perl script is ready for usage and contains explanations. The provided potential file may be used right away as well, as lammps input script.

Categories

Licence