Dear all,

I want to discretize the continuous Hamiltonian using kwant.continuum.discretize, but I have problems in inputing my hamiltonian. In the example, we only have elements in H(0,0),H(0,1),H(1,0),H(1,1),H(2,2),H(2,3),H(3,2),H(3,3), but in my system (for example, bilayer graphene) I have elements in H(0,2),H(1,3),H(2,0),H(3,1). How to input my hamiltonian into it.

For example,
hamiltonian = """
+ M * kron(sigma_z, sigma_0)
+A * k_x * kron(sigma_0, sigma_x)
+A * k_y * kron(sigma_0, sigma_y)+ Tc*[[0, 0, 1, 0], [0, 0, 0, 1],[1, 0, 0, 0] ],[0, 1, 0, 0]]

"""

The term Tc*[[0, 0, 1, 0], [0, 0, 0, 1],[1, 0, 0, 0] ],[0, 1, 0, 0]] will cause errors.

Thank you very much.

Khani Hosein