K-point integration: Difference between revisions

From VASP Wiki
No edit summary
Line 43: Line 43:
== Linear tetrahedron methods ==
== Linear tetrahedron methods ==


Within the linear tetrahedron method, the term $\eps_\nk$ is interpolated
Within the linear tetrahedron method, the term <math>\epsilon_{n\bold{k}}</math> is interpolated
linearly between two k-points. Bloechl~\cite{blo94} has recently revised
linearly between two k-points. Bloechl{{cite|bloechl:prb:1994}} has recently revised
the tetrahedron method to give effective weights $f(\{\eps_\nk\})$
the tetrahedron method to give effective weights <math>f(\{\epsilon_{n\bold{k}}\})</math>
for each band and k-point. In addition Bloechel was
for each band and k-point. In addition Bloechel was
able to derive a correction formula which removes the
able to derive a correction formula which removes the

Revision as of 10:55, 19 March 2019

In this section we discuss partial occupancies. A must for all readers.

First there is the question why to use partial occupancies at all. The answer is: partial occupancies help to decrease the number of k-points necessary to calculate an accurate band-structure energy. This answer might be strange at first sight. What we want to calculate is, the integral over the filled parts of the bands

where is the Dirac step function. Due to our finite computer resources this integral has to be evaluated using a discrete set of k-points[1][2][3]:

Keeping the step function we get a sum

which converges exceedingly slow with the number of k-points included. This slow convergence speed arises only from the fact that the occupancies jump form 1 to 0 at the Fermi-level. If a band is completely filled the integral can be calculated accurately using a low number of k-points (this is the case for semiconductors and insulators).

For metals the trick is now to replace the step function by a (smooth) function resulting in a much faster convergence speed without destroying the accuracy of the sum. Several methods have been proposed to solve this dazzling problem.

Linear tetrahedron methods

Within the linear tetrahedron method, the term is interpolated linearly between two k-points. Bloechl[4] has recently revised the tetrahedron method to give effective weights for each band and k-point. In addition Bloechel was able to derive a correction formula which removes the quadratic error inherent in the linear tetrahedron method (linear tetrahedron method with Bloechel corrections). The linear tetrahedron is more or less fool proof and requires a minimal interference by the user.

The main drawback is that the Bloechels method is not variational with respect to the partial occupancies if the correction terms are included, therefore the calculated forces might be wrong by a few percent. If accurate forces are required we recommend a finite temperature method.

References