Thank you Axel for your help,
So after looking further in the literature I found that the function governing the system can be described as :
V(r)=F(r) VBm(r) + (1- F(r)) VZbl(r)
VBm : Bornmayer Potential
VZbl: Ziegler Biersack potential
So I just want to confirm that I understood correctly how to create the tabulated potential please. Should i generate a code to evaluate the potential energy for each pair of atoms (ex: Zr-O, Zr-Zr)...etc for different "r" by substituting in the mentioned above relation and then use pair_style table to assign the potential to lammps and then use pair write to check on the evaluated potential?
Also will each pair of atoms have a separated potential file or the whole potential of the system will be in the same file?
I apologize if it is a vague question, I am just trying to get through this field from scratch
Thanks in advance
On Saturday, June 24, 2017 7:05 PM, Axel Kohlmeyer <akohlmey@...24...> wrote: