Hi, I guess, after my mail to the list, I made some headway to figure out the COULomb and LJ potentials. It seems in do_inbnf they are stored in the array ee[i] (grps->estat.ee[i]) instead of epot[i] where I was getting zeros for them. I will appreciate any correction, if I am thinking wrong. Pradip.