code: def get_mi_vec(slab):
      
      mvec = np.cross(slab.lattice.matrix[0], slab.lattice.matrix[1])
      return mvec / np.linalg.norm(mvec)
docs: Convenience function which returns the unit vector aligned
      with the miller index.