Strongly correlated physics arises from electron-electron scattering within partially filled orbitals. Organic molecules in open-shell configurations are therefore good candidates to exhibit many-body effects. We focus on electron transport in a two-terminal single-molecule junction setup, in which the molecular bridge consists of an organic radical with a molecular orbital hosting a single unpaired electron [singly occupied molecular orbital (SOMO)]. We perform beyond state-of-the-art numerical simulations combining an ab initio description of the chemical environment with quantum field-theoretical techniques that account for many-body effects. The key observation is that the SOMO resonance is prone to splitting, and we identify a giant electronic scattering rate as the driving many-body mechanism, akin to that of the Mott metal-to-insulator transition. By comparing linear and cyclic radicals, we show that the spatial distribution of the SOMO and its projection on the molecular backbone have dramatic consequences for the transport properties of the junction. We argue that the phenomenon and the underlying microscopic mechanism apply to a broad family of open-shell molecular systems, and can explain puzzling experimental observations such as suppressed conductance in radical junctions.

We investigate the stability of destructive quantum interference (DQI) in electron transport through graphene nanostructures connected to source and drain electrodes. The fingerprint of DQI is an antiresonance in the transmission function, and its origin is deeply connected to the topology of the atomic structure, which we discuss in terms of symmetry arguments supported by numerical simulations. A systematic analysis of the transmission function versus system size reveals that the DQI antiresonance persists for large systems in the ballistic regime and establishes the quantum confinement gap as the intrinsic resolution limit to detect QI effects. Furthermore, we consider the influence of disorder, electron–electron and electron–phonon interactions, and provide quantitative criteria for the robustness of DQI in their presence. We find that the conductance is quite sensitive to perturbations, and its value alone may not be sufficient to characterize the quantum interference properties of a junction. Instead, the characteristic behavior of the transmission function is more resilient, and we suggest it retains information on the presence of an antiresonance even if DQI is partially concealed or suppressed. At the same time, DQI results in a non-linear transport regime in the current-bias characteristics that can be possibly detected in transport experiments.

Graphene-based sensors are exceptionally sensitive with high carrier mobility and low intrinsic noise, and have been intensively investigated in the past decade. The detection of individual gas molecules has been reported, albeit the underlying sensing mechanism is not yet well understood. We focus on the adsorption of NO 2 , H 2 O, and NH 3 on a molecular junction with a pyrene core, which can be considered as a minimal graphene-like unit. We systematically investigate the chemiresistive response within the framework of density functional theory and non-equilibrium Green’s functions. We highlight the fundamental role of quantum interference (QI) in the sensing process, and we propose it as a paradigmatic mechanism for sensing. Owing to the open-shell character of NO 2 , its interaction with pyrene gives rise to a Fano resonance thereby triggering the strongest chemiresistive response, while the weaker interactions with H 2 O and NH 3 result in lower sensitivity. We demonstrate that by exploiting destructive QI arising in the meta-substituted pyrene, it is possible to calibrate the sensor to enhance both its sensitivity and chemical selectivity by almost two orders of magnitude so that individual molecules can be detected and distinguished. These results provide a fundamental strategy to design high-performance chemical sensors with graphene functional blocks.

Destructive quantum interference (QI) has been a source of interest as a new paradigm for molecular electronics as the electronic conductance is widely dependent on the occurrence or absence of destructive QI effects. In order to interpret experimentally observed transmission features, it is necessary to understand the effects of all components of the junction on electron transport. We perform non-equilibrium Green's function calculations within the framework of density functional theory to assess the structure-function relationship of transport through pyrene molecular junctions with distinct QI properties. The chemical nature of the anchor groups and the electrodes controls the Fermi level alignment, which determines the observability of destructive QI. A thorough analysis allows to disentangle the transmission features arising from the molecule and the electrodes. Interestingly, graphene electrodes introduce features in the low-bias regime, which can either mask or be misinterpreted as QI effects, while instead originating from the topological properties of the edges. Thus, this first principles analysis provides clear indications to guide the interpretation of experimental studies, which cannot be obtained from simple Hückel model calculations.

