The NEGF equations on a small molecule  (Datta Part B, Ch. 18, Eqs. 18.1–18.4, 18.30)

Give it a Hamiltonian H (hopping t = −1 eV between neighbours) and attach a source and a drain with broadening γ1, γ2. At each energy: GR = [E−H−Σ]−1, A = GRΓGA, Gn = GRΣinGA, T = Trace[Γ1GRΓ2GA]. Click the plots to move the energy marker.

T(E) at the marker–
Ĩ1 + Ĩ2 (Eq. 18.4, q/h units)–
Total current I (T=300 K)–
∫ Trace A dE / 2π (should be N)–
Eigenvalues of the isolated H (eV)–