American Chemical Society, The Journal of Physical Chemistry A, 21(118), p. 3827-3834, 2014
DOI: 10.1021/jp501485a
Full text: Download
Ortiz and Seminario (J. Chem. Phys. 2007, 127, 111106/1-3) proposed some years ago a simple and direct approach to obtain I/V profiles from the combination of ab initio equilibrium electronic structure calculations and the uncertainty principle, as an alternative or complementary tool to more sophisticated non-equilibrium Green's functions methods. In this work, we revisit the fundamentals of this approach and reformulate accordingly the expression of the electric current. By analogy with the spontaneous electron decay process in electron transitions, in our revision the current is calculated upon the relaxing process from the "polarized" state induced by the external electric field to the electronic ground state. The electric current is obtained from the total charge transferred through the molecule and the corresponding electronic energy relaxation. The electric current expression proposed is more general compared to the previous expression employed by Ortiz and Seminario, where the charge variation must be tested among different slabs of atoms at the contact. This new approach has been tested on benzene-1,4-dithiolate attached to different gold clusters that represent the contact with the electrodes. Analysis of the total electron deformation density induced by the external electric voltage and properties associated to the electron deformation orbitals support the conclusions obtained from the I/V profiles.