On the use of bond functions in molecular calculations