Saturation estimates taken from the QCD-EAM-SCM results at 0.95 eV for v=0 and 0.65 eV for v=1.
# KAUFMANN
fast(x,A,W,E0) = A/2 * ( 1 + erf( 1 * ( x - E0 ) / W ) )
slow(x,A,g) = A * exp( -x / g )
both(x,A1,W,E0,A2,g) = fast( x, A1, W, E0 ) + slow( x, A2, g )
#v=0  w=0.189                            Estimate based on theory
#J=0 E0=0.70813    As=0.00224 y=0.09656           A0=0.49918
#J=1 E0=0.72008    As=0.00121 y=0.14011           A0=0.51885
#J=2 E0=0.72103    As=0.00079 y=0.25531           A0=0.53378

#v=1  w=0.164
#J=0 E0=0.34240	   As=0.16153 y=0.03899           A0=0.51974
#J=1 E0=0.35042    As=0.09242 y=0.05097           A0=0.551028
#J=2 E0=0.35437    As=0.06495 y=0.05898           A0=0.59408


# Rettner E0 and W parameters come from the published polynomial fit.
# RETTNER
#v0
v0E0(J) = 0.592 + 0.0338*J - 0.0079*J*J + 0.00022*J*J*J
v0W(J)  = 0.172 + 0.00056*J - 0.00033*J*J
v0Ret(x,A,J) = fast( x, A, v0W(J), v0E0(J) )

#v1
v1E0(J) = 0.299 + 0.02000*J - 0.00407*J*J
v1W(J)  = 0.128 + 0.00940*J - 0.00090*J*J
v1Ret(x,A,J) = fast( x, A, v1W(J), v1E0(J) )




# From the estimate by Wijzenbroek et al.
# HODGSON  / MURPHY
# v0J1 E0=0.713  w=0.224