Linearized augmented-plane-wave method

The linearized augmented-plane-wave method (LAPW) is an implementation of Kohn-Sham density functional theory (DFT) adapted to periodic materials.[1][2][3] It typically goes along with the treatment of both valence and core electrons on the same footing in the context of DFT and the treatment of the full potential and charge density without any shape approximation. This is often referred to as the all-electron full-potential linearized augmented-plane-wave method (FLAPW).[4] It does not rely on the pseudopotential approximation and employs a systematically extendable basis set. These features make it one of the most precise implementations of DFT, applicable to all crystalline materials, regardless of their chemical composition. It can be used as a reference for evaluating other approaches.[5][6]

Introduction

At the core of density functional theory the Hohenberg-Kohn theorems state that every observable of an interacting many-electron system is a functional of its ground-state charge density and that this density minimizes the total energy of the system.[7] The theorems do not answer the question how to obtain such a ground-state density. A recipe for this is given by Walter Kohn and Lu Jeu Sham who introduce an auxiliary system of noninteracting particles constructed such that it shares the same ground-state density with the interacting particle system.[8] The Schrödinger-like equations describing this system are the Kohn-Sham equations. With these equations one can calculate the eigenstates of the system and with these the density. One contribution to the Kohn-Sham equations is the effective potential which itself depends on the density. As the ground-state density is not known before a Kohn-Sham DFT calculation and it is an input as well as an output of such a calculation, the Kohn-Sham equations are solved in an iterative procedure together with a recalculation of the density and the potential in every iteration. It starts with an initial guess for the density and after every iteration a new density is constructed as a mixture from the output density and previous densities. The calculation finishes as soon as a fixpoint of a self-consistent density is found, i.e., input and output density are identical. This is the ground-state density.

A method implementing Kohn-Sham DFT has to realize these different steps of the sketched iterative algorithm. The LAPW method is based on a partitioning of the material's unit cell into non-overlapping but nearly touching so-called muffin-tin (MT) spheres, centered at the atomic nuclei, and an interstitial region (IR) in between the spheres. The physical description and the representation of the Kohn-Sham orbitals, the charge density, and the potential is adapted to this partitioning. In the following this method design and the extraction of quantities from it are sketched in more detail. Variations and extensions are indicated.

Solving the Kohn-Sham equations

The central aspect of practical DFT implementations is the question how to solve the Kohn-Sham equations

${\displaystyle \left[{\hat {T}}_{\text{s}}+V_{\text{eff}}(\mathbf {r} )\right]\left|\Psi _{j}^{\mathbf {k} }(\mathbf {r} )\right\rangle =\epsilon _{j}^{\mathbf {k} }\left|\Psi _{j}^{\mathbf {k} }(\mathbf {r} )\right\rangle }$

with the single-electron kinetic energy operator ${\displaystyle {\hat {T}}_{\text{s}}}$, the effective potential ${\displaystyle V_{\text{eff}}(\mathbf {r} )}$, Kohn-Sham states ${\displaystyle \Psi _{j}^{\mathbf {k} }(\mathbf {r} )}$, energy eigenvalues ${\displaystyle \epsilon _{j}^{\mathbf {k} }}$, and position and Bloch vectors ${\displaystyle \mathbf {r} }$ and ${\displaystyle \mathbf {k} }$. While in abstract evaluations of Kohn-Sham DFT the model for the exchange-correlation contribution to the effective potential is the only fundamental approximation, in practice solving the Kohn-Sham equations is accompanied by the introduction of many additional approximations. These include the incompleteness of the basis set used to represent the Kohn-Sham orbitals, the choice of whether to use the pseudopotential approximation or to consider all electrons in the DFT scheme, the treatment of relativistic effects, and possible shape approximations to the potential. Beyond the partitioning of the unit cell, for the LAPW method the central design aspect is the use of the LAPW basis set ${\displaystyle \left\lbrace \phi _{\mathbf {k} ,\mathbf {G} }(\mathbf {r} )\right\rbrace }$ to represent the valence electron orbitals as

${\displaystyle \left|\Psi _{j}^{\mathbf {k} }(\mathbf {r} )\right\rangle =\sum \limits _{\mathbf {G} }c_{j}^{\mathbf {k} ,\mathbf {G} }\left|\phi _{\mathbf {k} ,\mathbf {G} }(\mathbf {r} )\right\rangle ,}$

where ${\displaystyle c_{j}^{\mathbf {k} ,\mathbf {G} }}$ are the expansion coefficients. The LAPW basis is designed to enable a precise representation of the orbitals and an accurate modelling of the physics in each region of the unit cell.

