Electronic structure of thin films by the self-consistent numerical-basis-set linear combination of atomic orbitals method

Ni(001)

C. S. Wang*, A. J. Freeman

*Corresponding author for this work

Research output: Contribution to journalArticle

95 Citations (Scopus)

Abstract

We present the self-consistent numerical-basis-set linear combination of atomic orbitals (LCAO) discrete variational method for treating the electronic structure of thin films. As in the case of bulk solids, this method provides for thin films accurate solutions of the one-particle local density equations with a non-muffintin potential. Hamiltonian and overlap matrix elements are evaluated accurately by means of a three-dimensional numerical Diophantine integration scheme. Application of this method is made to the self-consistent solution of one-, three-, and five-layer Ni(001) unsupported films. The LCAO Bloch basis set consists of valence orbitals (3d, 4s, and 4p states for transition metals) orthogonalized to the frozen-core wave functions. The self-consistent potential is obtained iteratively within the superposition of overlapping spherical atomic charge density model with the atomic configurations treated as adjustable parameters. Thus the crystal Coulomb potential is constructed as a superposition of overlapping spherically symmetric atomic potentials and, correspondingly, the local density Kohn-Sham (=23) potential is determined from a superposition of atomic charge densities. At each iteration in the self-consistency procedure, the crystal charge density is evaluated using a sampling of 15 independent k points in (18)th of the irreducible two-dimensional Brillouin zone. The total density of states (DOS) and projected local DOS (by layer plane) are calculated using an analytic linear energy triangle method (presented as an Appendix) generalized from the tetrahedron scheme for bulk systems. Distinct differences are obtained between the surface and central plane local DOS. The central plane DOS is found to converge rapidly to the DOS of bulk paramagnetic Ni obtained by Wang and Callaway. Unlike the result of earlier non-self-consistent calculations, only a very small surplus charge (0.03 electron/atom) is found on the surface planes"in agreement with jellium model calculations.

Original languageEnglish (US)
Pages (from-to)793-805
Number of pages13
JournalPhysical Review B
Volume19
Issue number2
DOIs
StatePublished - Jan 1 1979

Fingerprint

electronic structure
orbitals
thin films
Coulomb potential
Brillouin zones
triangles
tetrahedrons
crystals
iteration
transition metals
sampling
wave functions
valence
matrices
configurations
atoms
electrons

ASJC Scopus subject areas

  • Condensed Matter Physics

Cite this

@article{5ec45f2977fe48f282caee70ff08925c,
title = "Electronic structure of thin films by the self-consistent numerical-basis-set linear combination of atomic orbitals method: Ni(001)",
abstract = "We present the self-consistent numerical-basis-set linear combination of atomic orbitals (LCAO) discrete variational method for treating the electronic structure of thin films. As in the case of bulk solids, this method provides for thin films accurate solutions of the one-particle local density equations with a non-muffintin potential. Hamiltonian and overlap matrix elements are evaluated accurately by means of a three-dimensional numerical Diophantine integration scheme. Application of this method is made to the self-consistent solution of one-, three-, and five-layer Ni(001) unsupported films. The LCAO Bloch basis set consists of valence orbitals (3d, 4s, and 4p states for transition metals) orthogonalized to the frozen-core wave functions. The self-consistent potential is obtained iteratively within the superposition of overlapping spherical atomic charge density model with the atomic configurations treated as adjustable parameters. Thus the crystal Coulomb potential is constructed as a superposition of overlapping spherically symmetric atomic potentials and, correspondingly, the local density Kohn-Sham (=23) potential is determined from a superposition of atomic charge densities. At each iteration in the self-consistency procedure, the crystal charge density is evaluated using a sampling of 15 independent k points in (18)th of the irreducible two-dimensional Brillouin zone. The total density of states (DOS) and projected local DOS (by layer plane) are calculated using an analytic linear energy triangle method (presented as an Appendix) generalized from the tetrahedron scheme for bulk systems. Distinct differences are obtained between the surface and central plane local DOS. The central plane DOS is found to converge rapidly to the DOS of bulk paramagnetic Ni obtained by Wang and Callaway. Unlike the result of earlier non-self-consistent calculations, only a very small surplus charge (0.03 electron/atom) is found on the surface planes{"}in agreement with jellium model calculations.",
author = "Wang, {C. S.} and Freeman, {A. J.}",
year = "1979",
month = "1",
day = "1",
doi = "10.1103/PhysRevB.19.793",
language = "English (US)",
volume = "19",
pages = "793--805",
journal = "Physical Review B",
issn = "0163-1829",
number = "2",

}

