use eq from mdtraj