pair_style kolmogorov/crespi/z command
pair_style hybrid/overlay kolmogorov/crespi/z cutoff
pair_style hybrid/overlay kolmogorov/crespi/z 20.0 pair_coeff * * none pair_coeff 1 2 kolmogorov/crespi/z CC.KC C C pair_style hybrid/overlay rebo kolmogorov/crespi/z 14.0 pair_coeff * * rebo CH.airebo C C pair_coeff 1 2 kolmogorov/crespi/z CC.KC C C
The kolmogorov/crespi/z style computes the Kolmogorov-Crespi interaction potential as described in (Kolmogorov). An important simplification is made, which is to take all normals along the z-axis.
It is important to have a suffiently large cutoff to ensure smooth forces. Energies are shifted so that they go continously to zero at the cutoff assuming that the exponential part of Vij (first term) decays sufficiently fast. This shift is achieved by the last term in the equation for Vij above.
This potential is intended for interactions between two layers of graphene. Therefore, to avoid interaction between layers in multi-layered materials, each layer should have a separate atom type and interactions should only be computed between atom types of neighboring layers.
This potential must be used in combination with hybrid/overlay. Other interactions can be set to zero using pair_style none.
This fix is part of the USER-MISC package. It is only enabled if LAMMPS was built with that package. See the Build package doc page for more info.