Dear LAMMPS-community,
I need to implement a new potential in LAMMPS for a part of my project.
Since later on I want to make it available to the community I need your advice on how to do that properly.
The potential is the one described and deduced in the paper
Lazic I. and Thijsse B., Computational Materials Science 53, 483 (2012)
http://www.sciencedirect.com/science/article/pii/S0927025611004824
This is a potential for Al-O systems with:
-
a MEAM but Reference Free
-
a variable charge electrostatic potential (Charge Transfer Ionic Potential), that allows the atoms to determine their charge with a conjugate gradient algorithm on the basis of their local neighbourhood. Here the long range electrostatic potential is computed with a Particle-Particle Particle-Mesh method.
I guess that the problem of implementing this potential in LAMMPS can, therefore, be divided in first implementing the RF-MEAM and then insert the CTIP with all its details.
In the official release I have seen a MEAM potential, but from the documentation this is not reference free.
As for the CTIP, I couldn’t find anything with variable charges in the repository of LAMMPS, but I have heard that PPPM is already implemented.
My questions are:
-
Has anyone already implemented this potential (by chance in a not-yet-public release) ?
-
Is a version of the Reference-Free MEAM available somewhere (or is anybody developing it) ?
-
Where can I find the documentation for the PPPM method ?
-
Is anyone implementing a variable charge method similar to the one described in the paper above ?
-
Where can I find a “guide to LAMMPS programming” where I can see the style and sinthax I should conform to
(in view of making it available later, so not to rewrite it from scratch) ?
Or can somebody give me some advice on how to proceed ?
-
In what file the conjugate gradient algorithm is implemented
(so that I can use it and, in case, modify a copy of it for my purposes) ? -
Where can I find a flux diagram of LAMMPS source files so to understand what I need to modify ?
Does it exist ?
Thank you very much in advance for your answers
Best regards,
Daniele Scopece
EMPA, the Swiss Federal Laboratories for Materials Science and Technology, Duebendorf, Switzerland