An efficient step-forward way to solve the Schrödinger eigenvalue equation in self-consistent calculations

Ruqian Wu, Arthur J Freeman

Research output: Contribution to journalArticle

9 Citations (Scopus)

Abstract

A simple practical step-forward procedure is employed for the rapid diagonalization of large matrices occurring in ab initio pseudopotential (PP) and linearized augmented plane wave (LAPW) calculations in the process leading to the charge density self-consistency. Much smaller energy cut-offs can be adopted to reduce the dimensionality of the eigenvalue equation by transforming the Hamiltonian and overlap (for LAPW) matrices into a diagonally-dominant form with the eigenvectors obtained in previous iterations. For the FLAPW method, the possible variation of basis functions in the muffin-tin and vacuum regions is treated by adopting both the fixed basis approach and the direct step-forward approach. Test calculations carried out for an Fe monolayer and the MgO (001) surface indicate that this procedure works very well for both metallic and ionic systems.

Original languageEnglish
Pages (from-to)58-62
Number of pages5
JournalComputer Physics Communications
Volume76
Issue number1
DOIs
Publication statusPublished - 1993

Fingerprint

plane waves
eigenvalues
Hamiltonians
matrices
Charge density
Eigenvalues and eigenfunctions
Tin
pseudopotentials
iteration
Monolayers
tin
eigenvectors
Vacuum
vacuum
energy

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

@article{0c90f8b7a6b5481884bbdd0d6c81a706,
title = "An efficient step-forward way to solve the Schr{\"o}dinger eigenvalue equation in self-consistent calculations",
abstract = "A simple practical step-forward procedure is employed for the rapid diagonalization of large matrices occurring in ab initio pseudopotential (PP) and linearized augmented plane wave (LAPW) calculations in the process leading to the charge density self-consistency. Much smaller energy cut-offs can be adopted to reduce the dimensionality of the eigenvalue equation by transforming the Hamiltonian and overlap (for LAPW) matrices into a diagonally-dominant form with the eigenvectors obtained in previous iterations. For the FLAPW method, the possible variation of basis functions in the muffin-tin and vacuum regions is treated by adopting both the fixed basis approach and the direct step-forward approach. Test calculations carried out for an Fe monolayer and the MgO (001) surface indicate that this procedure works very well for both metallic and ionic systems.",
author = "Ruqian Wu and Freeman, {Arthur J}",
year = "1993",
doi = "10.1016/0010-4655(93)90120-2",
language = "English",
volume = "76",
pages = "58--62",
journal = "Computer Physics Communications",
issn = "0010-4655",
publisher = "Elsevier",
number = "1",

}

TY - JOUR

T1 - An efficient step-forward way to solve the Schrödinger eigenvalue equation in self-consistent calculations

AU - Wu, Ruqian

AU - Freeman, Arthur J

PY - 1993

Y1 - 1993

N2 - A simple practical step-forward procedure is employed for the rapid diagonalization of large matrices occurring in ab initio pseudopotential (PP) and linearized augmented plane wave (LAPW) calculations in the process leading to the charge density self-consistency. Much smaller energy cut-offs can be adopted to reduce the dimensionality of the eigenvalue equation by transforming the Hamiltonian and overlap (for LAPW) matrices into a diagonally-dominant form with the eigenvectors obtained in previous iterations. For the FLAPW method, the possible variation of basis functions in the muffin-tin and vacuum regions is treated by adopting both the fixed basis approach and the direct step-forward approach. Test calculations carried out for an Fe monolayer and the MgO (001) surface indicate that this procedure works very well for both metallic and ionic systems.

AB - A simple practical step-forward procedure is employed for the rapid diagonalization of large matrices occurring in ab initio pseudopotential (PP) and linearized augmented plane wave (LAPW) calculations in the process leading to the charge density self-consistency. Much smaller energy cut-offs can be adopted to reduce the dimensionality of the eigenvalue equation by transforming the Hamiltonian and overlap (for LAPW) matrices into a diagonally-dominant form with the eigenvectors obtained in previous iterations. For the FLAPW method, the possible variation of basis functions in the muffin-tin and vacuum regions is treated by adopting both the fixed basis approach and the direct step-forward approach. Test calculations carried out for an Fe monolayer and the MgO (001) surface indicate that this procedure works very well for both metallic and ionic systems.

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

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

U2 - 10.1016/0010-4655(93)90120-2

DO - 10.1016/0010-4655(93)90120-2

M3 - Article

VL - 76

SP - 58

EP - 62

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

IS - 1

ER -