I am reading the calculation of charge conductivity by Non-Equilibrium Green Function (NEGF) method in this following paper.
Van-Nam Do, Non-equilibirum Green function method: theory and application in simulation of nanometer electronic devices, Adv. Nat. Sci.: Nanosci. Nanotechnol. 5 (2014) 033001 (21pp)
According to it, the charge conductivity is computed by the following formula.
$$I=\frac{e_{0}}{\hbar}\int\frac{dE}{2\pi}T(E)[n_L(E)-n_R(E)]$$ Here, $T(E)$ is the transimission coefficient; $n_L(E)$ and $n_R(E)$ are the occupation number of carrier at the energy value $E$ on left and right leads. These occupation number is usually computed by the Fermi-Dirac distribution function.
If the left and right leads are the same materials; then, $n_L(E)$ and $n_R(E)$ should have the same value. This means $[n_L(E)-n_R(E)]$ should be zero and the total charge current would be zero as well.
Taking the Cu-benzene-Cu nanowire as example, there would be no charge current flowing through it, if the current conductivity is computed by this formula.
Obviously, this is not right.
Would anyone please tell me what is wrong with my understanding about this formula? Taking the Cu-benzene-Cu nanowire as example, would anyone please tell me how to compute $n_L(E)$ and $n_R(E)$ to make sure the final charge conductivity of the nanowire is not zero?
Thank you in advance.