Considering a unit cell of volume ${\displaystyle \Omega }$ covering atoms ${\displaystyle \alpha }$ at positions ${\displaystyle \mathbf {\tau } _{\alpha }}$, an LAPW basis function is characterized by a reciprocal lattice vector ${\displaystyle \mathbf {G} }$ and the considered Bloch vector ${\displaystyle \mathbf {k} }$. It is given as

${\displaystyle \phi _{\mathbf {k} ,\mathbf {G} }(\mathbf {r} )=\left\lbrace {\begin{array}{l l}{\frac {1}{\sqrt {\Omega }}}e^{i(\mathbf {k} +\mathbf {G} )\mathbf {r} }&{\text{for }}\mathbf {r} {\text{ in IR}}\\\sum \limits _{l=0}^{l_{{\text{max}},\alpha }}\sum \limits _{m=-l}^{l}\left[a_{l,m}^{\mathbf {k} ,\mathbf {G} ,\alpha }u_{l,\alpha }(r_{\alpha },E_{l,\alpha })+b_{l,m}^{\mathbf {k} ,\mathbf {G} ,\alpha }{\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right]Y_{l,m}(\mathbf {\hat {r}} _{\alpha })&{\text{for }}\mathbf {r} {\text{ in MT}}_{\alpha }\end{array}}\right.,}$

where ${\displaystyle \mathbf {r} _{\alpha }=\mathbf {r} -\mathbf {\tau } _{\alpha }}$ is the position vector relative to the position of atom nucleus ${\displaystyle \alpha }$. An LAPW basis function is thus a plane wave in the IR and a linear combination of the radial functions ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ and ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ multiplied by spherical harmonics ${\displaystyle Y_{l,m}}$ in each MT sphere. The radial function ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ is hereby the solution of the Kohn-Sham Hamiltonian for the spherically averaged potential with regular behavior at the nucleus for the given energy parameter ${\displaystyle E_{l,\alpha }}$. Together with its energy derivative ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ these augmentations of the plane wave in each MT sphere enable a representation of the Kohn-Sham orbitals at arbitrary eigenenergies linearized around the energy parameters. The coefficients ${\displaystyle a_{l,m}^{\mathbf {k} ,\mathbf {G} ,\alpha }}$ and ${\displaystyle b_{l,m}^{\mathbf {k} ,\mathbf {G} ,\alpha }}$ are automatically determined by enforcing the basis function to be continuously differentiable for the respective ${\displaystyle (l,m)}$ channel. The set of LAPW basis functions is defined by specifying a cutoff parameter ${\displaystyle K_{\text{max}}=|\mathbf {k} +\mathbf {G} |_{\text{max}}}$. In each MT sphere, the expansion into spherical harmonics is limited to a maximum number of angular momenta ${\displaystyle l_{{\text{max}},\alpha }\approx K_{\text{max}}R_{{\text{MT}}_{\alpha }}}$, where ${\displaystyle R_{{\text{MT}}_{\alpha }}}$ is the muffin-tin radius of atom ${\displaystyle \alpha }$. The choice of this cutoff is connected to the decay of expansion coefficients for growing ${\displaystyle l}$ in the Rayleigh expansion of plane waves into spherical harmonics.

While the LAPW basis functions are used to represent the valence states, core electron states, which are completely confined within a MT sphere, are calculated for the spherically averaged potential on radial grids, for each atom separately applying atomic boundary conditions. Semicore states, which are still localized but slightly extended beyond the MT sphere boundary, may either be treated as core electron states or as valence electron states. For the latter choice the linearized representation is not sufficient because the related eigenenergy is typically far away from the energy parameters. To resolve this problem the LAPW basis can be extended by additional basis functions in the respective MT sphere, so called local orbitals (LOs).[9] These are tailored to provide a precise representation of the semicore states.

The plane-wave form of the basis functions in the interstitial region makes setting up the Hamiltonian matrix

${\displaystyle H_{\mathbf {G'} ,\mathbf {G} }^{\mathbf {k} }=\left\langle \phi _{\mathbf {k} ,\mathbf {G'} }{\Big |}{\hat {H}}{\Big |}\phi _{\mathbf {k} ,\mathbf {G} }\right\rangle =\left\langle \phi _{\mathbf {k} ,\mathbf {G'} }{\Big |}{\hat {T}}_{\text{s}}+V_{\text{eff}}(\mathbf {r} ){\Big |}\phi _{\mathbf {k} ,\mathbf {G} }\right\rangle }$

for that region simple. In the MT spheres this setup is also simple and computationally inexpensive for the kinetic energy and the spherically averaged potential, e.g., in the muffin-tin approximation. The simplicity hereby stems from the connection of the radial functions to the spherical Hamiltonian in the spheres ${\displaystyle {\hat {H}}_{\text{sphr}}^{\alpha }}$, i.e., ${\displaystyle {\hat {H}}_{\text{sphr}}^{\alpha }\left|u_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right\rangle =E_{l,\alpha }\left|u_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right\rangle }$ and ${\displaystyle {\hat {H}}_{\text{sphr}}^{\alpha }\left|{\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right\rangle =E_{l,\alpha }\left|{\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right\rangle +\left|u_{l,\alpha }(r_{\alpha },E_{l,\alpha })\right\rangle }$. In comparison to the MT approximation, for the full-potential description (FLAPW) contributions from the non-spherical part of the potential are added to the Hamiltonian matrix in the MT spheres and in the IR contributions related to deviations from the constant potential.

After the Hamiltonian matrix ${\displaystyle H_{\mathbf {G'} ,\mathbf {G} }^{\mathbf {k} }}$ together with the overlap matrix ${\displaystyle S_{\mathbf {G'} ,\mathbf {G} }^{\mathbf {k} }=\left\langle \phi _{\mathbf {k} ,\mathbf {G'} }{\Big |}\phi _{\mathbf {k} ,\mathbf {G} }\right\rangle }$ is set up, the Kohn-Sham orbitals are obtained as eigenfunctions from the algebraic generalized dense Hermitian eigenvalue problem

${\displaystyle \sum \limits _{\mathbf {G} }H_{\mathbf {G'} ,\mathbf {G} }^{\mathbf {k} }c_{j}^{\mathbf {k} ,\mathbf {G} }=\epsilon _{j}^{\mathbf {k} }\sum \limits _{\mathbf {G} }S_{\mathbf {G'} ,\mathbf {G} }^{\mathbf {k} }c_{j}^{\mathbf {k} ,\mathbf {G} }~~,}$

where ${\displaystyle \epsilon _{j}^{\mathbf {k} }}$ is the energy eigenvalue of the j-th Kohn-Sham state at Bloch vector ${\displaystyle {\mathbf {k} }}$ and the state is given as indicated above by the expansion coefficients ${\displaystyle c_{j}^{\mathbf {k} ,\mathbf {G} }}$.

The considered degree of relativistic physics differs for core and valence electrons. The strong localization of core electrons due to the singularity of the effective potential at the atomic nucleus is connected to large kinetic energy contributions and thus a fully relativistic treatment is desirable and common. For the determination of the radial functions ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ and ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ the common approach is to make an approximation to the fully relativistic description. This may be the scalar-relativistic approximation[10][11] (SRA) or similar approaches.[12][13] The dominant effect neglected by these approximations is the spin-orbit coupling. As indicated above the construction of the Hamiltonian matrix within such an approximation is trivial. Spin-orbit coupling can additionally be included, though this leads to a more complex Hamiltonian matrix setup or a second variation scheme,[14][15] connected to increased computational demands. In the interstitial region it is reasonable and common to describe the valence electrons without considering relativistic effects.

Representation of the charge density and the potential

After calculating the Kohn-Sham eigenfunctions, the next step is to construct the electron charge density by occupying the lowest energy eigenstates up to the Fermi level with electrons. The Fermi level itself is determined in this process by keeping charge neutrality in the unit cell. The resulting charge density ${\displaystyle \rho (\mathbf {r} )}$ then has a region-specific form

${\displaystyle \rho (\mathbf {r} )=\left\{{\begin{array}{l l}\sum \limits _{\mathbf {G} }\rho _{\mathbf {G} }e^{i\mathbf {G} \mathbf {r} }&{\text{for }}\mathbf {r} {\text{ in IR}}\\\sum \limits _{l=0}^{l_{{\text{max}},\alpha }}\sum \limits _{m=-l}^{l}\rho _{l,m}^{\alpha }(r_{\alpha })Y_{l,m}(\mathbf {\hat {r}} _{\alpha })&{\text{for }}\mathbf {r} {\text{ in MT}}_{\alpha }\end{array}}\right.,}$

i.e., it is given as a plane-wave expansion in the interstitial region and as an expansion into radial functions times spherical harmonics in each MT sphere. The radial functions hereby are numerically given on a mesh.

The representation of the effective potential follows the same scheme. In its construction a common approach is to employ Weinert's method for solving the Poisson equation.[16] It efficiently and accurately provides a solution of the Poisson equation without shape approximation for an arbitrary periodic charge density based on the concept of multipole potentials and the boundary value problem for a sphere.

Postprocessing and extracting results

Because they are based on the same theoretical framework, different DFT implementations offer access to very similar sets of material properties. However, the variations in the implementations result in differences in the ease of extracting certain quantities and also in differences in their interpretation. In the following, these circumstances are sketched for some examples.

The most basic quantity provided by DFT is the ground-state total energy of an investigated system. To avoid the calculation of derivatives of the eigenfunctions in its evaluation, the common implementation[17] replaces the expectation value of the kinetic energy operator by the sum of the band energies of occupied Kohn-Sham states minus the energy due to the effective potential. The force exerted on an atom, which is given by the change of the total energy due to an infinitesimal displacement, has two major contributions. The first contribution is due to the displacement of the potential. It is known as Hellmann-Feynman force. The other, computationally more elaborate contribution, is due to the related change in the atom-position-dependent basis functions. It is often called Pulay force[18] and requires a method-specific implementation.[19][20] Beyond forces, similar method-specific implementations are also needed for further quantities derived from the total energy functional. For the LAPW method, formulations for the stress tensor[21] and for phonons[22] have been realized.

Independent of the actual size of an atom, evaluating atom-dependent quantities in LAPW is often interpreted as calculating the quantity in the respective MT sphere. This applies to quantities like charges at atoms, magnetic moments, or projections of the density of states or the band structure onto a certain orbital character at a given atom. Deviating interpretations of such quantities from experiments or other DFT implementations may lead to differences when comparing results. On a side note also some atom-specific LAPW inputs relate directly to the respective MT region. For example, in the DFT+U approach the Hubbard U only affects the MT sphere.[23]

A strength of the LAPW approach is the inclusion of all electrons in the DFT calculation, which is crucial for the evaluation of certain quantities. One of which are hyperfine interaction parameters like electric field gradients whose calculation involves the evaluation of the curvature of the all-electron Coulomb potential near the nuclei. The prediction of such quantities with LAPW is very accurate.[24]

Kohn-Sham DFT does not give direct access to all quantities one may be interested in. For example, most energy eigenvalues of the Kohn-Sham states are not directly related to the real interacting many-electron system. For the prediction of optical properties one therefore often uses DFT codes in combination with software implementing the GW approximation (GWA) to many-body perturbation theory and optionally the Bethe-Salpeter equation (BSE) to describe excitons. Such software has to be adapted to the representation used in the DFT implementation. Both the GWA and the BSE have been formulated in the LAPW context and several implementations of such tools are in use.[25][26][27] In other postprocessing situations it may be useful to project Kohn-Sham states onto Wannier functions. For the LAPW method such projections have also been implemented[28][29] and are in common use.

Variants and extensions of the LAPW method

• APW: The augmented-plane-wave method[30] is the predecessor of LAPW. It uses the radial solution to the spherically averaged potential for the augmentation in the MT spheres. The energy derivative of this radial function is not involved. This missing linearization implies that the augmentation has to be adapted to each Kohn-Sham state individually, i.e., it depends on the Bloch vector and the band index, which subsequently leads to a non-linear, energy-dependent eigenvalue problem. In comparison to LAPW this is a more complex problem to solve. A relativistic generalization of this approach, RAPW, has also been formulated.[31]
• Local orbitals extensions: The LAPW basis can be extended by local orbitals[9] (LOs). These are additional basis functions having nonvanishing values only in a single MT sphere. They are composed of the radial functions ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$, ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$, and a third radial function tailored to describe use-case-specific physics. LOs have originally been proposed for the representation of semicore states. Other uses involve the representation of unoccupied states[32][33] or the elimination of the linearization error for the valence states.[34]
• APW+lo: In the APW+lo method[35][36] the augmentation in the MT spheres only consists of the function ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$. It is matched to the plane wave in the interstitial region only in value. As an alternative implementation of the linearization the function ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ is included in the basis set as an additional local orbital. While the matching conditions result in an unphysical kink of the basis functions at the MT sphere boundaries, a careful consideration of the kink in the construction of the Hamiltonian matrix suppresses it in the Kohn-Sham eigenfunctions. In comparison to the classical LAPW method the APW+lo approach leads to a less stiff basis set. The outcome is a faster convergence of the DFT calculations with respect to the basis set size.
• Soler-Williams formulation of LAPW: In the Soler-Williams formulation of LAPW[37] the plane waves cover the whole unit cell. In the MT spheres the augmentation is implemented by replacing up to the angular momentum cutoff the plane waves by the functions ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$ and ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$. This yields basis functions continuously differentiable also in the ${\displaystyle (l,m)}$ channels above the angular momentum cutoff. As a consequence the Soler-Williams approach has reduced angular momentum cutoff requirements in comparison to the classical LAPW formulation.
• ELAPW: In the extended LAPW method[38][39] pairs of local orbitals introducing the functions ${\displaystyle u_{l,\alpha }(r_{\alpha },E_{l,\alpha }^{\text{lo}})}$ and ${\displaystyle {\dot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha }^{\text{lo}})}$ are added to the LAPW basis. The energy parameters ${\displaystyle E_{l,\alpha }^{\text{lo}}}$ are chosen to systematically extend the energy region in which Kohn-Sham states are accurately described by the linearization in LAPW.
• QAPW: In the quadratic APW method[40][41] the augmentation in the MT spheres additionally includes the second energy derivative ${\displaystyle {\ddot {u}}_{l,\alpha }(r_{\alpha },E_{l,\alpha })}$. The matching at the MT sphere boundaries is performed by enforcing continuity of the basis functions in value, slope, and curvature. This is similar to the super-linearized APW (SLAPW) method in which radial functions ${\displaystyle u_{l,\alpha }}$ and/or their derivatives ${\displaystyle {\dot {u}}_{l,\alpha }}$ at more than one energy parameter are used for the augmentation.[9] In comparison to a pure LAPW basis these approaches can precisely represent Kohn-Sham orbitals in a broader energy window around the energy parameters. The drawback is that the stricter matching conditions lead to a stiffer basis set.
• Lower-dimensional systems: The partitioning of the unit cell can be extended to explicitly include semi-infinite vacuum regions with their own augmentations of the plane waves.[42] This enables efficient calculations for lower-dimensional systems such as surfaces and thin films. For the treatment of atomic chains an extension to one-dimensional setups has been formulated.[43]

Software implementations

There are various software projects implementing the LAPW method and/or its variants. Examples for such codes are

• Elk
• Exciting
• Flair[44]
• FLEUR
• HiLAPW
• Wien2k

References

1. Marcus, Paul M. (January 1967). "Variational methods in the computation of energy bands". International Journal of Quantum Chemistry. 1 (S1): 567–588. doi:10.1002/qua.560010659.
2. Andersen, O. Krogh (15 October 1975). "Linear methods in band theory" (PDF). Physical Review B. 12 (8): 3060–3083. Bibcode:1975PhRvB..12.3060A. doi:10.1103/PhysRevB.12.3060. S2CID 122641701.
3. Koelling, D. D.; Arbman, G. O. (November 1975). "Use of energy derivative of the radial solution in an augmented plane wave method: application to copper". Journal of Physics F: Metal Physics. 5 (11): 2041–2054. Bibcode:1975JPhF....5.2041K. doi:10.1088/0305-4608/5/11/016. ISSN 0305-4608.
4. Wimmer, E.; Krakauer, H.; Weinert, M.; Freeman, A. J. (15 July 1981). "Full-potential self-consistent linearized-augmented-plane-wave method for calculating the electronic structure of molecules and surfaces: O 2 molecule". Physical Review B. 24 (2): 864–875. Bibcode:1981PhRvB..24..864W. doi:10.1103/PhysRevB.24.864.
5. Lejaeghere, K.; Bihlmayer, G.; Bjorkman, T.; Blaha, P.; Blugel, S.; Blum, V.; Caliste, D.; Castelli, I. E.; Clark, S. J.; Dal Corso, A.; de Gironcoli, S.; Deutsch, T.; Dewhurst, J. K.; Di Marco, I.; Draxl, C.; Dułak, M.; Eriksson, O.; Flores-Livas, J. A.; Garrity, K. F.; Genovese, L.; Giannozzi, P.; Giantomassi, M.; Goedecker, S.; Gonze, X.; Granas, O.; Gross, E. K. U.; Gulans, A.; Gygi, F.; Hamann, D. R.; Hasnip, P. J.; Holzwarth, N. A. W.; Iuşan, D.; Jochym, D. B.; Jollet, F.; Jones, D.; Kresse, G.; Koepernik, K.; Kucukbenli, E.; Kvashnin, Y. O.; Locht, I. L. M.; Lubeck, S.; Marsman, M.; Marzari, N.; Nitzsche, U.; Nordstrom, L.; Ozaki, T.; Paulatto, L.; Pickard, C. J.; Poelmans, W.; Probert, M. I. J.; Refson, K.; Richter, M.; Rignanese, G.-M.; Saha, S.; Scheffler, M.; Schlipf, M.; Schwarz, K.; Sharma, S.; Tavazza, F.; Thunstrom, P.; Tkatchenko, A.; Torrent, M.; Vanderbilt, D.; van Setten, M. J.; Van Speybroeck, V.; Wills, J. M.; Yates, J. R.; Zhang, G.-X.; Cottenier, S. (25 March 2016). "Reproducibility in density functional theory calculations of solids". Science. 351 (6280): aad3000. Bibcode:2016Sci...351.....L. doi:10.1126/science.aad3000. hdl:1854/LU-7191263. PMID 27013736. S2CID 206642768.
6. Bosoni, Emanuele; Beal, Louis; Bercx, Marnik; Blaha, Peter; Blügel, Stefan; Bröder, Jens; Callsen, Martin; Cottenier, Stefaan; Degomme, Augustin; Dikan, Vladimir; Eimre, Kristjan; Flage-Larsen, Espen; Fornari, Marco; Garcia, Alberto; Genovese, Luigi; Giantomassi, Matteo; Huber, Sebastiaan P.; Janssen, Henning; Kastlunger, Georg; Krack, Matthias; Kresse, Georg; Kühne, Thomas D.; Lejaeghere, Kurt; Madsen, Georg K. H.; Marsman, Martijn; Marzari, Nicola; Michalicek, Gregor; Mirhosseini, Hossein; Müller, Tiziano M. A.; Petretto, Guido; Pickard, Chris J.; Poncé, Samuel; Rignanese, Gian-Marco; Rubel, Oleg; Ruh, Thomas; Sluydts, Michael; Vanpoucke, Danny E. P.; Vijay, Sudarshan; Wolloch, Michael; Wortmann, Daniel; Yakutovich, Aliaksandr V.; Yu, Jusong; Zadoks, Austin; Zhu, Bonan; Pizzi, Giovanni (January 2024). "How to verify the precision of density-functional-theory implementations via reproducible and universal workflows". Nature Reviews Physics. 6 (1): 45–58. arXiv:2305.17274. doi:10.1038/s42254-023-00655-3.
7. Hohenberg, P.; Kohn, W. (9 November 1964). "Inhomogeneous Electron Gas". Physical Review. 136 (3B): B864–B871. Bibcode:1964PhRv..136..864H. doi:10.1103/PhysRev.136.B864. S2CID 117806275.
8. Kohn, W.; Sham, L. J. (15 November 1965). "Self-Consistent Equations Including Exchange and Correlation Effects". Physical Review. 140 (4A): A1133–A1138. Bibcode:1965PhRv..140.1133K. doi:10.1103/PhysRev.140.A1133.
9. Singh, David (15 March 1991). "Ground-state properties of lanthanum: Treatment of extended-core states". Physical Review B. 43 (8): 6388–6392. Bibcode:1991PhRvB..43.6388S. doi:10.1103/PhysRevB.43.6388. PMID 9998076.
10. Koelling, D D; Harmon, B N (28 August 1977). "A technique for relativistic spin-polarised calculations". Journal of Physics C: Solid State Physics. 10 (16): 3107–3114. Bibcode:1977JPhC...10.3107K. doi:10.1088/0022-3719/10/16/019.
11. Takeda, T. (March 1978). "The scalar relativistic approximation". Zeitschrift für Physik B. 32 (1): 43–48. Bibcode:1978ZPhyB..32...43T. doi:10.1007/BF01322185. S2CID 120097976.
12. van Lenthe, E.; Snijders, J. G.; Baerends, E. J. (15 October 1996). "The zero‐order regular approximation for relativistic effects: The effect of spin–orbit coupling in closed shell molecules". The Journal of Chemical Physics. 105 (15): 6505–6516. Bibcode:1996JChPh.105.6505V. doi:10.1063/1.472460. hdl:11370/2909724f-afb3-482e-a4e5-f8440d7dab09.
13. Dyall, Kenneth G.; van Lenthe, Erik (22 July 1999). "Relativistic regular approximations revisited: An infinite-order relativistic approximation". The Journal of Chemical Physics. 111 (4): 1366–1372. Bibcode:1999JChPh.111.1366D. doi:10.1063/1.479395.
14. MacDonald, A H; Picket, W E; Koelling, D D (20 May 1980). "A linearised relativistic augmented-plane-wave method utilising approximate pure spin basis functions". Journal of Physics C: Solid State Physics. 13 (14): 2675–2683. Bibcode:1980JPhC...13.2675M. doi:10.1088/0022-3719/13/14/009.
15. Li, Chun; Freeman, A. J.; Jansen, H. J. F.; Fu, C. L. (15 September 1990). "Magnetic anisotropy in low-dimensional ferromagnetic systems: Fe monolayers on Ag(001), Au(001), and Pd(001) substrates". Physical Review B. 42 (9): 5433–5442. doi:10.1103/PhysRevB.42.5433. PMID 9996126.
16. Weinert, M. (November 1981). "Solution of Poisson's equation: Beyond Ewald‐type methods". Journal of Mathematical Physics. 22 (11): 2433–2439. Bibcode:1981JMP....22.2433W. doi:10.1063/1.524800.
17. Weinert, M.; Wimmer, E.; Freeman, A. J. (15 October 1982). "Total-energy all-electron density functional method for bulk solids and surfaces". Physical Review B. 26 (8): 4571–4578. Bibcode:1982PhRvB..26.4571W. doi:10.1103/PhysRevB.26.4571.
18. Pulay, P. (January 1969). "Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules: I. Theory". Molecular Physics. 17 (2): 197–204. Bibcode:1969MolPh..17..197P. doi:10.1080/00268976900100941.
19. Yu, Rici; Singh, D.; Krakauer, H. (15 March 1991). "All-electron and pseudopotential force calculations using the linearized-augmented-plane-wave method". Physical Review B. 43 (8): 6411–6422. Bibcode:1991PhRvB..43.6411Y. doi:10.1103/PhysRevB.43.6411. PMID 9998079.
20. Klüppelberg, Daniel A.; Betzinger, Markus; Blügel, Stefan (5 January 2015). "Atomic force calculations within the all-electron FLAPW method: Treatment of core states and discontinuities at the muffin-tin sphere boundary". Physical Review B. 91 (3): 035105. Bibcode:2015PhRvB..91c5105K. doi:10.1103/PhysRevB.91.035105.
21. Belbase, Kamal; Tröster, Andreas; Blaha, Peter (29 November 2021). "Stress tensor in the linearized augmented plane wave method". Physical Review B. 104 (17): 174113. doi:10.1103/PhysRevB.104.174113.
22. Gerhorst, Christian-Roman; Neukirchen, Alexander; Klüppelberg, Daniel A; Bihlmayer, Gustav; Betzinger, Markus; Michalicek, Gregor; Wortmann, Daniel; Blügel, Stefan (1 March 2024). "Phonons from density-functional perturbation theory using the all-electron full-potential linearized augmented plane-wave method FLEUR *". Electronic Structure. 6 (1): 017001. arXiv:2309.14799. doi:10.1088/2516-1075/ad1614.
23. Shick, A. B.; Liechtenstein, A. I.; Pickett, W. E. (15 October 1999). "Implementation of the LDA+U method using the full-potential linearized augmented plane-wave basis". Physical Review B. 60 (15): 10763–10769. arXiv:cond-mat/9903439. Bibcode:1999PhRvB..6010763S. doi:10.1103/PhysRevB.60.10763. S2CID 119508105.
24. Blaha, P.; Dufek, P.; Schwarz, K.; Haas, H. (December 1996). "Calculation of electric hyperfine interaction parameters in solids". Hyperfine Interactions. 97–98 (1): 1–10. Bibcode:1996HyInt..97....1B. doi:10.1007/BF02150164. S2CID 95082151.
25. Friedrich, Christoph; Blügel, Stefan; Schindlmayr, Arno (3 March 2010). "Efficient implementation of the G W approximation within the all-electron FLAPW method". Physical Review B. 81 (12): 125102. arXiv:1003.0316. Bibcode:2010PhRvB..81l5102F. doi:10.1103/PhysRevB.81.125102. S2CID 43385321.
26. Jiang, Hong; Gómez-Abal, Ricardo I.; Li, Xin-Zheng; Meisenbichler, Christian; Ambrosch-Draxl, Claudia; Scheffler, Matthias (February 2013). "FHI-gap: A GW code based on the all-electron augmented plane wave method". Computer Physics Communications. 184 (2): 348–366. Bibcode:2013CoPhC.184..348J. doi:10.1016/j.cpc.2012.09.018.
27. Gulans, Andris; Kontur, Stefan; Meisenbichler, Christian; Nabok, Dmitrii; Pavone, Pasquale; Rigamonti, Santiago; Sagmeister, Stephan; Werner, Ute; Draxl, Claudia (10 September 2014). "exciting: a full-potential all-electron package implementing density-functional theory and many-body perturbation theory". Journal of Physics: Condensed Matter. 26 (36): 363202. Bibcode:2014JPCM...26J3202G. doi:10.1088/0953-8984/26/36/363202. PMID 25135665. S2CID 29889985.
28. Freimuth, F.; Mokrousov, Y.; Wortmann, D.; Heinze, S.; Blügel, S. (17 July 2008). "Maximally localized Wannier functions within the FLAPW formalism". Physical Review B. 78 (3): 035120. arXiv:0806.3213. Bibcode:2008PhRvB..78c5120F. doi:10.1103/PhysRevB.78.035120. S2CID 53133273.
29. Kuneš, Jan; Arita, Ryotaro; Wissgott, Philipp; Toschi, Alessandro; Ikeda, Hiroaki; Held, Karsten (November 2010). "Wien2wannier: From linearized augmented plane waves to maximally localized Wannier functions". Computer Physics Communications. 181 (11): 1888–1895. arXiv:1004.3934. Bibcode:2010CoPhC.181.1888K. doi:10.1016/j.cpc.2010.08.005. S2CID 27628281.
30. Slater, J. C. (15 May 1937). "Wave Functions in a Periodic Potential". Physical Review. 51 (10): 846–851. Bibcode:1937PhRv...51..846S. doi:10.1103/PhysRev.51.846.
31. Loucks, T. L. (16 August 1965). "Relativistic Electronic Structure in Crystals. I. Theory". Physical Review. 139 (4A): A1333–A1337. Bibcode:1965PhRv..139.1333L. doi:10.1103/PhysRev.139.A1333.
32. Betzinger, Markus; Friedrich, Christoph; Blügel, Stefan; Görling, Andreas (25 January 2011). "Local exact exchange potentials within the all-electron FLAPW method and a comparison with pseudopotential results". Physical Review B. 83 (4): 045105. arXiv:1102.2379. Bibcode:2011PhRvB..83d5105B. doi:10.1103/PhysRevB.83.045105. S2CID 118450370.
33. Friedrich, Christoph; Schindlmayr, Arno; Blügel, Stefan; Kotani, Takao (7 July 2006). "Elimination of the linearization error in GW calculations based on the linearized augmented-plane-wave method". Physical Review B. 74 (4): 045104. arXiv:cond-mat/0606605. Bibcode:2006PhRvB..74d5104F. doi:10.1103/PhysRevB.74.045104. S2CID 53596360.
34. Michalicek, Gregor; Betzinger, Markus; Friedrich, Christoph; Blügel, Stefan (December 2013). "Elimination of the linearization error and improved basis-set convergence within the FLAPW method". Computer Physics Communications. 184 (12): 2670–2679. arXiv:1302.3130. Bibcode:2013CoPhC.184.2670M. doi:10.1016/j.cpc.2013.07.002. S2CID 7540554.
35. Sjöstedt, E; Nordström, L; Singh, D.J (March 2000). "An alternative way of linearizing the augmented plane-wave method". Solid State Communications. 114 (1): 15–20. Bibcode:2000SSCom.114...15S. doi:10.1016/S0038-1098(99)00577-3.
36. Madsen, Georg K. H.; Blaha, Peter; Schwarz, Karlheinz; Sjöstedt, Elisabeth; Nordström, Lars (30 October 2001). "Efficient linearization of the augmented plane-wave method". Physical Review B. 64 (19): 195134. Bibcode:2001PhRvB..64s5134M. doi:10.1103/PhysRevB.64.195134.
37. Soler, José M.; Williams, Arthur R. (15 July 1989). "Simple formula for the atomic forces in the augmented-plane-wave method". Physical Review B. 40 (3): 1560–1564. Bibcode:1989PhRvB..40.1560S. doi:10.1103/PhysRevB.40.1560. PMID 9992009.
38. Krasovskii, E.E.; Yaresko, A.N.; Antonov, V.N. (May 1994). "Theoretical study of ultraviolet photoemission spectra of noble metals". Journal of Electron Spectroscopy and Related Phenomena. 68: 157–166. doi:10.1016/0368-2048(94)02113-9.
39. Krasovskii, E. E. (15 November 1997). "Accuracy and convergence properties of the extended linear augmented-plane-wave method". Physical Review B. 56 (20): 12866–12873. Bibcode:1997PhRvB..5612866K. doi:10.1103/PhysRevB.56.12866.
40. Smrčka, L. (July 1984). "Linearized augmented plane wave method utilizing the quadratic energy expansion of radial wave functions". Czechoslovak Journal of Physics. 34 (7): 694–704. Bibcode:1984CzJPh..34..694S. doi:10.1007/BF01589865. S2CID 121089065.
41. Petrů, J.; Smrčka, L. (January 1985). "Quadratic augmented plane wave method for self-consistent band structure calculations". Czechoslovak Journal of Physics. 35 (1): 62–71. Bibcode:1985CzJPh..35...62P. doi:10.1007/BF01590276. S2CID 120146231.
42. Krakauer, H.; Posternak, M.; Freeman, A. J. (15 February 1979). "Linearized augmented plane-wave method for the electronic band structure of thin films". Physical Review B. 19 (4): 1706–1719. Bibcode:1979PhRvB..19.1706K. doi:10.1103/PhysRevB.19.1706.
43. Mokrousov, Y.; Bihlmayer, G.; Blügel, S. (July 2005). "Full-potential linearized augmented plane-wave method for one-dimensional systems: Gold nanowire and iron monowires in a gold tube" (PDF). Physical Review B. 72 (4): 045402. Bibcode:2005PhRvB..72d5402M. doi:10.1103/PhysRevB.72.045402.
44. Weinert, M; Schneider, G; Podloucky, R; Redinger, J (25 February 2009). "FLAPW: applications and implementations". Journal of Physics: Condensed Matter. 21 (8): 084201. Bibcode:2009JPCM...21h4201W. doi:10.1088/0953-8984/21/8/084201. PMID 21817353. S2CID 13318752.