Editing pair_style MEAM

Hello LAMMPS users,

I am interested in using the pair_style mean command to test a variety of potentials in the literature for Lithium. One such potential found in the article “Development of modified embedded atom method for a bcc metal: lithium” (2003) by Yuan et al. contains a modified analytical form for the embedding function (by subtracting a rho term).

I was searching the mean_force.F file in the lib directory, however I could not find where the fp(i), what I assume to be the embedding energy (or at least derivative of embedding energy) was evaluated.

Would someone who is more code-savvy mind helping me out?

Thank you,
Joseph Vella
PhD Student
Department of Chemical & Biological Engineering
Princeton University

You can ask the author, Greg Wagner (CCd).