Power dissipation in nanoscale conductors: classical, semi-classical and quantum dynamics

Modelling Joule heating is a difficult problem because of the need to introduce correct correlations between the motions of the ions and the electrons. In this paper we analyse three different models of current induced heating (a purely classical model, a fully quantum model and a hybrid model in which the electrons are treated quantum mechanically and the atoms are treated classically). We find that all three models allow for both heating and cooling processes in the presence of a current, and furthermore the purely classical and purely quantum models show remarkable agreement in the limit of high biases. However, the hybrid model in the Ehrenfest approximation tends to suppress heating. Analysis of the equations of motion reveals that this is a consequence of two things: the electrons are being treated as a continuous fluid and the atoms cannot undergo quantum fluctuations. A means for correcting this is suggested.

Electron-phonon interaction in atomic-scale conductors: Einstein oscillators versus full phonon modes

Two extreme pictures of electron-phonon interactions in nanoscale conductors are compared: one in which the vibrations are treated as independent Einstein atomic oscillators, and one in which electrons are allowed to couple to the full, extended phonon modes of the conductor. It is shown that, under a broad range of conditions, the full-mode picture and the Einstein picture produce essentially the same net power at any given atom in the nanojunction. The two pictures begin to differ significantly in the limit of low lattice temperature and low applied voltages, where electron-phonon scattering is controlled by the detailed phonon energy spectrum. As an illustration of the behaviour in this limit, we study the competition between trapped vibrational modes and extended modes in shaping the inelastic current-voltage characteristics of one-dimensional atomic wires.

Inelastic current-voltage spectroscopy of atomic wires

A tight-binding model is developed to describe the electron-phonon coupling in atomic wires under an applied voltage and to model, their inelastic current-voltage spectroscopy. Particular longitudinal phonons are found to have greatly enhanced coupling to the electronic states of the system. This leads to a large drop in differential conductance at threshold energies associated with these phonons. It is found that with increasing tension these energies decrease, while the size of the conductance drops increases, in agreement with experiment.

Power dissipation in nanoscale conductors

A previous tight-binding model of power dissipation in a nanoscale conductor under an applied bias is extended to take account of the local atomic topology and the local electronic structure. The method is used to calculate the power dissipated at every atom in model nanoconductor geometries: a nanoscale constriction, a one-dimensional atomic chain between two electrodes with a resonant double barrier, and an irregular nanowire with sharp corners. The local power is compared with the local current density and the local density of states. A simple relation is found between the local power and the current density in quasiballistic geometries. A large enhancement in the power at special atoms is found in cases of resonant and anti-resonant transmission. Such systems may be expected to be particularly unstable against current-induced modifications.

Tight-binding simulation of current-carrying nanostructures

The tight-binding (TB) approach to the modelling of electrical conduction in small structures is introduced. Different equivalent forms of the TB expression for the electrical current in a nanoscale junction are derived. The use of the formalism to calculate the current density and local potential is illustrated by model examples. A first-principles time-dependent TB formalism for calculating current-induced forces and the dynamical response of atoms is presented. An earlier expression for current-induced forces under steady-state conditions is generalized beyond local charge neutrality and beyond orthogonal TB. Future directions in the modelling of power dissipation and local heating in nanoscale conductors are discussed.

Counterbalancing forces in electromigration

In electromigration (EM) experiments on metallic wires, a flux of atoms can lead to motion of the centre of mass (COM) of the wire. Hence, it may be tempting to assume that the flow of current produces a net force on the wire as a whole. We point out, on the basis of known momentum-balance arguments, that the net force on a metallic wire due to a passing steady-state current is zero. This is possible, because in addition to EM driving forces, acting on scattering centres, there are counterbalancing forces, acting on the rest of the system. Drift of the COM in EM experiments occurs inevitably because the substrate keeps the crystal lattice of the wire fixed, while allowing diffusion of defects in the bulk of the wire. This drift is not evidence for a net force on the wire.

Time-dependent tight binding

Starting from a Lagrangian mean-field theory, a set of time-dependent tight-binding equations is derived to describe dynamically and self-consistently an interacting system of quantum electrons and classical nuclei. These equations conserve norm, total energy and total momentum. A comparison with other tight-binding models is made. A previous tight-binding result for forces on atoms in the presence of electrical current flow is generalized to the time-dependent domain and is taken beyond the limit of local charge neutrality.

Non-linear conductance of disordered quantum wires

The self-consistent electron potential in a current-carrying disordered quantum wire is spatially inhomogeneous due to the formation of resistivity dipoles across scattering centres. In this paper it is argued that these inhomogeneities in the potential result in a suppression of the differential conductance of such a wire at finite applied voltage. A semi-classical argument allows this suppression, quadratic in the voltage, to be related directly to the amount of intrinsic defect scattering in the wire. This result is then tested against numerical calculations.

The transfer of energy between electrons and ions in solids

In this review we consider those processes in condensed matter that involve the irreversible flow of energy between electrons and nuclei that follows from a system being taken out of equilibrium. We survey some of the more important experimental phenomena associated with these processes, followed by a number of theoretical techniques for studying them. The techniques considered are those that can be applied to systems containing many nonequivalent atoms. They include both perturbative approaches (Fermi's Golden Rule and non-equilibrium Green's functions) and molecular dynamics based (the Ehrenfest approximation, surface hopping, semi-classical Gaussian wavefunction methods and correlated electron-ion dynamics). These methods are described and characterized, with indications of their relative merits.

Effects of plastic coating on K-alpha yield from ultra-short pulse

A study of the K-alpha radiation emitted from Ti foils irradiated with intense, similar to0.2 J, 67 fs, 800 nm laser pulses, scanning a range of intensities (similar to10(15)-10(18) W cm(-2)), is reported. The brightness of single-shot K-alpha line emission from the front of the targets is recorded. The yield from bare titanium (Ti) is compared to that from plastic (parylene-E) coated Ti. It is demonstrated that, for a defocused beam, a thin layer of plastic increases the yield.

Time dependence of rotational state populations of excited hydrogen molecules in an RF excited plasma reactor

We report on time-dependent population distributions of excited rotational states of hydrogen in a capacitively coupled RF discharge. The common model to obtain the gas temperature from the rotational distribution is not applicable at all times during the discharge cycle due to the time dependence of the EEDF. The apparent temperature within a cycle assumes values between 350 K and 450 K for the discharge parameters of this experiment. We discuss the optimum time window within the discharge cycle that yields the best approximation to the actual temperature. Erroneous results can be obtained, in principle, with time-integrated measurements; we find, however, that in the present case the systematic error amounts to only approximately 20 K. This is due to the fact that the dominant contribution to the average intensity arises during that time window for which the assumptions underlying the analysis are best fulfilled. A similar analysis can be performed for N+2 rotational bands with a small amount of nitrogen added to the discharge gas. These populations do not exhibit the time variations found in the case of H2.