First-principles theory of correlated transport through nanojunctions

Robust conductance of dumbbell molecular junctions with fullerene anchoring groups

The conductance of a molecular wire connected to metallic electrodes is known to be sensitive to the atomic structure of the molecule-metal contact. This contact is to a large extent determined by the anchoring group linking the molecular wire to the metal. It has been found experimentally that a dumbbell construction with C(60) molecules acting as anchors yields more well-defined conductances as compared to the widely used thiol anchoring groups. Here, we use density functional theory to investigate the electronic properties of this dumbbell construction. The conductance is found to be stable against variations in the detailed bonding geometry and in good agreement with the experimental value of G=3×10(-4) G(0). Electron tunneling across the molecular bridge occurs via the lowest unoccupied orbitals of C(60) which are pinned close to the Fermi energy due to partial charge transfer. Our findings support the original motivation to achieve conductance values more stable towards changes in the structure of the molecule-metal contact leading to larger reproducibility in experiments.

Critical comparison of electrode models in density functional theory based quantum transport calculations

We study the performance of two different electrode models in quantum transport calculations based on density functional theory: parametrized Bethe lattices and quasi-one-dimensional wires or nanowires. A detailed account of implementation details in both the cases is given. From the systematic study of nanocontacts made of representative metallic elements, we can conclude that the parametrized electrode models represent an excellent compromise between computational cost and electronic structure definition as long as the aim is to compare with experiments where the precise atomic structure of the electrodes is not relevant or defined with precision. The results obtained using parametrized Bethe lattices are essentially similar to the ones obtained with quasi-one-dimensional electrodes for large enough cross-sections of these, adding a natural smearing to the transmission curves that mimics the true nature of polycrystalline electrodes. The latter are more demanding from the computational point of view, but present the advantage of expanding the range of applicability of transport calculations to situations where the electrodes have a well-defined atomic structure, as is the case for carbon nanotubes, graphene nanoribbons, or semiconducting nanowires. All the analysis is done with the help of codes developed by the authors which can be found in the quantum transport toolbox ALACANT and are publicly available.

Are Kohn-Sham conductances accurate?

We use Fermi-liquid relations to address the accuracy of conductances calculated from the single-particle states of exact Kohn-Sham (KS) density functional theory. We demonstrate a systematic failure of this procedure for the calculation of the conductance, and show how it originates from the lack of renormalization in the KS spectral function. In certain limits this failure can lead to a large overestimation of the true conductance. We also show, however, that the KS conductances can be accurate for single-channel molecular junctions and systems where direct Coulomb interactions are strongly dominant.

Influence of CO in the structural and electrical properties of Pt nanocontacts: comparison with H2 molecule addition

Density-functional theory is used to determine the atomic structure and the zero-bias conductance of Pt nanocontacts upon CO addition. Three stable atomic configurations of the nanocontacts, in conjunction with a statistical analysis of their electrical response, are enough to explain the main features of the experimental conductance histograms. Remarkably, the role of nonlocal effects when approaching the exchange-correlation potential within density-functional theory turns out to be crucial for CO molecular-based systems. Finally, a comparison with H(2) molecule addition on Pt nanocontacts is presented.

Nonequilibrium steady state transport via the reduced density matrix operator

We present a very simple model for numerically describing the steady state dynamics of a system interacting with continua of states representing a bath. Our model can be applied to equilibrium and nonequilibrium problems. For a one-state system coupled to two free electron reservoirs, our results match the Landauer formula for current traveling through a molecule. More significantly, we can also predict the nonequilibrium steady state population on a molecule between two out-of-equilibrium contacts. While the method presented here is for one-electron Hamiltonians, we outline how this model may be extended to include electron-electron interactions and correlations, an approach which suggests a connection between the conduction problem and the electronic structure problem.

Kondo effect and conductance of nanocontacts with magnetic impurities

