This code calculates the MSD from LAMMPS trajectories using ASE to get atom positiions - View it on GitHub