Electronic structure of thin films by the self-consistent numerical-basis-set linear combination of atomic orbitals method : Ni(001). / Wang, C. S.; Freeman, A. J.

In: Physical Review B, Vol. 19, No. 2, 01.01.1979, p. 793-805.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Electronic structure of thin films by the self-consistent numerical-basis-set linear combination of atomic orbitals method

T2 - Ni(001)

AU - Wang, C. S.

AU - Freeman, A. J.

PY - 1979/1/1

Y1 - 1979/1/1

N2 - We present the self-consistent numerical-basis-set linear combination of atomic orbitals (LCAO) discrete variational method for treating the electronic structure of thin films. As in the case of bulk solids, this method provides for thin films accurate solutions of the one-particle local density equations with a non-muffintin potential. Hamiltonian and overlap matrix elements are evaluated accurately by means of a three-dimensional numerical Diophantine integration scheme. Application of this method is made to the self-consistent solution of one-, three-, and five-layer Ni(001) unsupported films. The LCAO Bloch basis set consists of valence orbitals (3d, 4s, and 4p states for transition metals) orthogonalized to the frozen-core wave functions. The self-consistent potential is obtained iteratively within the superposition of overlapping spherical atomic charge density model with the atomic configurations treated as adjustable parameters. Thus the crystal Coulomb potential is constructed as a superposition of overlapping spherically symmetric atomic potentials and, correspondingly, the local density Kohn-Sham (=23) potential is determined from a superposition of atomic charge densities. At each iteration in the self-consistency procedure, the crystal charge density is evaluated using a sampling of 15 independent k points in (18)th of the irreducible two-dimensional Brillouin zone. The total density of states (DOS) and projected local DOS (by layer plane) are calculated using an analytic linear energy triangle method (presented as an Appendix) generalized from the tetrahedron scheme for bulk systems. Distinct differences are obtained between the surface and central plane local DOS. The central plane DOS is found to converge rapidly to the DOS of bulk paramagnetic Ni obtained by Wang and Callaway. Unlike the result of earlier non-self-consistent calculations, only a very small surplus charge (0.03 electron/atom) is found on the surface planes"in agreement with jellium model calculations.

AB - We present the self-consistent numerical-basis-set linear combination of atomic orbitals (LCAO) discrete variational method for treating the electronic structure of thin films. As in the case of bulk solids, this method provides for thin films accurate solutions of the one-particle local density equations with a non-muffintin potential. Hamiltonian and overlap matrix elements are evaluated accurately by means of a three-dimensional numerical Diophantine integration scheme. Application of this method is made to the self-consistent solution of one-, three-, and five-layer Ni(001) unsupported films. The LCAO Bloch basis set consists of valence orbitals (3d, 4s, and 4p states for transition metals) orthogonalized to the frozen-core wave functions. The self-consistent potential is obtained iteratively within the superposition of overlapping spherical atomic charge density model with the atomic configurations treated as adjustable parameters. Thus the crystal Coulomb potential is constructed as a superposition of overlapping spherically symmetric atomic potentials and, correspondingly, the local density Kohn-Sham (=23) potential is determined from a superposition of atomic charge densities. At each iteration in the self-consistency procedure, the crystal charge density is evaluated using a sampling of 15 independent k points in (18)th of the irreducible two-dimensional Brillouin zone. The total density of states (DOS) and projected local DOS (by layer plane) are calculated using an analytic linear energy triangle method (presented as an Appendix) generalized from the tetrahedron scheme for bulk systems. Distinct differences are obtained between the surface and central plane local DOS. The central plane DOS is found to converge rapidly to the DOS of bulk paramagnetic Ni obtained by Wang and Callaway. Unlike the result of earlier non-self-consistent calculations, only a very small surplus charge (0.03 electron/atom) is found on the surface planes"in agreement with jellium model calculations.

UR - http://www.scopus.com/inward/record.url?scp=0642374954&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0642374954&partnerID=8YFLogxK

U2 - 10.1103/PhysRevB.19.793

DO - 10.1103/PhysRevB.19.793

M3 - Article

VL - 19

SP - 793

EP - 805

JO - Physical Review B

JF - Physical Review B

SN - 0163-1829

IS - 2

ER -