Hi dears

Based on the code for zincblende structure inside the Kwant website, I want to build a diamond lattice, but I don't know what I need to change in this code to create a diamond structure, also I need to create one lead for this structure.

can you help me?

Thank you in advance.

>>>import kwant
>>>from matplotlib import pyplot

>>>lat = kwant.lattice.general([(0, 0.5, 0.5), (0.5, 0, 0.5), (0.5, 0.5, 0)],
                                              [(0, 0, 0), (0.25, 0.25, 0.25)])
>>>a, b = lat.sublattices

>>>def make_cuboid(a=15, b=10, c=5):
              def cuboid_shape(pos):
                    x, y, z = pos
              return 0 <= x < a and 0 <= y < b and 0 <= z < c

>>>syst = kwant.Builder()
>>>syst[lat.shape(cuboid_shape, (0, 0, 0))] = None
>>>syst[lat.neighbors()] = None

>>>return syst

>>>def main():
>>># the standard plotting style for 3D is mainly useful for
>>># checking shapes:
>>>syst = make_cuboid()

>>>kwant.plot(syst)

>>># visualize the crystal structure better for a very small system
>>>syst = make_cuboid(a=1.5, b=1.5, c=1.5)

>>>def family_colors(site):
              return 'r' if site.family == a else 'g'

>>>kwant.plot(syst, site_size=0.18, site_lw=0.01, hop_lw=0.05)#,site_color=family_colors)

>>>if __name__ == '__main__':
       main()