Adding a new dimension to DFT calculations of solids ...
The Full Potential APW methods
Recently, the development of the Augmented Plane Wave (APW) methods from Slater's APW, to LAPW and the new APW+lo was described by Schwarz et al. 2001.
The LAPW method
The linearized augmented plane wave (LAPW) method is among the most accurate methods for performing electronic structure calculations for crystals. It is based on the density functional theory for the treatment of exchange and correlation and uses e.g. the local spin density approximation (LSDA). Several forms of LSDA potentials exist in the literature , but recent improvements using the generalized gradient approximation (GGA) are available too. For valence states relativistic effects can be included either in a scalar relativistic treatment (Koelling and Harmon 77) or with the second variational method including spin-orbit coupling (Macdonald 80, Novak 97). Core states are treated fully relativistically (Desclaux 69).
A description of this method to linearize Slater's old APW method (i.e. the LAPW formalism) and further programming hints are found in many references: Andersen 73, 75, Koelling 72, Koelling and Arbman 75, Wimmer et al. 81, Weinert 81, Weinert et al. 82, Blaha and Schwarz 83, Blaha et al. 85, Wei et al. 85, Mattheiss and Hamann 86, Jansen and Freeman 84, Schwarz and Blaha 96). An excellent book by D. Singh (Singh 94) describes all the details of the LAPW method and is highly recommended to the interested reader. Here only the basic ideas are summarized; details are left to those references.
Like most ``energy-band methods``, the LAPW method is a procedure for solving the Kohn-Sham equations for the ground state density, total energy, and (Kohn-Sham) eigenvalues (energy bands) of a many-electron system (here a crystal) by introducing a basis set which is especially adapted to the problem.
This adaptation is achieved by dividing the unit cell into (I) non-overlapping atomic spheres (centered at the atomic sites) and (II) an interstitial region. In the two types of regions different basis sets are used:
The solutions to the Kohn-Sham equations are expanded in this combined basis set of LAPW's according to the linear variation method
and the coefficients cn are determined by the Rayleigh-Ritz variational principle. The convergence of this basis set is controlled by a cutoff parameter RmtKmax= 6 - 9, where Rmt is the smallest atomic sphere radius in the unit cell and Kmax is the magnitude of the largest K vector.
In order to improve upon the linearization (i.e. to increase the flexibility of the basis) and to make possible a consistent treatment of semicore and valence states in one energy window (to ensure orthogonality) additional (kn independent) basis functions can be added. They are called "local orbitals" (Singh 91) and consist of a linear combination of 2 radial functions at 2 different energies (e.g. at the 3s and 4s energy) and one energy derivative (at one of these energies):
The coefficients Alm, Blm, and Clm, are determined by the requirements that should be normalized and has zero value and slope at the sphere boundary.
The APW+lo method
Sjöstedt, Nordström and Singh (2000) have shown that the standard LAPW method with the additional constraint on the PWs of matching in value AND slope to the solution inside the sphere is not the most efficient way to linearize Slater's APW method. It can be made much more efficient when one uses the standard APW basis, but of course with ul(r,El) at a fixed energy El in order to keep the linear eigenvalue problem. One then adds a new local orbital (lo) to have enough variational flexibility in the radial basisfunctions:
Thus we construct basis functions that have ``kinks'' at the sphere boundary, which makes it necessary to include surface terms in the kinetic energy part of the Hamiltonian. Note, however, that the total wavefunction is of course smooth and differentiable.
As shown by Madsen et al. (2001) this new scheme converges practically to identical results as the LAPW method, but allows to reduce ``RKmax'' by about one, leading to significantly smaller basis sets (up to 50 %) and thus the corresponding the computational time is drastically reduced (up to an order of magnitude). Within one calculation a mixed ``LAPW and APW+lo'' basis can be used for different atoms and even different l-values for the same atom (Madsen et al. 2001). In general one describes by APW+lo those orbitals which converge most slowly with the number of PWs (such as TM 3d states) or the atoms with a small sphere size, but the rest with ordinary LAPWs. One can also add a second ``lo'' at a different energy so that both, semicore and valence states, can be described simultaneously.
In its general form the LAPW method expands the potential in the following form
and the charge densities analogously. Thus no shape approximations are made, a procedure frequently called the ``full-potential`` method.
The ``muffin-tin`` approximation used in early band calculations corresponds to retaining only the L=0 and M=0 component in the first expression of equ. (9) and only the K=0 component in the second. This (much older) procedure corresponds to taking the spherical average inside the spheres and the volume average in the interstitial region.
The total energy is computed according to Weinert et al. 82.
Rydberg atomic units are used except internally in the atomic-like programs (LSTART and LCORE) or in subroutine outwin (LAPW1, LAPW2), where Hartree units are used. The output is always given in Rydberg units.
The forces at the atoms are calculated according to Yu et al (91). For the implementation of this formalism in WIEN see Kohler et al (94) and Madsen et al. 2001. An alternative formulation by Soler and Williams (89) has also been tested and found to be equivalent, both in computationally efficiency and numerical accuracy and the respective code is available from M.Fähnle (Krimmel et al 94).
The Fermi energy and the weights of each band state can be calculated using a modified tetrahedron method (Bloechl et al. 94), a Gaussian or a temperature broadening scheme.
Spin-orbit interactions can be considered via a second variational step using the scalar-relativistic eigenfunctions as basis. (See MacDonald 80, Singh 94 and Novak 97). In order to overcome the problems due to the missing p1/2 radial basis function in the scalar-relativistic basis (which corresponds to p3/2, we have recently extended the standard LAPW basis by an additional p1/2-local orbital'', i.e. a LO with a p1/2 basis function, which is added in the second-variational SO calculation (Kunes et al. 2001).
It is well known that for localized electrons (like the 4f states in lanthanides or 3d states in some TM-oxides) the LDA (GGA) method is not accurate enough for a proper description. Thus we have implemented various forms of the LDA+U method as well as the ``Orbital polarization method'' (OP) (see Novak 2001 and references therein).
One can also consider interactions with an external magnetic (see Novak 2001) or electric field (via a supercell approach, see Stahn et al. 2000).
The density of states (DOS) can be calculated using the modified tetrahedron method of Blöchl et al. 94.
X-ray absorption and emission spectra are determined using Fermi's golden rule and dipole matrix elements (between a core and valence or conduction band state respectively). (Neckel et al. 75)
X-ray structure factors are obtained by Fourier Transformation of the charge density.
Optical properties are obtained using the ``Joint density of states'' modified with the respective dipole matrix elements according to Ambrosch et al. 95, Abt et al. 94, Abt 97. A Kramers-Kronig transformation is also possible.
An analysis of the electron density according to Bader's ``atoms in molecules'' theory can be made using a program by J. Sofo and J. Fuhr (2001)
©2001 by P. Blaha and K. Schwarz