We describe a first-principles method for calculating electronic structure, vibrational modes and frequencies, electron-phonon couplings, and inelastic electron transport properties of an atomic-scale device bridging two metallic contacts under nonequilibrium conditions. The method extends the density-functional codes SIESTA and TRANSIESTA that use atomic basis sets. The inelastic conductance characteristics are calculated using the nonequilibrium Green’s function formalism, and the electron-phonon interaction is addressed with perturbation theory up to the level of the self-consistent Born approximation. While these calculations often are computationally demanding, we show how they can be approximated by a simple and efficient lowest order expansion. Our method also addresses effects of energy dissipation and local heating of the junction via detailed calculations of the power flow. We demonstrate the developed procedures by considering inelastic transport through atomic gold wires of various lengths, thereby extending the results presented in Frederiksen et al. [Phys. Rev. Lett. 93, 256601 (2004)]. To illustrate that the method applies more generally to molecular devices, we also calculate the inelastic current through different hydrocarbon molecules between gold electrodes. Both for the wires and the molecules our theory is in quantitative agreement with experiments, and characterizes the system-specific mode selectivity and local heating.

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 (Ne≥1) with individual chemical potentials and electronic temperatures. We describe its novel methods for electrostatic gating, contour optimizations, and assertion of charge conservation, as well as the newly implemented algorithms for optimized and scalable matrix inversion, performance-critical pivoting, and hybrid parallelization. Additionally, a generic NEGF “post-processing” code (tbtrans/phtrans) for electron and phonon transport is presented with several novelties such as Hamiltonian interpolations, Ne≥1 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 106 atoms on workstation computers. The new features of both codes are demonstrated and bench-marked for relevant test systems

Computationally inexpensive approximations describing electron-phonon scattering in molecular-scale conductors are derived from the nonequilibrium Green’s function method. The accuracy is demonstrated with a first-principles calculation on an atomic gold wire. Quantitative agreement between the full nonequilibrium 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 nonequilibrium 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.

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.

Graphene can develop large magnetic moments in custom-crafted open-shell nanostructures such as triangulene, a triangular piece of graphene with zigzag edges. Current methods of engineering graphene nanosystems on surfaces succeeded in producing atomically precise open-shell structures, but demonstration of their net spin remains elusive to date. Here, we fabricate triangulenelike graphene systems and demonstrate that they possess a spin $S=1$ ground state. Scanning tunneling spectroscopy identifies the fingerprint of an underscreened $S=1$ Kondo state on these flakes at low temperatures, signaling the dominant ferromagnetic interactions between two spins. Combined with simulations based on the meanfield Hubbard model, we show that this $S=1$ $\pi$ paramagnetism is robust and can be turned into an $S=1/2$ state by additional H atoms attached to the radical sites. Our results demonstrate that $\pi$ paramagnetism of high-spin graphene flakes can survive on surfaces, opening the door to study the quantum behavior of interacting $\pi$ spins in graphene systems.

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 symmetries 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.