Ist mit wenigen -Punkten ein Potential selbstkonsistent
bestimmt, wird mit diesem festen Potential weitergearbeitet: durch
Lösung der Matrixgleichung für eine größere Menge von -Punkten
werden Dispersionsrelationen für verschiedenen Symmetrierichtungen der
Brillouinzone bestimmt. Durch numerische Integration ergeben sich
totale und projizierte
Zustandsdichten.
Die totale Energie wird berechnet als Summe aus der kinetischen
Energie, die aus den Wellenfunktionen bestimmt wird, und verschiedenen
Beiträgen der potentiellen Energie, in die die
Ladungsdichteverteilung eingeht: innerhalb jeder atomaren Kugel sind
dies die Coulombenergie der Elektronen im Kernfeld, die
Hartree-Energie zwischen den Elektronen und die Austausch- und
Korrelationsenergie des Dichtefunktionals. Außerdem sind Beiträge
zwischen den abgeschirmten Kernen (Madelungterme) und zwischen den
elektronischen Ladungsdichten aus verschiedenen Kugeln zu
berücksichtigen. Die
Zuverlässigkeit der Energiewerte genügt jedoch meiner Erfahrung
nach nicht für die Zwecke von Strukturrelaxationen.