Dear all,

I have a code for the Mobius structure, using python scripts. However, I want to construct a structure by adding hexagonal lattices in a Mobius structure. How can I do this using Kwant.

Thank you in advance.
Best,
Nuwan,
MSU, Starkville

P.S.: The python script for the mobius structure is as follows:


import numpy as np
import matplotlib.pyplot as plt
import matplotlib.tri as mtri

from mpl_toolkits.mplot3d import Axes3D  

fig = plt.figure(figsize=plt.figaspect(0.5))

u = np.linspace(0, 2.0 * np.pi, endpoint=True, num=50)
v = np.linspace(-0.5, 0.5, endpoint=True, num=10)
u, v = np.meshgrid(u, v)
u, v = u.flatten(), v.flatten()

x = (1 + 0.5 * v * np.cos(u / 2.0)) * np.cos(u)
y = (1 + 0.5 * v * np.cos(u / 2.0)) * np.sin(u)
z = 0.5 * v * np.sin(u / 2.0)

tri = mtri.Triangulation(u, v)

ax = fig.add_subplot(1, 2, 1, projection='3d')
ax.plot_trisurf(x, y, z, triangles=tri.triangles)
ax.set_zlim(-1, 1)

plt.show()