Molecule-lead coupling at molecular junctions: relation between the real- and state-space perspectives

Electron Dynamics in Open Quantum Systems: The Driven Liouville-von Neumann Methodology within Time-Dependent Density Functional Theory

A first-principles approach to describe electron dynamics in open quantum systems driven far from equilibrium via external time-dependent stimuli is introduced. Within this approach, the driven Liouville-von Neumann methodology is used to impose open boundary conditions on finite model systems whose dynamics is described using time-dependent density functional theory. As a proof of concept, the developed methodology is applied to simple spin-compensated model systems, including a hydrogen chain and a graphitic molecular junction. Good agreement between steady-state total currents obtained via direct propagation and those obtained from the self-consistent solution of the corresponding Sylvester equation indicates the validity of the implementation. The capability of the new computational approach to analyze, from first principles, non-equilibrium dynamics of open quantum systems in terms of temporally and spatially resolved current densities is demonstrated. Future extensions of the approach toward the description of dynamical magnetization and decoherence effects are briefly discussed.

A single molecule diode based on gold electrodes and benzene molecule: conductivity and coupling analysis

CONTEXT

In this paper, we simulate a single-molecule diode to calculate the effective coupling and investigate the conductivity, as well as the effect of the electric field on these two parameters. First, we obtain the molecule states and energies at 0 V. The next step is to calculate the electrode/molecule coupling using the obtained electrode and molecule Hamiltonian. The electrode/molecule coupling depends on distance. By increasing the distance from 5 to 5.5 angstroms, the coupling decreases from 0.004 to 0.0002 eV. After calculating the electrode/molecule coupling, which is the most significant parameter in electron transfer, the results can be used to obtain the current-voltage and conductivity curves of the device. Simulation results demonstrate that externally applied electric field to the benzene molecule (isolated molecule) can cause a reduction in the effective coupling between the Au electrode and benzene, leading to decreased current and conductivity. Additionally, the applied electric field narrows the gap between the HOMO and LUMO energy levels.

METHODS

We conducted this computational work using Gaussian 09 software and a MATLAB code, both of which are based on the density functional theory (DFT) approach and the self-consistent field (SCF) method. For DFT calculations, we employed the three-parameter Beck hybrid exchange functional (B3), hybridized with the nonlocal correlation functional developed by Lee, Yang, and Parr (LYP). All optimizations were performed with triple-zeta polarized (TZP) split-valence 6-311G basis sets. The final step involved calculating the electrode/molecule coupling using the Huckel method and integrating the site-to-state transformation with Huckel parameters and the Fermi golden rule. After this calculation, we obtained the current-voltage and conductivity curves using MATLAB software.

A Green's Function Approach for Determining Surface Induced Broadening and Shifting of Molecular Energy Levels

Upon adsorption of a molecule onto a surface, the molecular energy levels (MELs) broaden and change their alignment. This phenomenon directly affects electron transfer across the interface and is, therefore, a fundamental observable that influences electrochemical device performance. Here, we propose a rigorous parameter-free framework, built upon the theoretical construct of Green's functions, for studying the interface between a molecule and a bulk surface and its effect on MELs. The method extends beyond the usual wide-band limit approximation, and its generality allows its use with any level of electronic structure theory. We demonstrate its ability to predict the broadening and shifting of MELs as a function of intramolecular coupling, molecule/surface coupling, and the surface density of states for a molecule with two MELs adsorbed on a one-dimensional model metal surface. The new approach could help provide guidelines for the design and experimental characterization of electrochemical devices with optimal electron transport.

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.

Interface-Induced Conservation of Momentum Leads to Chiral-Induced Spin Selectivity

We study the nonequilibrium dynamics of electron transmission from a straight waveguide to a helix with spin-orbit coupling. Transmission is found to be spin-selective and can lead to large spin polarizations of the itinerant electrons. The degree of spin selectivity depends on the width of the interface region, and no polarization is found for single-point couplings. We show that this is due to momentum conservation conditions arising from extended interfaces. We therefore identify interface structure and conservation of momentum as crucial ingredients for chiral-induced spin selectivity, and we confirm that this mechanism is robust against static disorder.

Performance of reservoir discretizations in quantum transport simulations

