Hi,
> ------------------------------------------------------------------------ > In my case, I have a system with more than one mode. I have found some > lines that consider all the modes (at the end of this page : > https://kwant-project.org/doc/dev/tutorial/operators) > 2.7. Computing local quantities: densities and currents ... > <https://kwant-project.org/doc/dev/tutorial/operators> > kwant-project.org > 2.7. Computing local quantities: densities and currents¶ In the > previous tutorials we have mainly concentrated on calculating global > properties such as conductance ... > > > ------------------------------------------------------------------------ > wf_left=wf(0) > J_m_left=sum(J_m_bound(p)forpinwf_left) > J_z_left=sum(J_z_bound(p)forpinwf_left) > ------------------------------------------------------------------------ > But this is equivalent to calculate the expression > > > which is simply the sum over the current coming from each mode. > However, I would have rather use the following expression to find a > more general expression of the current density > > > > What is therefore the good way to find J? I have not found more > informations in the Kwant documentation. > > Best regards, > > Nicolas If you want to calculate the total current at a given energy you really should sum the contributions from each mode independently! You say that the "more general" expression for the current density is: but I'm not sure what you mean by that. In the above expression you are making some weird kind of mix of different modes. This is certainly not the expectation value of the current operator; are you maybe trying to calculate some kind of noise properties? Regards, Joe