implemented pre calculated epsilon geometric mean
reflection boundarys are now working with type autonomous lennard-jones-potential