tersoff_zbl potential for SiC

Dear LAMMS users,

I’m trying to reproduce SiC Tersoff ZBL potential energy dependence from distance presentd in Fig. 1 of original publication by Devanathan et al. in JNM 253 (1998) 47 using metal units and original code pair_tersoff_zdl.cpp

I need coefficient 4 in front of ZBL potential to have an agreement with the figure.

lammps version is lammps-25Jul12.

Did anybody try to reproduce this dependence?

I would really appreciate any help!


German Samolyuk