Time-dependent view of sequential transport through molecules with rapidly fluctuating bridges

Coherent and incoherent contributions to molecular electron transport

We numerically isolate the limits of validity of the Landauer approximation to describe charge transport along molecular junctions in condensed phase environments. To do so, we contrast Landauer with exact time-dependent non-equilibrium Green’s function quantum transport computations in a two-site molecular junction subject to exponentially correlated noise. Under resonant transport conditions, we find Landauer accuracy to critically depend on intramolecular interactions. By contrast, under nonresonant conditions, the emergence of incoherent transport routes that go beyond Landauer depends on charging and discharging processes at the electrode–molecule interface. In both cases, decreasing the rate of charge exchange between the electrodes and molecule and increasing the interaction strength with the thermal environment cause Landauer to become less accurate. The results are interpreted from a time-dependent perspective where the noise prevents the junction from achieving steady-state and from a fully quantum perspective where the environment introduces dephasing in the dynamics. Using these results, we analyze why the Landauer approach is so useful to understand experiments, isolate regimes where it fails, and propose schemes to chemically manipulate the degree of transport coherence.

Quantum Phonon Transport in Nanomaterials: Combining Atomistic with Non-Equilibrium Green's Function Techniques

A crucial goal for increasing thermal energy harvesting will be to progress towards atomistic design strategies for smart nanodevices and nanomaterials. This requires the combination of computationally efficient atomistic methodologies with quantum transport based approaches. Here, we review our recent work on this problem, by presenting selected applications of the PHONON tool to the description of phonon transport in nanostructured materials. The PHONON tool is a module developed as part of the Density-Functional Tight-Binding (DFTB) software platform. We discuss the anisotropic phonon band structure of selected puckered two-dimensional materials, helical and horizontal doping effects in the phonon thermal conductivity of boron nitride-carbon heteronanotubes, phonon filtering in molecular junctions, and a novel computational methodology to investigate time-dependent phonon transport at the atomistic level. These examples illustrate the versatility of our implementation of phonon transport in combination with density functional-based methods to address specific nanoscale functionalities, thus potentially allowing for designing novel thermal devices.

Effect of Structure and Disorder on the Charge Transport in Defined Self-Assembled Monolayers of Organic Semiconductors

Self-assembled monolayer field-effect transistors (SAMFETs) are not only a promising type of organic electronic device but also allow detailed analyses of structure-property correlations. The influence of the morphology on the charge transport is particularly pronounced, due to the confined monolayer of 2D-π-stacked organic semiconductor molecules. The morphology, in turn, is governed by relatively weak van-der-Waals interactions and is thus prone to dynamic structural fluctuations. Accordingly, combining electronic and physical characterization and time-averaged X-ray analyses with the dynamic information available at atomic resolution from simulations allows us to characterize self-assembled monolayer (SAM) based devices in great detail. For this purpose, we have constructed transistors based on SAMs of two molecules that consist of the organic p-type semiconductor benzothieno[3,2-b][1]benzothiophene (BTBT), linked to a C₁₁ or C₁₂ alkylphosphonic acid. Both molecules form ordered SAMs; however, our experiments show that the size of the crystalline domains and the charge-transport properties vary considerably in the two systems. These findings were confirmed by molecular dynamics (MD) simulations and semiempirical molecular-orbital electronic-structure calculations, performed on snapshots from the MD simulations at different times, revealing, in atomistic detail, how the charge transport in organic semiconductors is influenced and limited by dynamic disorder.

Influence of Coherent Tunneling and Incoherent Hopping on the Charge Transfer Mechanism in Linear Donor-Bridge-Acceptor Systems

The mechanism of charge transfer has been observed to change from tunneling to hopping with increasing numbers of DNA base pairs in polynucleotides and with the length of molecular wires. The aim of this paper is to investigate this transition by examining the population dynamics using a tight-binding Hamiltonian with model parameters to describe a linear donor-bridge-acceptor (D-B-A) system. The model includes a primary vibration and an electron-vibration coupling at each site. A further coupling of the primary vibration with a secondary phonon bath allows the system to dissipate energy to the environment and reach a steady state. We apply the quantum master equation (QME) approach, based on second-order perturbation theory in a quantum dissipative system, to examine the dynamical processes involved in charge-transfer and follow the population transfer rate at the acceptor, ka, to shed light on the transition from tunneling to hopping. With a small tunneling parameter, V, the on-site population tends to localize and form polarons, and the hopping mechanism dominates the transfer process. With increasing V, the population tends to be delocalized and the tunneling mechanism dominates. The competition between incoherent hopping and coherent tunneling governs the mechanism of charge transfer. By varying V and the total number of sites, we also examine the onset of the transition from tunneling to hopping with increasing length.

Time-dependent quantum transport through an interacting quantum dot beyond sequential tunneling: second-order quantum rate equations

