Dear Gusthavo, The script ldos = kwant.ldos(syst, energy=0.8, params=dict(k_x = 1.7)), provides you the ldos for a given energy but not for specific kpoint. To do so you need to get the indices for those allowed_kpoints you are looking for (condictions). For instance, you might use prop_modes, _ = flead.modes(energy=0.8, args=[params]) momentas_allowed = prop_modes.momenta[np.where(condictions)[0]] inco_wf = kwant.wave_function(fsyst, energy=0.8, args=[params])(0) wf_for_related_allowed_kpoints = inco_wf[momentas_allowed, :] kwant.plotter.density(fsyst, wf_for_related_allowed_kpoints, cmap='gist_heat_r') As Abbout said, If you share the script, it would be better. I hope this help best Le mer. 14 févr. 2024 à 22:18, Gusthavo Miranda < gusthavomsbrizolla@gmail.com> a écrit :

Dear Kwant users,

I am trying to compute localized density of states and the wavefunction for a specific energy and kpoint.

Looking at the tutorials, I understand that the LDOS can be computed as

ldos = kwant.ldos(syst, energy=0.8)

However, when I try to compute for a specific kpoint in the following manner

ldos = kwant.ldos(syst, energy=0.8, params=dict(k_x = 1.7))

the plot of the LDOS looks the same. Am I missing something? Is this the correct way of plotting for a specific kpoint?

Thanks in advance

Gusthavo