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
def main(): # the standard plotting style for 3D is mainly useful for # checking shapes: syst = make_cuboid()
# visualize the crystal structure better for a very small system syst = make_cuboid(a=1.5, b=1.5, c=1.5)
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__':