“…When a given voltage is applied, the current is allowed to flow across the system. The electric current ( I ) is further calculated using the Landauer approach, and this can be obtained from the integration of the transmission curve as where I ( V b ) represents the electric current under applied bias voltages, e is the electron charge, h is Planck’s constant, f ( E – μ L/R ) is the Fermi–Dirac distribution function of the L and R electrodes, μ L/R is the chemical potential, which can move up and down according to the Fermi energy E F , and T ( E , V b ) is the quantum transmission probability of the electrons, which can be given as follows where the coupling matrices are given as Γ L/R = i [∑ L/R – ∑ L/R † ] and the NEGFs for the scattering region given as G ( E , V b ) = [ E × S s – H s [ρ] – ∑ L ( E , V b ) – ∑ R ( E , V b )] −1 , where S s is the overlap matrix, H s is the Hamiltonian matrix, ∑ L/R = V s L/R g L/R V L/R s is the self-interaction energy, ∑ L/R is a molecule electrode that takes into account L/R electrodes in the central scattering region, g L/R is the surface L/R Green’s function, and V L/R s = V S L/R † is the coupling matrix between L/R electrodes and the scattering region. ,,, …”