Plato: a localised orbital based density functional theory code

The Plato package allows both orthogonal and non-orthogonal tight-binding as well as density functional theory (DFT) calculations to be performed within a single framework. The package also provides extensive tools for analysing the results of simulations as well as a number of tools for creating input files. The code is based upon the ideas first discussed in Sankey and Niklewski (1989) [1] with extensions to allow high-quality DFT calculations to be performed. DFT calculations can utilise either the local density approximation or the generalised gradient approximation. Basis sets from minimal basis through to ones containing multiple radial functions per angular momenta and polarisation functions can be used. Illustrations of how the package has been employed are given along with instructions for its utilisation

Density-functional study of adsorption of Co on Si(100)

We have studied the stable sites for Co both on the surface of Si(100) and subsurface by using ab initio methods. We show that the most stable surface site for Co is situated in the dimer trenches (the low site). The subsurface sites that we study are all found to be more stable than the most stable surface site. The most stable subsurface site is the under-dimer site. The most stable site of all is, however, the dimer vacancy site formed by removing the dimer above the cobalt in the under-dimer site

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

Inelastic quantum transport: the self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared with a steady-state method based on non-equilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born Approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations, where differences begin to emerge.Comment: 12 pages, 7 figure

Analog of Rabi oscillations in resonant electron-ion systems

Quantum coherence between electron and ion dynamics, observed in organic semiconductors by means of ultrafast spectroscopy, is the object of recent theoretical and computational studies. To simulate this kind of quantum coherent dynamics, we have introduced in a previous article [L. Stella, M. Meister, A. J. Fisher, and A. P. Horsfield, J. Chem. Phys. 127, 214104 (2007)] an improved computational scheme based on Correlated Electron-Ion Dynamics (CEID). In this article, we provide a generalization of that scheme to model several ionic degrees of freedom and many-body electronic states. To illustrate the capability of this extended CEID, we study a model system which displays the electron-ion analog of the Rabi oscillations. Finally, we discuss convergence and scaling properties of the extended CEID along with its applicability to more realistic problems. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3589165

Numerical atomic orbitals for linear scaling

The performance of basis sets made of numerical atomic orbitals is explored in density-functional calculations of solids and molecules. With the aim of optimizing basis quality while maintaining strict localization of the orbitals, as needed for linear-scaling calculations, several schemes have been tried. The best performance is obtained for the basis sets generated according to a new scheme presented here, a flexibilization of previous proposals. The basis sets are tested versus converged plane-wave calculations on a significant variety of systems, including covalent, ionic and metallic. Satisfactory convergence (deviations significantly smaller than the accuracy of the underlying theory) is obtained for reasonably small basis sizes, with a clear improvement over previous schemes. The transferability of the obtained basis sets is tested in several cases and it is found to be satisfactory as well.Comment: 9 pages with 2 encapsulated postscript figures, submitted to Phys. Rev.