Neighbor list now return r[j]-r[i] as distance vector (and corresponding shift vector). This corresponds to the definition used in libAtoms: dr_ij = pos(:,j) - pos(:,i) + (lattice .mult. shift)