The ability to control the charge state of individual molecules wired in two-terminal single-molecule junctions is a key challenge in molecular electronics, particularly in relation to the development of molecular memory and other computational componentry. Here we demonstrate that single porphyrin molecular junctions can be reversibly charged and discharged at elevated biases under ambient conditions due to the presence of a localised molecular eigenstate close to the Fermi edge of the electrodes. In particular, we can observe long-lived charge-states with lifetimes upwards of 1-10 seconds after returning to low bias and large changes in conductance, in excess of 100-fold at low bias. Our theoretical analysis finds charge-state lifetimes within the same time range as the experiments. The ambient operation demonstrates that special conditions such as low temperatures or ultra-high vacuum are not essential to observe hysteresis and stable charged molecular junctions.

Localized basis sets in the projector augmented wave formalism allow for computationally efficient calculations within density functional theory (DFT). However, achieving high numerical accuracy requires an extensive basis set, which also poses a fundamental problem for the interpretation of the results. We present a way to obtain a reduced basis set of atomic orbitals through the subdiagonalization of each atomic block of the Hamiltonian. The resulting local orbitals (LOs) inherit the information of the local crystal field. In the LO basis, it becomes apparent that the Hamiltonian is nearly block-diagonal, and we demonstrate that it is possible to keep only a subset of relevant LOs that provide an accurate description of the physics around the Fermi level. This reduces to some extent the redundancy of the original basis set, and at the same time, it allows one to perform post-processing of DFT calculations, ranging from the interpretation of electron transport to extracting effective tight-binding Hamiltonians, very efficiently and without sacrificing the accuracy of the results.

In our theoretical study where we combine a nonequilibrium Green's function (NEGF) approach with density functional theory (DFT) we investigate compounds containing a ferrocene moiety which is connected to i) thiol anchor groups on both sides in a para-connection, ii) a pyridyl anchor group on one side in a meta-connection and a thiol group on the other side in a para-connection, in both cases with acetylenic spacers in between the Ferrocene and the anchors. We predict possible single molecule junction geometries within a scanning tunneling microscopy (STM) setup where we find that the conductance trend for the set of conformations are intriguing in the sense that the conductance does not decrease while the junction length increases which we analyze and explain in terms of Fermi level alignment. We also find a pattern for the current-voltage (IV) curves within the linear-response regime for both molecules we study, where the conductance variation with the molecular configurations is surprisingly small.

In our theoretical study in which we combine a nonequilibrium Green's function approach with density functional theory we investigate branched compounds containing Ru or Os metal complexes in two branches, which due to their identical or different metal centers are symmetric or asymmetric. In these compounds the metal atoms are connected to pyridyl anchor groups via acetylenic and phenyl spacer groups in a metaconnection. We find there is no destructive quantum interference (DQI) feature in the transmission function near the Fermi level for the investigated molecules regardless of their symmetry, either in their neutral states or in their charged states. We map the structural characteristics of the range of molecules onto a simplified tight-binding model in order to identify the main differences between the molecules in this study and previously investigated ferrocene compounds in order to clarify the structural sources for DQI, which we found for the latter but not for the former. We also find that local charging on one of the branches changes only the conductance by about one order of magnitude, which we explain in terms of the spatial distributions and charge-induced energy shifts of the relevant molecular orbitals for the branched compounds.

In our theoretical study where we combine a nonequilibrium Green's function approach with density functional theory we investigate branched compounds containing ferrocene moieties in both branches which, due to their metal centers, are designed to allow for asymmetry induced by local charging. In these compounds the ferrocene moieties are connected to pyridyl anchor groups either directly or via acetylenic spacers in a metaconnection, where we also compare our results with those obtained for the respective single-branched molecules with both meta-and paraconnections between the metal center and the anchors. We find a destructive quantum interference (DQI) feature in the transmission function slightly below the lowest unoccupied molecular orbital, which dominates the conductance even for the uncharged branched compound with spacer groups inserted. In an analysis based on mapping the structural characteristics of the range of molecules in our article onto tight-binding models, we identify the structural source of the DQI minimum as the through-space coupling between the pyridyl anchor groups. We also find that local charging in one of the branches changes the conductance only by about one order of magnitude, which we explain in terms of the spatial distributions of the relevant molecular orbitals for the branched compounds.

