!! STILL BUGS TO FIX !!
periodic boundaries