For detailed quantum mechanical concepts, definitions and approaches see:

J. K. Tomfohr and O. F. Sankey, Phys. Status Solidi B 226, 115 (2001).
P. Orellana and F. Claro, Phys. Rev. Lett. 90, 178302 (2003).
A. P. Horsfield, D. R. Bowler, and A. J. Fisher, J. Phys.: Condens. Matter 16, L65 (2004).
G. Stefanucci and C.-O. Almbladh, Phys. Rev. B 69, 195318 (2004).
M. Di Ventra and T. N. Todorov, J. Phys.: Condens. Matter 16, 8025 (2004).
K. Burke, R. Car, and R. Gebauer, Phys. Rev. Lett. 94, 146803 (2005).
S. Kurth, G. Stefanucci, C.-O. Almbladh, A. Rubio, and E. K. U. Gross, Phys. Rev. B 72, 035308 (2005).
N. Bushong, N. Sai, and M. Di Ventra, Nano Lett. 5, 2569 (2005).
J. Pedersen and A. Wacker, Phys. Rev. B 72, 195330 (2005).
Y. Zhu, J. Maciejko, T. Ji, H. Guo, and J. Wang, Phys. Rev. B 71, 075317 (2005).
R. Baer, T. Seideman, S. Ilani, and D. Neuhauser, J. Chem. Phys. 120, 3387 (2004).
Copyright © 2019, Femto Program. All Rights Reserved.
Theoretical background
Dynamical quantum mechanical phenomena in nature are governed by Time-Dependent Schrödinger Equation (TDSE), which plays the same role as the Newton’s second law in classical mechanics
in mathematical sense,
There exist too many approach to solve TDSE such as time-dependent perturbation theory, Green function and evolution operator technique.
One may find the conduction phenomenon in a molecular wire a time-dependent process, i.e., the electrons are transported from one end to the other end. Therefore, TDSE can be employed to simulate this process, utilizing evolution of the system’s wave function from its initial state to final state:
For infinitesimal time interval, evolution operator is defined as
In practice, matrix representation of evolution operator can be used in some finite basis set, e.g., atom-centered Gaussian type orbitals. The below equation describes superposition  of the time-dependent wave function over eigenstates of the time-indepndent Hamiltonian, i.e., stationary states of the unperturbed Hamiltonian in the absence of applied voltage:
Using time-dependent perturbation theory, it can be shown that the coefficients vary with time as:
Once evolution operator is determined, the system’s wave function at any time can be computed and all physical properties will be available. On the other hand, the time dependence of observable A can be directly computed by Heisenberg equation: