code: def normal(self):
          
          normal = np.cross(self.lattice.matrix[0], self.lattice.matrix[1])
          normal /= np.linalg.norm(normal)
          return normal
docs: Calculates the surface normal vector of the slab