skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Full potential multiple scattering theory

Abstract

A practical method for performing self-consistent electronic structure calculations based upon full-potential multiple-scattering theory is presented. Solutions to the single site Schroedinger equation are obtained by solving coupled channel integral equations for a potential which is analytically continued out to the circumscribing sphere. This potential coincides with the full cell potential inside each atomic cell. Scattering matrices and wavefunctions for the full cell potential are obtained from surface Wronskian relations. The charge density is obtained from the single particle Green`s function. This Green`s function is computed using the cell scattering matrices and wavefunctions using the layer multiple scattering theory. Self consistent solutions require a solution at each iteration to the Poisson equation. The Poisson equation is solved using a variational cellular method. In the approach a local solution to each cell is augmented by adding a series of regular harmonics (solutions to Laplace`s equation). Minimizing the coulomb energy, subject to continuity of the potential across all cell boundary provides an expression for the coefficients of the regular harmonics. This method is applied to BCC Nb. Calculated properties converge well in angular momentum and show comparable accuracy to full potential linearized muffin-tin orbital calculations.

Authors:
 [1]
  1. Tulane Univ., New Orleans, LA (United States). Dept. of Physics
Publication Date:
Research Org.:
Lawrence Livermore National Lab., CA (United States); Tulane Univ., New Orleans, LA (United States). Dept. of Physics
Sponsoring Org.:
USDOE, Washington, DC (United States)
OSTI Identifier:
10195709
Report Number(s):
UCRL-CR-116369
ON: DE95003197;; TRN: AHC29430%%76
DOE Contract Number:
W-7405-ENG-48
Resource Type:
Technical Report
Resource Relation:
Other Information: PBD: 20 Oct 1994
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE; NIOBIUM; LATTICE PARAMETERS; ELASTICITY; BCC LATTICES; ELECTRONIC STRUCTURE; CALCULATION METHODS; MUFFIN-TIN POTENTIAL; SELF-CONSISTENT FIELD; ALGORITHMS; THEORETICAL DATA; 360102; 360103; STRUCTURE AND PHASE STUDIES; MECHANICAL PROPERTIES

Citation Formats

MacLaren, J.M. Full potential multiple scattering theory. United States: N. p., 1994. Web. doi:10.2172/10195709.
MacLaren, J.M. Full potential multiple scattering theory. United States. doi:10.2172/10195709.
MacLaren, J.M. Thu . "Full potential multiple scattering theory". United States. doi:10.2172/10195709. https://www.osti.gov/servlets/purl/10195709.
@article{osti_10195709,
title = {Full potential multiple scattering theory},
author = {MacLaren, J.M.},
abstractNote = {A practical method for performing self-consistent electronic structure calculations based upon full-potential multiple-scattering theory is presented. Solutions to the single site Schroedinger equation are obtained by solving coupled channel integral equations for a potential which is analytically continued out to the circumscribing sphere. This potential coincides with the full cell potential inside each atomic cell. Scattering matrices and wavefunctions for the full cell potential are obtained from surface Wronskian relations. The charge density is obtained from the single particle Green`s function. This Green`s function is computed using the cell scattering matrices and wavefunctions using the layer multiple scattering theory. Self consistent solutions require a solution at each iteration to the Poisson equation. The Poisson equation is solved using a variational cellular method. In the approach a local solution to each cell is augmented by adding a series of regular harmonics (solutions to Laplace`s equation). Minimizing the coulomb energy, subject to continuity of the potential across all cell boundary provides an expression for the coefficients of the regular harmonics. This method is applied to BCC Nb. Calculated properties converge well in angular momentum and show comparable accuracy to full potential linearized muffin-tin orbital calculations.},
doi = {10.2172/10195709},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Oct 20 00:00:00 EDT 1994},
month = {Thu Oct 20 00:00:00 EDT 1994}
}

Technical Report:

Save / Share:
  • The single scattering optical potential was in the multiple scattering approaches of Watson and KMT. Since the kinematics of single scattering is three-body in nature, one builds a three-body model of this term. This approach can include the proper kinematics for the struck nucleon, the identity of the target nucleons, and the binding interaction of the struck nucleon. Integral equations of the Faddeev type are derived for both the Watson and KMT single scattering optical potentials. Unitarity relations are investigated and one observes that these relations can be expanded in order to identify the intermediate states responsible for the absorptivemore » parts. The transition amplitudes to the inelastic states implicit in the model are extracted and evaluated. This permits one to understand the physical meaning of the imaginary part in precise terms. The same procedure is applied to the closure and impulse approximations for the single scattering term and their implicit inelastic states and reaction amplitudes are identified. These approximations are evaluated by analyzing the inelastic data. It is concluded that the impulse approximation to the Watson single scattering term should provide the best two-body approximation to a single-scattering optical potential.« less
  • One-electron Green s functions play a central role in multiple-scattering theory (MST) based electronic- structure methods. Robust methods exist for calculating the Green s function for crystal potentials that are spherically symmetric about atomic centers. When applied to potentials of general shape, these same techniques result in pathologies in the small-r behavior of the electronic charge density because a portion of the Green s function can become singular at the origin for that case. We propose an algebraic method that eliminates the singular behavior by making use of the equivalence of two terms that involve poles in the inverse ofmore » the sine matrix. Our accurate calculations illustrate the limitations of previous methods for treating this problem that rely on extrapolating the solutions near the origin.« less
  • The Green function plays an essential role in the Korringa–Kohn–Rostoker(KKR) multiple scattering method. In practice, it is constructed from the regular and irregular solutions of the local Kohn–Sham equation and robust methods exist for spherical potentials. However, when applied to a non-spherical potential, numerical errors from the irregular solutions give rise to pathological behaviors of the charge density at small radius. Here we present a full-potential implementation of the fully-relativistic KKR method to perform ab initio self-consistent calculation by directly solving the Dirac differential equations using the generalized variable phase (sine and cosine matrices) formalism Liu et al. (2016). Themore » pathology around the origin is completely eliminated by carrying out the energy integration of the single-site Green function along the real axis. Here, by using an efficient pole-searching technique to identify the zeros of the well-behaved Jost matrices, we demonstrated that this scheme is numerically stable and computationally efficient, with speed comparable to the conventional contour energy integration method, while free of the pathology problem of the charge density. As an application, this method is utilized to investigate the crystal structures of polonium and their bulk properties, which is challenging for a conventional real-energy scheme. The noble metals are also calculated, both as a test of our method and to study the relativistic effects.« less