for i in range(2, 6):
            particle = self.system.storage.getParticle(i)
            diff = []
            for j in range(3):
                x_i = particle.pos[j] - dmy_p[i - 2][j]