Lagrange multiplier based transport theory for quantum wires

We discuss how a Lagrange multiplier method of non-equilibrium steady state statistical mechanics can be applied to describe the electronic transport in a quantum wire. We describe a theoretical scheme using tight-binding model. The Hamiltonian of the wire is extended via a Lagrange multiplier to ``open'' the quantum system and to drive the current through it. Diagonalization of the extended Hamiltonian yields transport properties of the wire. We show that the Lagrange multiplier method is equivalent to the Landauer approach within the considered model

Schroedinger equation for current carrying states

Schr\"odinger equation with given, {\it a priori} known current is formulated. A non-zero current density is maintained in the quantum system via a subsidiary condition imposed by vector, local Lagrange multiplier. Constrained minimization of the total energy on the manifold of an arbitrary current density topology results into a non-linear self-consistent Schr\"odinger equation. The applications to electronic transport in two-terminal molecular devices are developed and new macroscopic definition of a molecular current-voltage characteristic is proposed. The Landauer formula for the conductance of an ideal one-dimensional lead is obtained within the approach. The method is examined by modeling of current carrying states of one-dimensional harmonic oscillator

Nature of well-defined conductance of amine anchored molecular junctions

Amine terminated molecules show well behaved conductance in the scanning tunneling microscope break-junction experimental measurements. We performed density functional theory based electron transport calculations to explain the nature of this phenomenon. We find that amines can be adsorbed only on apex Au atom, while thiolate group can be attached equally well to undercoordinated and clean Au surfaces. Our calculations show that only one adsorption geo metry is sterically and energetically possible for amine anchored junction whereas three different adsorption geometries with very distinct transport pro perties are almost equally probable for thiolate anchored junction. We calculated the conductance as a function of the junction stretching when the molecules are pulled by the scanning tunneling microscope tip from the Au electrode. Our calculations show that the stretching of the thiolate anchored junction during its formation is accompanied by significant electrode geometry distortio n. The amine anchored junctions exhibit very different behavior -- the electrode remains intact when the scan ning tunneling microscope tip stretches the junction

Calculations of canonical averages from the grand canonical ensemble

Grand canonical and canonical ensembles become equivalent in the thermodynamic limit, but when the system size is finite the results obtained in the two ensembles deviate from each other. In many important cases, the canonical ensemble provides an appropriate physical description but it is often much easier to perform the calculations in the corresponding grand canonical ensemble. We present a method to compute averages in canonical ensemble based on calculations of the expectation values in grand canonical ensemble. The number of particles, which is fixed in the canonical ensemble, is not necessarily the same as the average number of particles in the grand canonical ensemble

Stability analysis of multiple nonequilibrium fixed points in self-consistent electron transport calculations

We present a method to perform stability analysis of nonequilibrium fixed points appearing in self-consistent electron transport calculations. The nonequilibrium fixed points are given by the self-consistent solution of stationary, nonlinear kinetic equation for single-particle density matrix. We obtain the stability matrix by linearizing the kinetic equation around the fixed points and analyze the real part of its spectrum to assess the asymptotic time behavior of the fixed points. We derive expressions for the stability matrices within Hartree-Fock and linear response adiabatic time-dependent density functional theory. The stability analysis of multiple fixed points is performed within the nonequilibrium Hartree-Fock approximation for the electron transport through a molecule with a spin-degenerate single level with local Coulomb interaction

Calculation of semiclassical free energy differences along non-equilibrium classical trajectories

We have derived several relations, which allow the evaluation of the system free energy changes in the leading order in $\hbar^{2}$ along classically generated trajectories. The results are formulated in terms of purely classical Hamiltonians and trajectories, so that semiclassical partition functions can be computed, e.g., via classical molecular dynamics simulations. The Hamiltonians, however, contain additional potential-energy terms, which are proportional to $\hbar^{2}$ and are temperature-dependent. We discussed the influence of quantum interference on the nonequilibrium work and problems with unambiguous definition of the semiclassical work operator

Kohn-Sham equations for nanowires with direct current

The paper describes the derivation of the Kohn-Sham equations for a nanowire with direct current. A value of the electron current enters the problem as an input via a subsidiary condition imposed by pointwise Lagrange multiplier. Using the constrained minimization of the Hohenberg-Kohn energy functional, we derive a set of self-consistent equations for current carrying orbitals of the molecular wire