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)


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
Issue number1
Publication statusPublished - 1993


ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this