Playing with my code (graphene lattice), I see that if I double my hopping parameter, the density in my colourplot is reduced by 2, which makes sense according to the units. However I'm not finding any changes when varying either the scattering region size
or the lattice parameter. I thought that when we said volume (or area in my case) it was related to the unit cell volume. Is it or is for example the area of the scattering region or something else?
Thanks again for your time.
I've found in other threads in the mailing list that the units of current is for example (unit of charge)/(hbar/unit of energy) (https://firstname.lastname@example.org/msg01100.html). Also, the local density of states has units of energy/volume (https://email@example.com/msg00169.html).
My question is, what is the units of the output of the density operator? Is it energy/volume as well?
I ask this because intuitively I view it as the square of the wavefunction, but it gives me values larger than 1 for each site when there is only 1 mode involved (see attached picture) so it is not just the probability of findinge the electron at that site because this should be maximum 1. I have also noticed that the values I get in the colorbar depend on the value of my hopping (e.g. case of graphene), but overall I'm not so sure of the units.