This is an old revision of the document!
The implementation use the shape of tersoff potential published in PRB37(1988)6991, with aij=1.