Destructive quantum interference (DQI) in single molecule electronics is a purely quantum mechanical effect and is entirely defined by the inherent properties of the molecule in the junction such as its structure and symmetry. This definition of DQI by molecular properties alone suggests its relation to other more general concepts in chemistry as well as the possibility of deriving simple models for its understanding and molecular device design. Recently, two such models have gained a wide spread attention, where one was a graphical scheme based on visually inspecting the connectivity of the carbon sites in conjugated π systems in an atomic orbital (AO) basis and the other one puts the emphasis on the amplitudes and signs of the frontier molecular orbitals (MOs). There have been discussions on the range of applicability for these schemes, but ultimately conclusions from topological molecular Hamiltonians should not depend on whether they are drawn from an AO or a MO representation, as long as all the orbitals are taken into account. In this article, we clarify the relation between both models in terms of the zeroth order Green’s function and compare their predictions for a variety of systems. From this comparison, we conclude that for a correct description of DQI from a MO perspective, it is necessary to include the contributions from all MOs rather than just those from the frontier orbitals. The cases where DQI effects can be successfully predicted within a frontier orbital approximation we show them to be limited to alternant even-membered hydrocarbons, as a direct consequence of the Coulson-Rushbrooke pairing theorem in quantum chemistry.

The paper provides a comprehensive theoretical description of electron transport through transition metal complexes in single molecule junctions, where the main focus is on an analysis of the structural parameters responsible for bias-induced conductance switching as found in recent experiments, where an interpretation was provided by our simulations. The switching could be theoretically explained by a two-channel model combining coherent electron transport and electron hopping, where the underlying mechanism could be identified as a charging of the molecule in the junction made possible by the presence of a localized electronic state on the transition metal center. In this article, we present a framework for the description of an electron hopping-based switching process within the semi-classical Marcus–Hush theory, where all relevant quantities are calculated on the basis of density functional theory (DFT). Additionally, structural aspects of the junction and their respective importance for the occurrence of irreversible switching are discussed.

Azulene (Az) is a non-alternating, aromatic hydrocarbon composed of a five-membered, electron-rich and a seven-membered, electron-poor ring; an electron distribution that provides intrinsic redox activity. By varying the attachment points of the two electrode-bridging substituents to the Az center, the influence of the redox functionality on charge transport is evaluated. The conductance of the 1,3 Az derivative is at least one order of magnitude lower than those of the 2,6 Az and 4,7 Az derivatives, in agreement with density functional theory (DFT) calculations. In addition, only 1,3 Az exhibits pronounced nonlinear current-voltage characteristics with hysteresis, indicating a bias-dependent conductance switching. DFT identifies the LUMO to be nearest to the Fermi energy of the electrodes, but to be an active transport channel only in the case of the 2,6 and the 4,7 Az derivatives, whereas the 1,3 Az derivative uses the HOMO at low and the LUMO+1 at high bias. In return, the localized, weakly coupled LUMO of 1,3 Az creates a slow electron-hopping channel responsible for the voltage-induced switching due to the occupation of a single molecular orbital (MO).

AbstractAzulen (Az) ist ein nicht‐alternierender aromatischer Kohlenwasserstoff, der aus einem elektronenreichen Ring mit fünf und einem elektronenarmen Ring mit sieben Gliedern besteht. Diese Elektronenverteilung führt zu einer intrinsischen Redoxaktivitiät. Durch Variieren der Befestigungspunkte der beiden elektrodenverbindenden Substituenten am Az‐Zentrum wird der Einfluss der Redoxfunktionalität auf den Ladungstransport untersucht. Die Leitfähigkeit des 1,3‐Az‐Derivats ist dabei mindestens eine Größenordnung geringer als die der 2,6‐Az‐ und 4,7‐Az‐Derivate, in Übereinstimmung mit Resultaten der Dichtefunktionaltheorie (DFT). Außerdem weist nur 1,3‐Az eine ausgeprägt nichtlineare Strom‐Spannungs‐Kennlinie mit Hysterese auf, was durch einen spannungsabhängigen Schaltvorgang ausgelöst wird. Die DFT‐Ergebnisse zeigen weiter, dass das LUMO einerseits am nächsten bei der Fermi‐Energie der Elektroden liegt, andererseits aber nur bei den 2,6‐ und 4,7‐Az‐Derivaten den aktiven Ladungskanal darstellt, während das 1,3‐Az‐Derivat bei niedriger Spannung das HOMO und bei hoher Spannung das LUMO+1 verwendet. Im Gegenzug schafft das lokalisierte, nur schwach gekoppelte LUMO des 1,3‐Az einen zusätzlichen, langsamen Transportkanal für Elektronenhopping, das für den spannungsinduzierte Schaltvorgang aufgrund der Besetzung eines einzelnen Molekülorbitals (MO) verantwortlich ist.

