Quantum conductance of multiwall carbon nanotubes

An analysis on the conductance of multiwall carbon nanotubes (MWNT's) is presented. Recent experiment indicated that MWNT's are good quantum conductors. Our theory shows that tunneling current between states on different walls of a defect-free, infinitely long MWNT is vanishingly small in general, which leads to the quantization of the conductance of the MWNT's. With a reasonable simple model, we explicitly show that the conductance of a capped MWNT can be determined by the outermost wall for an infinitely long nanotube. We apply the theory to finite MWNT's and estimate the generic interwall conductance to be negligible compared to the intrawall ballistic conductance.

Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Monte Carlo simulation of cluster growth in surface defects induced by the tip of a scanning tunnelling microscope

First steps are taken to model the electrochemical deposition of metals in nanometer-sized cavities. In the present work, the electrochemical deposition of Cu atoms in nanometer-sized holes dug on Au(111) is investigated through Monte Carlo simulations using the embedded atom method to represent particle interactions. By sweeping the chemical potential of Cu, a cluster is allowed to grow within the hole rising four atomic layers above the surface. Its lateral extension remains confined to the area defined by the borders of the original defect. (C) 2004 Elsevier B.V. All rights reserved.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

Compton profiles of Si: Pseudopotential calculation and reconstruction effects

Using an ab initio pseudopotential calculation, we compute Compton profiles of silicon along the (100), (110), and (111) directions, and then reconstruct the pseudo-wave-functions to regain the oscillatory behavior of the all-electron valence wave functions near the atomic cores. We study the effect that this reconstruction has on the Compton profiles and their anisotropies. We find a decrease in the magnitude of the profiles at small wave vectors and in their anisotropies. These changes bring the theoretical predictions closer to experimental results.

Broken symmetry and pseudogaps in ropes of carbon nanotubes

We investigate the influence of tube-tube interactions in ropes of (10,10) carbon nanotubes, and find that these effects induce a pseudogap in the density of state (DOS) of the rope of width 0.1 eV at the Fermi level. In an isolated (n,n) carbon nanotube there are two bands that cross in a linear fashion at the Fermi level, making the nanotube metallic with a DOS that is constant in a 1.5 eV wide window around the Fermi energy. The presence of the neighbouring tubes causes these two bands to repel, opening up a band gap that can be as large as 0.3 eV. The small dispersion in the plane perpendicular to the rope smears out this gap for a rope with a large cross-sectional area, and we see a pseudogap at the Fermi energy in the DOS where the DOS falls to one third of its value for the isolated tube. This phenomenon should affect many properties of the behavior of ropes of (n,n) nanotubes, which should display a more semimetallic character than expected in transport and doping experiments, with the existence of both hole and electron carriers leading to qualitatively different thermopower and Hall-effect behaviors from those expected for a normal metal. Band repulsion like this can be expected to occur for any tube perturbed by a sufficiently strong interaction, for example, from contact with a surface or with other tubes.

On the generation of metal clusters with the electrochemical scanning tunneling microscope

In the present work we consider two aspects of the deposition of metal clusters on an electrode surface. The formation of such clusters with the tip of a scanning tunneling microscope is simulated by atom dynamics. Subsequently the stability of these clusters is investigated by Monte Carlo simulations in a grand-canonical ensemble. In particular, the following systems were considered explicitly: Pd clusters on Au(111), Cu on Au(111), Ag on Au(111), Pb on Au(111) and Cu on Ag(111). The analysis of the results obtained for the different systems leads to the conclusion that optimal systems for nanostructuring are those where the metals participating have similar cohesive energies and negative heats of alloy formation. In this respect, the system Cu-Pd(111) is predicted as a good candidate for the formation of stable clusters. (c) 2005 Elsevier B.V. All rights reserved.

Current-driven magnetic rearrangements in spin-polarized point contacts

A method for investigating the dynamics of atomic magnetic moments in current-carrying magnetic point contacts under bias is presented. This combines the nonequilibrium Green's function (NEGF) method for evaluating the current and the charge density with a description of the dynamics of the magnetization in terms of quasistatic thermally activated transitions between stationary configurations. This method is then implemented in a tight-binding (TB) model with parameters chosen to simulate the main features of the electronic structures of magnetic transition metals. We investigate the domain wall (DW) migration in magnetic monoatomic chains sandwiched between magnetic leads, and for realistic parameters find that collinear arrangement of the magnetic moments of the chain is always favorable. Several stationary magnetic configurations are identified, corresponding to a different number of Bloch walls in the chain and to a different current. The relative stability of these configurations depends on the geometrical details of the junction and on the bias; however, we predict transitions between different configurations with activation barriers of the order of a few tens of meV. Since different magnetic configurations are associated with different resistances, this suggests an intrinsic random telegraph noise at microwave frequencies in the I-V curves of magnetic atomic point contacts at room temperature. Finally, we investigate whether or not current-induced torques are conservative.

Correlated electron-ion dynamics with open boundaries: formalism

We extend a new formalism, which allows correlated electron-ion dynamics to be applied to the problem of open boundary conditions. We implement this at the first moment level (allowing heating of ions by electrons) and observe the expected cooling in the classical part of the ionic kinetic energy and current-induced heating in the quantum contribution. The formalism for open boundaries should be easily extended to higher moments of the correlated electron-ion fluctuations.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Transport in nanoscale systems: the microcanonical versus grand-canonical picture

We analyse a picture of transport in which two large but finite charged electrodes discharge across a nanoscale junction. We identify a functional whose minimization, within the space of all bound many-body wavefunctions, defines an instantaneous steady state. We also discuss factors that favour the onset of steady-state conduction in such systems, make a connection with the notion of entropy, and suggest a novel source of steady-state noise. Finally, we prove that the true many-body total current in this closed system is given exactly by the one-electron total current, obtained from time-dependent density-functional theory.

Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

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.