I am performing an NEB calculation with VASP and VTST. I would like to visualize the atoms along the path and to calculate their displacements. When I have a single relaxation's XDATCAR, I can unwrap the trajectories using ASE and get the displacements: steps = read(xdatcar, ':', format='vasp-xda...