No Arabic abstract
Highly accurate experimental structure factors of silicon are available in the literature, and these provide the ideal test for any emph{ab initio} method for the construction of the all-electron charge density. In a recent paper [J. R. Trail and D. M. Bird, Phys. Rev. B {bf 60}, 7863 (1999)] a method has been developed for obtaining an accurate all-electron charge density from a first principles pseudopotential calculation by reconstructing the core region of an atom of choice. Here this method is applied to bulk silicon, and structure factors are derived and compared with experimental and Full-potential Linear Augmented Plane Wave results (FLAPW). We also compare with the result of assuming the core region is spherically symmetric, and with the result of constructing a charge density from the pseudo-valence density + frozen core electrons. Neither of these approximations provide accurate charge densities. The aspherical reconstruction is found to be as accurate as FLAPW results, and reproduces the residual error between the FLAPW and experimental results.
The structure and properties of vacancies in a 2 nm Si nano-crystal are studied using a real space density functional theory/pseudopotential method. It is observed that a vacancys electronic properties and energy of formation are directly related to the local symmetry of the vacancy site. The formation energy for vacancies and Frenkel pair are calculated. It is found that both defects have lower energy in smaller crystals. In a 2 nm nano-crystal the energy to form a Frenkel pair is 1.7 eV and the energy to form a vacancy is no larger than 2.3 eV. The energy barrier for vacancy diffusion is examined via a nudged elastic band algorithm.
Quadratic-response theory is shown to provide a conceptually simple but accurate approximation for the self-consistent one-electron potential of semiconductor nanostructures. Numerical examples are presented for GaAs/AlAs and InGaAs/InP (001) superlattices using the local-density approximation to density-functional theory and norm-conserving pseudopotentials without spin-orbit coupling. When the reference crystal is chosen to be the virtual-crystal average of the two bulk constituents, the absolute error in the quadratic-response potential for Gamma(15) valence electrons is about 2 meV for GaAs/AlAs and 5 meV for InGaAs/InP. Low-order multipole expansions of the electron density and potential response are shown to be accurate throughout a small neighborhood of each reciprocal lattice vector, thus providing a further simplification that is confirmed to be valid for slowly varying envelope functions. Although the linear response is about an order of magnitude larger than the quadratic response, the quadratic terms are important both quantitatively (if an accuracy of better than a few tens of meV is desired) and qualitatively (due to their different symmetry and long-range dipole effects).
A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {bf 14}, 3795 (1981)]. In this method the core region is emph{reconstructed}, and none of the simplifying approximations (such as spherical symmetry of the charge density/potential or frozen core electrons) that previous solutions to this problem have required are made. The embedding method requires an accurate real space Green function, and an analysis of the errors introduced in constructing this from a set of numerical eigenstates is given. Results are presented for an all-electron reconstruction of bulk aluminium, for both the charge density and the density of states.
Linear scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high performance computers. The linear scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas is then discussed. The applications of linear scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear scaling methods are discussed.
We discuss an efficiency of various band structure algorithms in determining the Fermi surface (FS) of the paramagnetic ErGa3. The linear muffin-tin orbital (LMTO) in the atomic sphere approximation (ASA) method and three full potential (FP) codes: FP-LMTO, FP linear augmented plane wave (FLAPW), and FP local orbitals (FPLO) methods are employed. Results are compared with electron-positron (e-p) momentum densities reconstructed from two dimensional angular correlation of annihilation radiation (2D ACAR). Unexpectedly, none of used modern FP codes is able to give satisfying description of the experimental data that are in perfect agreement with LMTO-ASA results. We suspect that it can be connected with a different choice of the linearization energy.