Thank you for your help. The new potential function is as follows.
where, f(p)=exp(-C5*(p-C6)^2). r is the distance between atoms. C1, C2, C3, C4, C5, C6 are the constant values.
How can I use this potential function in Lammps?
I know that "pair style hybrid" command can make different potential functions add together, but I do not know how to multiply one function by another function.
University of Wisconsin-Madison, USA