Materials Theory

The ground-state properties of Fe, Co, and Ni are studied with the linear-augmented-plane-wave (LAPW) method and norm-conserving pseudopotentials. The calculated lattice constant, bulk modulus, and magnetic moment with both the local-spin-density approximation (LSDA) and the generalized gradient approximation (GGA) are in good agreement with those of all-electron calculations, respectively. The GGA results show a substantial improvement over the LSDA results, i.e., better agreement with experiment. The accurate treatment of the nonlinear core-valence exchange and correlation interaction is found to be essential for the determination of the magnetic properties of 3d transition metals. The present study demonstrates the successful application of the LAPW pseudopotential approach to the calculation of ground-state properties of magnetic 3d transition metals.

Liquid NaSn alloys in five different compositions (20, 40, 50, 57 and 80% sodium) are studied using density functional calculations combined with molecular dynamics(Car-Parrinello method). The frequency-dependent electric conductivities for the systems are calculated by means of the Kubo-Greenwood formula. The extrapolated DC conductivities are in good agreement with the experimental data and reproduce the strong variation with the concentration. The maximum of conductivity is obtained, in agreement with experiment, near the equimolar composition. The strong variation of conductivity, ranging from almost semiconducting up to metallic behaviour, can be understood by an analysis of the densities-of-states.

An effective Hamiltonian for the ferroelectric transition in PbTiO3 is constructed from first-principles density-functional-theory total-energy and linear-response calculations through the use of a localized, symmetrized basis set of ``lattice Wannier functions.'' Explicit parametrization of the polar lattice Wannier functions is used for subspace projection, addressing the issues of LO-TO splitting and coupling to the complementary subspace. In contrast with ferroelectric BaTiO3 and KNbO3, we find significant involvement of the Pb atom in the lattice instability. Monte Carlo simulations for this Hamiltonian show a first-order cubic-tetragonal transition at 660 K. Resulting temperature dependence of spontaneous polarization, c/a ratio and unit-cell volume near the transition are in good agreement with experiment. Comparison of Monte Carlo results with mean field theory analysis shows that both strain and fluctuations are necessary to produce the first order character of this transition.

Density-functional simulations are used to calculate structural properties and high-symmetry phonons of the hypothetical cubic phase, the stable orthorhombic phase and an intermediate tetragonal phase of magnesium silicate perovskite. We show that the structure of the stable phase is well described by freezing in a small number of unstable phonons into the cubic phase. We use the frequencies of these unstable modes to estimate transition temperatures for cubic--tetragonal and tetragonal--orthorhombic phase transitions. These are investigated further to find that the coupling with the strain suggests that phonons give a better representation than rigid unit modes. The phonons of an intermediate tetragonal phase were found to be stable except for two rotational modes. The eigenvectors of the most unstable mode of each of the cubic and tetragonal phases account for all the positional parameters of the orthorhombic phase. The phase boundary for the orthorhombic--tetragonal transition intersects possible mantle geotherms, suggesting that the tetragonal phase may be present in the lower mantle.

Motivated by the negative thermal expansion observed for silicon between 20 K and 120 K, we present first an ab initio study of the volume dependence of interatomic force constants, phonon frequencies of TA(X) and TA(L) modes, and of the associated mode Gruneisen parameters. The influence of successive nearest neighbors shells is analysed. Analytical formulas, taking into account interactions up to second nearest neighbors, are developped for phonon frequencies of TA(X) and TA(L) modes and the corresponding mode Gruneisen parameters. We also analyze the volume and pressure dependence of various thermodynamic properties (specific heat, bulk modulus, thermal expansion), and point out the effect of the negative mode Gruneisen parameters of the acoustic branches on these properties. Finally, we present the evolution of the mean square atomic displacement and of the atomic temperature factor with the temperature for different volumes, for which the anomalous effects are even greater.

Using density-functional theory we investigate several properties of Al(111), Al(100), Al(110), and stepped Al(111) surfaces. We report results of formation energies of surfaces, steps, adatoms, and vacancies. For the adsorption and diffusion of Al on flat regions of Al(111) surfaces we find the hcp site energetically slightly preferred over the fcc site. The energy barrier for self-diffusion on Al(111) is very low (0.04eV). Coming close to one of the two sorts of close packed, monoatomic steps on Al(111), labeled according to their {111} and {100} micro-facets, Al adatoms experience an attraction of <~ 0.1eV already before direct contact with the edge of the step. This attraction has a range of several atomic spacings and is of electronic origin. Upon arrival at the lower step edge, the adatom attaches with no barrier at a low energy five-fold coordinated site. Coming from the upper terrace, it incorporates into the step by an atomic exchange process, which has a barrier below 0.1eV for both sorts of close packed steps. The barrier for diffusion ...

Ab initio molecular dynamics simulation is used to investigate the structure and dynamics of liquid Se at temperatures of 870 and 1370~K. The calculated static structure factor is in excellent agreement with experimental data. The calculated radial distribution function gives a mean coordination number close to 2, but we find a significant fraction of one-fold and three-fold atoms, particularly at 1370~K, so that the chain structure is considerably disrupted. The self-diffusion coefficient has values ( ∼1× 10 −8 ~m~s −1 ) typical of liquid metals.

Helium atom scattering (HAS) studies of the H-covered Mo(110) and W(110) surfaces reveal a twofold anomaly in the respective dispersion curves. In order to explain this unusual behavior we performed density functional theory calculations of the atomic and electronic structure, the vibrational properties, and the spectrum of electron-hole excitations of those surfaces. Our work provides evidence for hydrogen adsorption induced Fermi surface nesting. The respective nesting vectors are in excellent agreement with the HAS data and recent angle resolved photoemission experiments of the H-covered alloy system Mo_0.95Re_0.05(110). Also, we investigated the electron-phonon coupling and discovered that the Rayleigh phonon frequency is lowered for those critical wave vectors. Moreover, the smaller indentation in the HAS spectra can be clearly identified as a Kohn anomaly. Based on our results for the susceptibility and the recently improved understanding of the He scattering mechanism we argue that the larger anomalous dip is due to a direct interaction of the He atoms with electron-hole excitations at the Fermi level.

Using density-functional-theory we study the electronic and structural properties of a monolayer of Cu on the fcc (100) and (111) surfaces of the late 4d transition metals, as well as a monolayer of Pd on Mo bcc(110). We calculate the ground states of these systems, as well as the difference of the ionization energies of an adlayer core electron and a core electron of the clean surface of the adlayer metal. The theoretical results are compared to available experimental data and discussed in a simple physical picture; it is shown why and how adlayer core-level binding energy shifts can be used to deduce information on the adlayer's chemical reactivity.

The growth morphology of clean silver exhibits a profound anisotropy: The growing surface of Ag(111) is typically very rough while that of Ag(100) is smooth and flat. This serious and important difference is unexpected, not understood, and hitherto not observed for any other metal. Using density functional theory calculations of self-diffusion on flat and stepped Ag(100) we find, for example, that at flat regions a hopping mechanism is favored, while across step edges diffusion proceeds by an exchange process. The calculated microscopic parameters explain the experimentally reported growth properties.