Quantum transport simulations often use explicit, yet finite, electronic reservoirs. These should converge to the correct continuum limit, albeit with a trade-off between discretization and computational cost. Here, we study this interplay for extended reservoir simulations, where relaxation maintains a bias or temperature drop across the system. Our analysis begins in the non-interacting limit, where we parameterize different discretizations to compare them on an even footing. For many-body systems, we develop a method to estimate the relaxation that best approximates the continuum by controlling virtual transitions in Kramers turnover for the current. While some discretizations are more efficient for calculating currents, there is little benefit with regard to the overall state of the system. Any gains become marginal for many-body, tensor network simulations, where the relative performance of discretizations varies when sweeping other numerical controls. These results indicate that typical reservoir discretizations have little impact on numerical costs for certain computational tools. The choice of a relaxation parameter is nonetheless crucial, and the method we develop provides a reliable estimate of the optimal relaxation for finite reservoirs.

A projection-based reduced-order method for electron transport problems with long-range interactions

Long-range interactions play a central role in electron transport. At the same time, they present a challenge for direct computer simulations since sufficiently large portions of the bath have to be included in the computation to accurately compute the Coulomb potential. This article presents a reduced-order approach by deriving an open quantum model for the reduced density matrix. To treat the transient dynamics, the problem is placed in a reduced-order framework. The dynamics described by the Liouville-von Neumann equation is projected to subspaces using a Petrov-Galerkin projection. In order to recover the global electron density profile as a vehicle to compute the Coulomb potential, we propose a domain decomposition approach, where the computational domain also includes segments of the bath that are selected using logarithmic grids. This approach leads to a multi-component self-energy that enters the effective Hamiltonian. We demonstrate the accuracy of the reduced model using a molecular junction built from lithium chains.

Analytic expressions for the steady-state current with finite extended reservoirs

Open-system simulations of quantum transport provide a platform for the study of true steady states, Floquet states, and the role of temperature, time dynamics, and fluctuations, among other physical processes. They are rapidly gaining traction, especially techniques that revolve around "extended reservoirs," a collection of a finite number of degrees of freedom with relaxation that maintains a bias or temperature gradient, and have appeared under various guises (e.g., the extended or mesoscopic reservoir, auxiliary master equation, and driven Liouville-von Neumann approaches). Yet, there are still a number of open questions regarding the behavior and convergence of these techniques. Here, we derive general analytical solutions, and associated asymptotic analyses, for the steady-state current driven by finite reservoirs with proportional coupling to the system/junction. In doing so, we present a simplified and unified derivation of the non-interacting and many-body steady-state currents through arbitrary junctions, including outside of proportional coupling. We conjecture that the analytic solution for proportional coupling is the most general of its form for isomodal relaxation (i.e., relaxing proportional coupling will remove the ability to find compact, general analytical expressions for finite reservoirs). These results should be of broad utility in diagnosing the behavior and implementation of extended reservoir and related approaches, including the convergence to the Landauer limit (for non-interacting systems) and the Meir-Wingreen formula (for many-body systems).

Quantum transport with electronic relaxation in electrodes: Landauer-type formulas derived from the driven Liouville-von Neumann approach

We derive the exact steady-state solutions for the simplest model systems of resonant tunneling and tunneling with destructive quantum interference from the driven Liouville-von Neumann (DLvN) approach. Under the finite-state lead condition (the two electrodes have finite states), we analyze the asymptotic behavior of the steady-state current in the two limits of electronic relaxation. Under the infinite-state lead condition, the steady-state solutions of the two model systems can be cast as Landauer-type current formulas. According to the formulas, we show that the transmission functions near the resonant peak and the antiresonant dip can be significantly influenced by electronic relaxation in the electrodes. Moreover, under intermediate and strong electronic relaxation conditions, we analytically show that the steady-state current of the DLvN approach dramatically deviates from the Landauer current when destructive quantum interference occurs. In the regime of zero electronic relaxation, our results are reduced to the Landauer formula, indicating that the DLvN approach is equivalent to the Landauer approach when the leads have infinite states without any electronic relaxation.

Open System Tensor Networks and Kramers' Crossover for Quantum Transport

