Help needed for potentials/Au_u3.eam

Hello everyone!
I have some questions about the well-known eam potential file Au_u3.eam (SM Foiles et al, PRB, 33, 7983 (1986)) in potentials package. Can anyone help me?

  1. I drew a picture of the effective charge function Z® versus distance r. And I think value of Z® is listed between line 104 and line 203 in Au_u3.eam potential file. I am confused with the value of Z® which ranges from 0 to 10.8. But according to FIG.1 of the original paper PRB, 33, 7983 (1986), Z® seems between 0 and 1. Is this due to the different units of Z®? My picture of Z® is shown below. Value of Z® is the Y axis and value of r is the X axis. Z_Auu3_1.jpg
  2. I also want to get the pair potential term phi® through Z®. According to the doc page of pair_eam ( Should I compute phi® like this: phi(i)=27.20.529Z(i)*Z(i)/r(i) and i=1,2,3, …, Nr, Nr is the number of tabulated values of Z® which is 500 in Au_u3.eam file. Two pictures of phi® versus distance r are shown below. I am not sure whether my calculation of phi® is correct. Value of phi® is the Y axis and value of r is the X axis. phi_Auu3_1.jpgphi_Auu3_2.jpg
    Any help would be appreciated!

Best Regards

I can’t tell you anything further than the doc page
and the code (src/pair_eam.cpp). The doc page
is quite explicit about what the units in the various

kinds of EAM file tabulations are.