Added interaction potential calculation