Materials Theory

The equilibrium shape of strained InAs quantum dots grown epitaxially on a GaAs(001) substrate is derived as a function of volume. InAs surface energies are calculated within density-functional theory, and a continuum approach is applied for the elastic relaxation energies.

Using first-principles total-energy pseudopotential calculations, we have studied the properties of chains of potassium and aluminum in nanotubes. For BN tubes, there is little interaction between the metal chains and the tubes, and the conductivity of these tubes is through carriers located at the inner part of the tube. In contrast, for small radius carbon nanotubes, there are two types of interactions: charge-transfer (dominant for alkali atoms) leading to strong ionic cohesion, and hybridization (for multivalent metal atoms) resulting in a smaller cohesion. For Al-atomic chains in carbon tubes, we show that both effects contribute. New electronic properties related to these confined atomic chains of metal are analyzed.

Recent experiments performed at high pressures indicate that ruthenium can support unusually high concentrations of oxygen at the surface. To investigate the structure and stability of high coverage oxygen structures, we performed density functional theory calculations, within the generalized gradient approximation, for O adlayers on Ru(0001) from low coverage up to a full monolayer. We achieve quantitative agreement with previous low energy electron diffraction intensity analyses for the (2x2) and (2x1) phases and predict that an O adlayer with a (1x1) periodicity and coverage of 1 monolayer can form on Ru(0001), where the O adatoms occupy hcp-hollow sites.

LEED experiments show that Li adsorbed at Cu(100) surfaces at room temperature induces a (2x1) missing row substrate reconstruction while adsorption at lower temperatures, T=180 K, results in an unreconstructed Cu(100)+c(2x2)--Li overlayer structure. Substrate reconstruction has not been observed for Na nor for K adsorption. In order to study the specific reconstruction behavior of the Li adsorbate ab initio DFT calculations have been performed on Cu(100)+Ad, Ad = Li, Na, K systems at coverages Theta_Ad=0.25-0.5 with and without reconstruction. The calculations show that the (2x1) MR reconstructed surface lies energetically above the ideal (1x1) surface by 0.2 eV per unit cell. However, alkali binding is stronger in the MR geometry as compared to that of the ideal surface where the increase in bond strength becomes smaller in going from Li to Na to K. As a result, the MR reconstructed and the overlayer adsorbate systems are energetically very close for Cu(100)+Li while for Na and K the overlayer geometry is always favored.

While the small sticking coefficient for molecular hydrogen on the Si(001) surface apparently requires a large energy barrier of adsorption, no such barrier is observed in desorption experiments. We have calculated the potential-energy surface of an H_2 molecule in front of a Si(001) surface. If we relax the Si substrate, we find an optimum desorption path with a low (\lesssim 0.3 eV) adsorption energy barrier. While molecules impinging on the surface will mostly be reflected at the larger barrier of some frozen-substrate, molecules adsorbed on the surface can desorb along the low-barrier path.

Results of recent density functional theory calculations for alkali metal adsorbates on close-packed metal surfaces are discussed. Single adatoms on the (111) surface of Al and Cu are studied with the self-consistent surface Green-function method by which the pure adsorbate-substrate interaction may be analyzed. Higher coverage ordered adlayers of K on Al(111), Na on Al(111), and Na on Al(001) are treated using the <em>ab-initio</em> pseudopotential plane wave method which affords the prediction of coverage dependent stable and metastable adsorbate geometries and phase transitions of the adsorbate layers. Together, these studies give insight and understanding into current key issues in alkali metal adsorption, namely, the nature of the adsorbate-substrate bond at low coverage and the occurrence of hitherto unanticipated adsorbate geometries, and the associated electronic properties.

In the context of perovskite quantum paraelectrics, we study the effects of a quadrupolar interaction J q , in addition to the standard dipolar one J d . We concentrate here on the nonlinear dielectric response χ (3) P , as the main response function sensitive to quadrupolar (in our case antiquadrupolar) interactions. We employ a 3D quantum four-state lattice model and mean-field theory. The results show that inclusion of quadrupolar coupling of moderate strength ( J q ∼ 1 4 J d ) is clearly accompanied by a double change of sign of χ (3) P from negative to positive, near the quantum temperature T Q where the quantum paraelectric behaviour sets in. We fit our χ (3) P to recent experimental data for SrTiO 3 , where the sign change is identified close to T Q ∼37K .

The low-temperature structure of A1C60 (A=K, Rb) is an ordered array of polymerized C60 chains, with magnetic properties that suggest a non-metallic ground state. We study the paramagnetic state of this phase using first-principles electronic-structure methods, and examine the magnetic fluctuations around this state using a model Hamiltonian. The electronic and magnetic properties of even this polymerized phase remain strongly three dimensional, and the magnetic fluctuations favor an unusual three-dimensional antiferromagnetically ordered structure with a semi-metallic electronic spectrum.

An energy functional for orbital based O(N) calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based O(N) methods; it therefore makes it possible to perform O(N) calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed.

A recently proposed linear-scaling scheme for density-functional pseudopotential calculations is described in detail. The method is based on a formulation of density functional theory in which the ground state energy is determined by minimization with respect to the density matrix, subject to the condition that the eigenvalues of the latter lie in the range [0,1]. Linear-scaling behavior is achieved by requiring that the density matrix should vanish when the separation of its arguments exceeds a chosen cutoff. The limitation on the eigenvalue range is imposed by the method of Li, Nunes and Vanderbilt. The scheme is implemented by calculating all terms in the energy on a uniform real-space grid, and minimization is performed using the conjugate-gradient method. Tests on a 512-atom Si system show that the total energy converges rapidly as the range of the density matrix is increased. A discussion of the relation between the present method and other linear-scaling methods is given, and some problems that still require solution are indicated.