Physics, Python, Python Programming, Coding, Programming.
0]) a2 = np.array([a/2, a*np.sqrt(3)/2]) # Nearest-neighbor vectors delta1 = a1 / 2 delta2 = (a1 + a2) / 2 delta3 = (a2 - a1) / 2 Next, we will define the Hamiltonian matrix for the tight-binding model: def graphene_hamiltonian(kx, ky, t=2.8): """ Calculate the Hamiltonian matrix for graphene using