def METHOD_NAME(self):
    module_dir = os.path.dirname(os.path.abspath(__file__))
    with open(f"{module_dir}/surface_samples.json") as data_file:
        surface_properties = json.load(data_file)
    surface_energies, miller_indices = {}, {}
    for mpid in surface_properties:
        e_surf_list, miller_list = [], []
        for surface in surface_properties[mpid]["surfaces"]:
            e_surf_list.append(surface["surface_energy"])
            miller_list.append(surface["miller_index"])
        surface_energies[mpid] = e_surf_list
        miller_indices[mpid] = miller_list
    # In the case of a high anisotropy material
    # Nb: mp-8636
    latt_Nb = Lattice.cubic(2.992)
    # In the case of an fcc material
    # Ir: mp-101
    latt_Ir = Lattice.cubic(3.8312)
    # In the case of a hcp material
    # Ti: mp-72
    latt_Ti = Lattice.hexagonal(4.6000, 2.8200)
    self.ucell_Nb = Structure(
        latt_Nb,
        ["Nb", "Nb", "Nb", "Nb"],
        [[0, 0, 0], [0, 0.5, 0.5], [0.5, 0, 0.5], [0.5, 0.5, 0]],
    )
    self.wulff_Nb = WulffShape(latt_Nb, miller_indices["mp-8636"], surface_energies["mp-8636"])
    self.ucell_Ir = Structure(
        latt_Nb,
        ["Ir", "Ir", "Ir", "Ir"],
        [[0, 0, 0], [0, 0.5, 0.5], [0.5, 0, 0.5], [0.5, 0.5, 0]],
    )
    self.wulff_Ir = WulffShape(latt_Ir, miller_indices["mp-101"], surface_energies["mp-101"])
    self.ucell_Ti = Structure(
        latt_Ti,
        ["Ti", "Ti", "Ti"],
        [[0, 0, 0], [0.333333, 0.666667, 0.5], [0.666667, 0.333333, 0.5]],
    )
    self.wulff_Ti = WulffShape(latt_Ti, miller_indices["mp-72"], surface_energies["mp-72"])
    self.cube = WulffShape(Lattice.cubic(1), [(1, 0, 0)], [1])
    self.hex_prism = WulffShape(Lattice.hexagonal(2.63, 5.21), [(0, 0, 1), (1, 0, 0)], [0.35, 0.53])
    self.surface_properties = surface_properties