Modeling inelastic phonon scattering in atomic- and molecular-wire junctions

Computationally inexpensive approximations describing electron-phonon scattering in molecular-scale conductors are derived from the non-equilibrium Green's function method. The accuracy is demonstrated with a first principles calculation on an atomic gold wire. Quantitative agreement between the full non-equilibrium Green's function calculation and the newly derived expressions is obtained while simplifying the computational burden by several orders of magnitude. In addition, analytical models provide intuitive understanding of the conductance including non-equilibrium heating and provide a convenient way of parameterizing the physics. This is exemplified by fitting the expressions to the experimentally observed conductances through both an atomic gold wire and a hydrogen molecule.Comment: 5 pages, 3 figure

From tunneling to contact: Inelastic signals in an atomic gold junction

The evolution of electron conductance in the presence of inelastic effects is studied as an atomic gold contact is formed evolving from a low-conductance regime (tunneling) to a high-conductance regime (contact). In order to characterize each regime, we perform density functional theory (DFT) calculations to study the geometric and electronic structures, together with the strength of the atomic bonds and the associated vibrational frequencies. The conductance is calculated by first evaluating the transmission of electrons through the system, and secondly by calculating the conductance change due to the excitation of vibrations. As found in previous studies [Paulsson et al., Phys. Rev. B. 72, 201101(R) (2005)] the change in conductance due to inelastic effects permits to characterize the crossover from tunneling to contact. The most notorious effect being the crossover from an increase in conductance in the tunneling regime to a decrease in conductance in the contact regime when the bias voltage matches a vibrational threshold. Our DFT-based calculations actually show that the effect of vibrational modes in electron conductance is rather complex, in particular when modes localized in the contact region are permitted to extend into the electrodes. As an example, we find that certain modes can give rise to decreases in conductance when in the tunneling regime, opposite to the above mentioned result. Whereas details in the inelastic spectrum depend on the size of the vibrational region, we show that the overall change in conductance is quantitatively well approximated by the simplest calculation where only the apex atoms are allowed to vibrate. Our study is completed by the application of a simplified model where the relevant parameters are obtained from the above DFT-based calculations.Comment: 8 pages, 5 figure

Inelastic fingerprints of hydrogen contamination in atomic gold wire systems

We present series of first-principles calculations for both pure and hydrogen contaminated gold wire systems in order to investigate how such impurities can be detected. We show how a single H atom or a single H2 molecule in an atomic gold wire will affect forces and Au-Au atom distances under elongation. We further determine the corresponding evolution of the low-bias conductance as well as the inelastic contributions from vibrations. Our results indicate that the conductance of gold wires is only slightly reduced from the conductance quantum G0=2e^2/h by the presence of a single hydrogen impurity, hence making it difficult to use the conductance itself to distinguish between various configurations. On the other hand, our calculations of the inelastic signals predict significant differences between pure and hydrogen contaminated wires, and, importantly, between atomic and molecular forms of the impurity. A detailed characterization of gold wires with a hydrogen impurity should therefore be possible from the strain dependence of the inelastic signals in the conductance.Comment: 5 pages, 3 figures, Contribution to ICN+T2006, Basel, Switzerland, July-August 200

Improvements on non-equilibrium and transport Green function techniques: the next-generation transiesta

We present novel methods implemented within the non-equilibrium Green function code (NEGF) transiesta based on density functional theory (DFT). Our flexible, next-generation DFT-NEGF code handles devices with one or multiple electrodes ($N_e\ge1$) with individual chemical potentials and electronic temperatures. We describe its novel methods for electrostatic gating, contour opti- mizations, and assertion of charge conservation, as well as the newly implemented algorithms for optimized and scalable matrix inversion, performance-critical pivoting, and hybrid parallellization. Additionally, a generic NEGF post-processing code (tbtrans/phtrans) for electron and phonon transport is presented with several novelties such as Hamiltonian interpolations, $N_e\ge1$ electrode capability, bond-currents, generalized interface for user-defined tight-binding transport, transmission projection using eigenstates of a projected Hamiltonian, and fast inversion algorithms for large-scale simulations easily exceeding $10^6$ atoms on workstation computers. The new features of both codes are demonstrated and bench-marked for relevant test systems.Comment: 24 pages, 19 figure

Unified description of inelastic propensity rules for electron transport through nanoscale junctions

We present a method to analyze the results of first-principles based calculations of electronic currents including inelastic electron-phonon effects. This method allows us to determine the electronic and vibrational symmeties in play, and hence to obtain the so-called propensity rules for the studied systems. We show that only a few scattering states -- namely those belonging to the most transmitting eigenchannels -- need to be considered for a complete description of the electron transport. We apply the method on first-principles calculations of four different systems and obtain the propensity rules in each case.Comment: 4 pages, 4 figures, 1 table http://link.aps.org/abstract/PRL/v100/e22660

Inelastic scattering and local heating in atomic gold wires

We present a method for including inelastic scattering in a first-principles density-functional computational scheme for molecular electronics. As an application, we study two geometries of four-atom gold wires corresponding to two different values of strain, and present results for nonlinear differential conductance vs. device bias. Our theory is in quantitative agreement with experimental results, and explains the experimentally observed mode selectivity. We also identify the signatures of phonon heating.Comment: 4 pages, 3 figures; minor changes, updated figures, final version published in Phys. Rev. Let

Search for a Metallic Dangling-Bond Wire on nn-doped H-passivated Semiconductor Surfaces

We have theoretically investigated the electronic properties of neutral and $n$-doped dangling bond (DB) quasi-one-dimensional structures (lines) in the Si(001):H and Ge(001):H substrates with the aim of identifying atomic-scale interconnects exhibiting metallic conduction for use in on-surface circuitry. Whether neutral or doped, DB lines are prone to suffer geometrical distortions or have magnetic ground-states that render them semiconducting. However, from our study we have identified one exception -- a dimer row fully stripped of hydrogen passivation. Such a DB-dimer line shows an electronic band structure which is remarkably insensitive to the doping level and, thus, it is possible to manipulate the position of the Fermi level, moving it away from the gap. Transport calculations demonstrate that the metallic conduction in the DB-dimer line can survive thermally induced disorder, but is more sensitive to imperfect patterning. In conclusion, the DB-dimer line shows remarkable stability to doping and could serve as a one-dimensional metallic conductor on $n$-doped samples.Comment: 8 pages, 5 figure

Magnetic frustration and fractionalization in oligo(indenoindenes)

Poly(indenoindenes) are {\pi}-conjugated ladder carbon polymers with alternating hexagons and pentagons hosting one unpaired electron for each five-membered ring in the open-shell limit. Here we study the main magnetic interactions that are present in finite oligo(indenoindenes) (OInIn), classifying the six possible isomers in two different classes of three isomers each. One class can be rationalized by frustrated S = 1/2 Heisenberg chains, with ferromagnetic interactions between neighbour sites and antiferromagnetic interactions between the next neighbours. The other class is characterized by the more trivial antiferromagnetic order. Employing several levels of theory we further show that the ground state of one of the isomers is a valence-bond solid (VBS) of ferromagnetic dimers (S = 1). This is topologically similar to that of the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, which is known to show fractional S = 1/2 states at the edges