For defining the conductance of single molecule junctions with a redox functionality in an electrochemical cell, two conceptually different electron transport mechanisms, namely coherent tunnelling and vibrationally induced hopping compete with each other, where implicit parameters of the setup such as the length of the molecule and the applied gate voltage decide which mechanism is the dominant one. Although coherent tunnelling is most efficiently described within Landauer theory, while the common theoretical treatment of electron hopping is based on Marcus theory, both theories are adequate for the processes they describe without introducing accuracy limiting approximations. For a direct comparison, however, it has to be ensured that the crucial quantities obtained from electronic structure calculations, i.e. the transmission function T(E) in Landauer theory, and the transfer integral V, the reorganisation energy $\lambda$ and the driving force $\Delta G^0$ in Marcus theory, are derived from similar grounds as pointed out by Nitzan and co-workers in a series of publications. In this article our framework is a single particle picture, where we perform density functional theory calculations for the conductance corresponding to both transport mechanisms for junctions with the central molecule containing one, two or three Ruthenium centers, respectively, where we extrapolate our results in order to define the critical length for the transition point of the two regimes which we identify at 5.76 nm for this type of molecular wire. We also discuss trends in dependence on an electrochemically induced gate potential.

Charge transport through single molecules can be influenced by the charge and spin states of redox-active metal centres placed in the transport pathway. These intrinsic properties are usually manipulated by varying the molecule's electrochemical and magnetic environment, a procedure that requires complex setups with multiple terminals. Here we show that oxidation and reduction of organometallic compounds containing either Fe, Ru or Mo centres can solely be triggered by the electric field applied to a two-terminal molecular junction. Whereas all compounds exhibit bias-dependent hysteresis, the Mo-containing compound additionally shows an abrupt voltage-induced conductance switching, yielding high-to-low current ratios exceeding 1,000 at bias voltages of less than 1.0 V. Density functional theory calculations identify a localized, redox-active molecular orbital that is weakly coupled to the electrodes and closely aligned with the Fermi energy of the leads because of the spin-polarized ground state unique to the Mo centre. This situation provides an additional slow and incoherent hopping channel for transport, triggering a transient charging effect in the entire molecule with a strong hysteresis and large high-to-low current ratios.

Nonequilibrium Green's function techniques (NEGF) combined with density functional theory (DFT) calculations have become a standard tool for the description of electron transport through single molecule nanojunctions in the coherent tunneling (CT) regime. However, the applicability of these methods for transport in the Coulomb blockade (CB) regime is questionable. For a molecular assembly model, with multideterminant calculations as a benchmark, we show how a closed-shell ansatz, the usual ingredient of mean-field methods, fails to properly describe the step like electron-transfer characteristic in weakly coupled systems. Detailed analysis of this misbehavior allows us to propose a practical scheme to extract the addition energies in the CB regime for single-molecule junctions from NEGF DFT within the local-density approximation (closed shell). We show also that electrostatic screening effects are taken into account within this simple approach.

There are various quantum chemical approaches for an ab initio description of transfer integrals within the framework of Marcus theory in the context of electron transfer reactions. In our paper, we aim to calculate transfer integrals in redox-active single molecule junctions, where we focus on the coherent tunneling limit with the metal leads taking the position of donor and acceptor and the molecule acting as a transportmediating bridge. This setup allows us to derive a conductance, which can be directly compared with recent results from a nonequilibrium Green's function approach. Compared with purely molecular systems we face additional challenges due to the metallic nature of the leads, which rules out some of the common techniques, and due to their periodicity, which requires k-space integration. We present three different methods, all based on density functional theory, for calculating the transfer integral under these constraints, which we benchmark on molecular test systems from the relevant literature. We also discuss many-body effects and apply all three techniques to a junction with a Ruthenium complex in different oxidation states.