Tensor networks are a powerful tool for many-body ground states with limited entanglement. These methods can nonetheless fail for certain time-dependent processes-such as quantum transport or quenches-where entanglement growth is linear in time. Matrix-product-state decompositions of the resulting out-of-equilibrium states require a bond dimension that grows exponentially, imposing a hard limit on simulation timescales. However, in the case of transport, if the reservoir modes of a closed system are arranged according to their scattering structure, the entanglement growth can be made logarithmic. Here, we apply this ansatz to open systems via extended reservoirs that have explicit relaxation. This enables transport calculations that can access steady states, time dynamics and noise, and periodic driving (e.g., Floquet states). We demonstrate the approach by calculating the transport characteristics of an open, interacting system. These results open a path to scalable and numerically systematic many-body transport calculations with tensor networks.

Imaging Time-Dependent Electronic Currents through a Graphene-Based Nanojunction

To assist the design of efficient molecular junctions, a precise understanding of the charge transport mechanisms through nanoscaled devices is of prime importance. In the present contribution, we present time- and space-resolved electron transport simulations through a nanojunction under time-dependent potential biases. We use the driven Liouville-von Neumann approach to simulate the time evolution of the one-electron density matrix under nonequilibrium conditions, which allows us to capture the ultrafast scattering dynamics, the electronic relaxation process, and the quasi-stationary current limit from the same simulation. Using local projection techniques, we map the coherent electronic current density, unraveling insightful mechanistic details of the transport on time scales ranging from atto- to picoseconds. Memory effects dominate the early time transport process, and they reveal different current patterns on short time scales in comparison to those in the long-time regime. For nanotransistors with high switching rates, the scattering perspective on electron transport should thus be favored.

Communication: Master equations for electron transport: The limits of the Markovian limit

Master equations are increasingly popular for the simulation of time-dependent electronic transport in nanoscale devices. Several recent Markovian approaches use "extended reservoirs"-explicit degrees of freedom associated with the electrodes-distinguishing them from many previous classes of master equations. Starting from a Lindblad equation, we develop a common foundation for these approaches. Due to the incorporation of explicit electrode states, these methods do not require a large bias or even "true Markovianity" of the reservoirs. Nonetheless, their predictions are only physically relevant when the Markovian relaxation is weaker than the thermal broadening and when the extended reservoirs are "sufficiently large," in a sense that we quantify. These considerations hold despite complete positivity and respect for Pauli exclusion at any relaxation strength.

Electron transport in real time from first-principles

While the vast majority of calculations reported on molecular conductance have been based on the static non-equilibrium Green's function formalism combined with density functional theory (DFT), in recent years a few time-dependent approaches to transport have started to emerge. Among these, the driven Liouville-von Neumann equation [C. G. Sánchez et al., J. Chem. Phys. 124, 214708 (2006)] is a simple and appealing route relying on a tunable rate parameter, which has been explored in the context of semi-empirical methods. In the present study, we adapt this formulation to a density functional theory framework and analyze its performance. In particular, it is implemented in an efficient all-electron DFT code with Gaussian basis functions, suitable for quantum-dynamics simulations of large molecular systems. At variance with the case of the tight-binding calculations reported in the literature, we find that now the initial perturbation to drive the system out of equilibrium plays a fundamental role in the stability of the electron dynamics. The equation of motion used in previous tight-binding implementations with massive electrodes has to be modified to produce a stable and unidirectional current during time propagation in time-dependent DFT simulations using much smaller leads. Moreover, we propose a procedure to get rid of the dependence of the current-voltage curves on the rate parameter. This method is employed to obtain the current-voltage characteristic of saturated and unsaturated hydrocarbons of different lengths, with very promising prospects.

Driven Liouville von Neumann Equation in Lindblad Form

The Driven Liouville von Neumann approach [J. Chem. Theory Comput. 2014, 10, 2927-2941] is a computationally efficient simulation method for modeling electron dynamics in molecular electronics junctions. Previous numerical simulations have shown that the method can reproduce the exact single-particle dynamics while avoiding density matrix positivity violation found in previous implementations. In this study we prove that in the limit of infinite lead models the underlying equation of motion can be cast in Lindblad form. This provides a formal justification for the numerically observed density matrix positivity conservation.