We study the impact of dynamical correlations on the electronic structure and coherent transport properties of Cu nanocontacts hosting a single magnetic impurity (Ni, Co, Fe) in the contact region. The strong dynamical correlations of the impurity 3d electrons are fully taken into account by combining density-functional calculations with a dynamical treatment of the impurity 3d shell in the one-crossing approximation. We find that dynamical correlations give rise to the Kondo effect and lead to Fano features in the coherent transport characteristics similar to those observed in related experiments.

Electron transport in open systems from finite-size calculations: Examination of the principal layer method applied to linear gold chains

We describe the occurrence of computational artifacts when the principal layer method is used in combination with the cluster approximation for the calculation of electronic transport properties of nanostructures. For a one-dimensional gold chain, we observe an unphysical band in the band structure. The artificial band persists for large principal layers and for large buffer sizes. We demonstrate that the assumption of equality between Hamiltonian elements of neighboring layers is no longer valid and that a discontinuity is introduced in the potential at the layer transition. The effect depends on the basis set. When periodic boundary conditions are imposed and the k-space sampling is converged, the discontinuity disappears and the principal layer method can be correctly applied by using a linear combination of atomic orbitals as basis set.

Amine-gold linked single-molecule circuits: experiment and theory

A combination of theory and experiment is used to quantitatively understand the conductance of single-molecule benzenediamine-gold junctions. A newly developed analysis is applied to a measured junction conductance distribution, based on 59 000 individual conductance traces, which has a clear peak at 0.0064 G0 and a width of +/-47%. This analysis establishes that the distribution width originates predominantly from variations in conductance across different junctions rather than variations in conductance during junction elongation. Conductance calculations based on density functional theory (DFT) for 15 distinct junction geometries show a similar spread. We show explicitly that differences in local structure have a limited influence on conductance because the amine-Au bonding motif is well-defined and flexible, explaining the narrow distributions seen in the experiments. The minimal impact of junction structure on conductance permits an unambiguous comparison of calculated and measured conductance values and a direct assessment of the widely used DFT theoretical framework. The average calculated conductance (0.046 G0) is found to be seven times larger than experiment. This discrepancy is explained quantitatively in terms of electron correlation effects to the molecular level alignments in the junction.

Ab initio correlation effects on the electronic and transport properties of metal(II)-phthalocyanine-based devices

Using first-principles calculations in the framework of density functional theory, we investigated the electronic and transport properties of metal(II)-phthalocyanine (M(II)Pc) systems, both in a single-molecule configuration and in a model device geometry. In particular, using copper(II)-Pc and manganese(II)-Pc as prototypical examples, we studied how electronic correlations on the central metal ion influence the analysis of the electronic structure of the system and we demonstrated that the choice of the exchange-correlation functional, also beyond the standard local or gradient corrected level, is of crucial importance for a correct interpretation of the data. Finally, our electronic transport simulations have shown that M(II)Pc-based devices can act selectively as molecular conductors, as in the case of copper, or as spin valves, as in the case of manganese, demonstrating once more the great potential of these systems for molecular nanoelectronics applications.

Nonequilibrium GW approach to quantum transport in nano-scale contacts

Correlation effects within the GW approximation have been incorporated into the Keldysh nonequilibrium transport formalism. We show that GW describes the Kondo effect and the zero-temperature transport properties of the Anderson model fairly well. Combining the GW scheme with density functional theory and a Wannier function basis set, we illustrate the impact of correlations by computing the I-V characteristics of a hydrogen molecule between two Pt chains. Our results indicate that self-consistency is fundamental for the calculated currents, but that it tends to wash out satellite structures in the spectral function.

Nonequilibrium dynamics of correlated electron transfer in molecular chains

The relaxation dynamics of correlated electron transport along molecular chains is studied based on a substantially improved numerically exact path integral Monte Carlo approach. As an archetypical model, we consider a Hubbard chain containing two interacting electrons coupled to a bosonic bath. For this generalization of the ubiquitous spin-boson model, non-Boltzmann equilibrium distributions are found for many-body states. By mapping the multiparticle dynamics onto an isomorphic single particle motion, this phenomenon is shown to be sensitive to particle statistics and, due to its robustness, allows for new control schemes in designed quantum aggregates.