Fixed the potential bug. D_ij had the wrong numerical value in the
paper.
Prepared benchmark.
Simplified the potential expressions from python.