Hi Ranjani, Defining the system using wraparound, you have access to the momentum dependent Hamiltonian as: syst.hamiltonian_submatrix(params=dict(kx=0, ky=0)). Then, you can apply various functionalities to it at different momenta. Regards, Ousmane