[Kwant] Eigenvectors of H(k) for different k values

Hello :) I'm trying to calculate the Berry curvature. Is there a way to get access to the eigenvectors of H (k) for different k values just diagonalizing the following hamiltonian H_k? -----------Code------------------------------------------------------------------------------- wrapped = kwant.wraparound.wraparound(sys).finalized() def ham(sys,k_x,k_y=None, **params): k= momentum_to_lattice(sys, [k_x] if k_y is None else [k_x,k_y]) p=dict(zip(sys._momentum_names,k),**params) return sys.hamiltonian_submatrix(params=p,sparse=False) H_k= ham(wrapped,kx,ky) ------------------------------------------------------------------------------------------------------ Thanks for your time and help! Nayra