A general theoretical formulation for the effect of a strong on-site Coulomb interaction on the time-dependent electron transport through a quantum dot under the influence of arbitrary time-varying bias voltages and/or external fields is presented, based on slave bosons and the Keldysh nonequilibrium Green's function (GF) techniques. To avoid the difficulties of computing double-time GFs, we generalize the propagation scheme recently developed by Croy and Saalmann to combine the auxiliary-mode expansion with the celebrated Lacroix's decoupling approximation in dealing with the second-order correlated GFs and then establish a closed set of coupled equations of motion, called second-order quantum rate equations (SOQREs), for an exact description of transient dynamics of electron correlated tunneling. We verify that the stationary solution of our SOQREs is able to correctly describe the Kondo effect on a qualitative level. Moreover, a comparison with other methods, such as the second-order von Neumann approach and Hubbard-I approximation, is performed. As illustrations, we investigate the transient current behaviors in response to a step voltage pulse and a harmonic driving voltage, and linear admittance as well, in the cotunneling regime.

Optically induced transport through semiconductor-based molecular electronics

A tight binding model is used to investigate photoinduced tunneling current through a molecular bridge coupled to two semiconductor electrodes. A quantum master equation is developed within a non-Markovian theory based on second-order perturbation theory with respect to the molecule-semiconductor electrode coupling. The spectral functions are generated using a one dimensional alternating bond model, and the coupling between the molecule and the electrodes is expressed through a corresponding correlation function. Since the molecular bridge orbitals are inside the bandgap between the conduction and valence bands, charge carrier tunneling is inhibited in the dark. Subject to the dipole interaction with the laser field, virtual molecular states are generated via the absorption and emission of photons, and new tunneling channels open. Interesting phenomena arising from memory are noted. Such a phenomenon could serve as a switch.

Using the Chebychev expansion in quantum transport calculations

Irradiation by laser pulses and a fluctuating surrounding liquid environment can, for example, lead to time-dependent effects in the transport through molecular junctions. From the theoretical point of view, time-dependent theories of quantum transport are still challenging. In one of these existing transport theories, the energy-dependent coupling between molecule and leads is decomposed into Lorentzian functions. This trick has successfully been combined with quantum master approaches, hierarchical formalisms, and non-equilibrium Green's functions. The drawback of this approach is, however, its serious limitation to certain forms of the molecule-lead coupling and to higher temperatures. Tian and Chen [J. Chem. Phys. 137, 204114 (2012)] recently employed a Chebychev expansion to circumvent some of these latter problems. Here, we report on a similar approach also based on the Chebychev expansion but leading to a different set of coupled differential equations using the fact that a derivative of a zeroth-order Bessel function can again be given in terms of Bessel functions. Test calculations show the excellent numerical accuracy and stability of the presented formalism. The time span for which this Chebychev expansion scheme is valid without any restrictions on the form of the spectral density or temperature can be determined a priori.

Dynamics of the excitonic coupling in organic crystals

We show that the excitonic coupling in molecular crystals undergoes a very large fluctuation at room temperature as a result of the combined thermal motions of the nuclei. This observation dramatically affects the description of exciton transport in organic crystals and any other phenomenon (like singlet fission or exciton dissociation) that originates from an exciton in a molecular crystal or thin film. This unexpected result is due to the predominance of the short-range excitonic coupling mechanisms (exchange, overlap, and charge-transfer mediated) over the Coulombic excitonic coupling for molecules in van der Waals contact. To quantify this effect we develop a procedure to evaluate accurately the short-range excitonic coupling (via a diabatization scheme) along a molecular dynamics trajectory of the representative molecular crystals of anthracene and tetracene.

The influence of polaron formation on exciton dissociation

The influence of the competition between polaron formation and population injection on exciton dissociation.

Charge transport in desolvated DNA

The conductivity of DNA in molecular junctions is often probed experimentally under dry conditions, but it is unclear how much of the solvent remains attached to the DNA and how this impacts its structure, electronic states, and conductivity. Classical MD simulations show that DNA is unstable if the solvent is removed completely, while a micro-hydrated system with few water molecules shows similar charge transport properties as fully solvated DNA does. This surprising effect is analyzed in detail by mapping the density functional theory-based electronic structure to a tight-binding Hamiltonian, allowing for an estimate of conductivity of various DNA sequences with snapshot-averaged Landauer's approach. The characteristics of DNA charge transport turn out to be determined by the nearest hydration shell(s), and the removal of bulk solvent has little effect on the transport.

A hybrid approach to simulation of electron transfer in complex molecular systems

Electron transfer (ET) reactions in biomolecular systems represent an important class of processes at the interface of physics, chemistry and biology. The theoretical description of these reactions constitutes a huge challenge because extensive systems require a quantum-mechanical treatment and a broad range of time scales are involved. Thus, only small model systems may be investigated with the modern density functional theory techniques combined with non-adiabatic dynamics algorithms. On the other hand, model calculations based on Marcus's seminal theory describe the ET involving several assumptions that may not always be met. We review a multi-scale method that combines a non-adiabatic propagation scheme and a linear scaling quantum-chemical method with a molecular mechanics force field in such a way that an unbiased description of the dynamics of excess electron is achieved and the number of degrees of freedom is reduced effectively at the same time. ET reactions taking nanoseconds in systems with hundreds of quantum atoms can be simulated, bridging the gap between non-adiabatic ab initio simulations and model approaches such as the Marcus theory. A major recent application is hole transfer in DNA, which represents an archetypal ET reaction in a polarizable medium. Ongoing work focuses on hole transfer in proteins, peptides and organic semi-conductors.