Pair potential (Including Biersack and Ziegler)

V(R)=
(9.7342365892908E+03/r)*(0.1818*exp(-2.8616724320005E+01*r)+0.5099*exp(-8.4267310396064E+00*r)+0.2802*exp(-3.6030244464156E+00*r)+0.02817*exp(-1.8028536
321603E+00*r))*H0(1.0000-r)
+exp(6.4265260576348E+00+1.7900488524286E+00*r-4.5108316729807E+00*r^2+1.0866199373306E+00*r^3)*H(r-1.0000)*H0(2.0000-r)+(-2.4028204854115E+01*(H(2.2-r)*(2.2-r)^3)+1.1300691696477E+01*(H(2.3-r)*(2.3-r)^3)

+5.3144495820462E+00*(H(2.4-r)*(2.4-r)^3)
-4.6659532856049E+00*(H(2.5-r)*(2.5-r)^3)
+5.9637758529194E+00*(H(2.6-r)*(2.6-r)^3)
-1.7710262006061E+00*(H(2.7-r)*(2.7-r)^3)
+8.5913830768731E-01*(H(2.8-r)*(2.8-r)^3)
-2.1845362968261E+00*(H(3.0-r)*(3.0-r)^3)
+2.6424377007466E+00*(H(3.3-r)*(3.3-r)^3)
-1.0358345370208E+00*(H(3.7-r)*(3.7-r)^3)
+3.3548264951582E-01*(H(4.2-r)*(4.2-r)^3)
-4.6448582149334E-02*(H(4.7-r)*(4.7-r)^3)
-7.0294963048689E-03*(H(5.3-r)*(5.3-r)^3))*H(r-2.0000)


cutoff radius for (1-1) = 5.30 A

Density function:

phi(r)=
1.1686859407970E+01*(H(2.4-r)*(2.4-r)^3)-
1.4710740098830E-02*(H(3.2-r)*(3.2-r)^3)+
4.7193527075943E-01*(H(4.2-r)*(4.2-r)^3)

cutoff radius for (1-1) = 4.20 A


Embedding Function:

F(rho) = (-rho^0.5)-3.5387096579929E-04*(rho^2)