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