The small-bias conductance of the C-6 molecule, stretched between two metallic leads, is studied using time-dependent density functional theory within the adiabatic local density approximation. The leads are modeled by jellium slabs, the electronic density and the current density are described on a grid, whereas the core electrons and the highly oscillating valence orbitals are approximated using standard norm-conserving pseudopotentials. The jellium leads are supplemented by a complex absorbing potential that serves to absorb charge reaching the edge of the electrodes and hence mimic irreversible flow into the macroscopic metal. The system is rapidly exposed to a ramp potential directed along the C-6 axis, which gives rise to the onset of charge and current oscillations. As time progresses, a fast redistribution of the molecular charge is observed, which translates into a direct current response. Accompanying the dc signal, alternating current fluctuations of charge and currents within the molecule and the metallic leads are observed. These form the complex impedance of the molecule and are especially strong at the plasmon frequency of the leads and the lowest excitation peak of C-6. We study the molecular conductance in two limits: the strong coupling limit, where the edge atoms of the chain are submerged in the jellium and the weak coupling case, where the carbon atoms and the leads do not overlap spatially. (C) 2004 American Institute of Physics.

Nanoshells have been previously shown to have tunable absorption frequencies that are dependent on the ratio of their inner and outer radii. Inspired by this, we ask: can a nanoshell increase the absorption of a small core system embedded within it? A theoretical model is constructed to answer this question. A core, composed of a “jellium” ball of the density of gold is embedded within a jellium nanoshell of nanometric diameter. The shell plasmon frequency is tuned to the core absorption line. A calculation based the time-dependent density functional theory was performed showing a 10 fold increase in core excitation yield.

We present a linear-response approach for time-dependent density-functional theories using time-adiabatic functionals. The resulting theory can be performed both in the time and in the frequency domain. The derivation considers an impulsive perturbation after which the Kohn-Sham orbitals develop in time autonomously. The equation describing the evolution is not strictly linear in the wave function representation. Only after going into a symplectic real-spinor representation does the linearity make itself explicit. For performing the numerical integration of the resulting equations, yielding the linear response in time, we develop a modified Chebyshev expansion approach. The frequency domain is easily accessible as well by changing the coefficients of the Chebyshev polynomial, yielding the expansion of a formal symplectic Green's operator. (C) 2004 American Institute of Physics.

A numerical method is given for effecting nonlinear local density functional evolution. Within a given time interval, Chebyshev quadrature points are used to sample the evolving orbitals. An implicit equation coupling wave functions at the different time points is then set up. The equation is solved iteratively using the ‘‘direct inversion in iterative space’’ acceleration technique. Spatially, the orbitals are represented on a Fourier grid combined with soft pseudopotentials. The method is first applied to the computation of the 3Pg adiabatic potential energy curves of Al2 . Next, the electronic dynamics of a toy molecular wire is studied. The wire consists of a C2H4 molecule connected via sulfur atoms to two gold atoms, the ‘‘electrodes.’’ The molecule is placed in a homogeneous electric field and a dynamical process of charge transfer is observed. By comparing the transient with that of a resistance-capacitance circuit, an effective Ohmic resistance and capacitance is estimated for the system.

The representation of a quantum system by an evenly spaced Fourier grid is examined. This grid faithfully represents wave functions whose projection is contained in a rectangular phase space. This is mathematically equivalent to a band limited function with finite support. In general, have packets decay exponentially in classically forbidden regions of phase space. This idea is then used first to optimize the rectangular shape of the Fourier grid, leading to exponential convergence. Nevertheless, in most cases the representation is suboptimal. The representation efficiency can then be extremely enhanced by mapping the coordinates. The mapping procedure reshapes the wave function to fit into the rectangular Fourier shape such that the wasted phase space area is minimal. It is shown that canonical transformations, which rescale the coordinates, improve the representation dramatically. A specific scaling transformation enables the representation of the notoriously difficult Coulomb potentials. The scaling transformation enables one to extract almost as many converged eigenstate energies as there are grid points. The method is extendible to more than one dimension, which is demonstrated by the study of the H + 2 problem. This scaling transformation can bridge the gap between quantum chemistry and quantum molecular dynamics by enabling the treatment of electronic problems in the vicinity of Coulomb potentials by grid methods developed for molecular dynamics.