Maximally localized Wannier functions within the (L)APW plus LO method

We present a robust algorithm that computes (maximally localized) Wannier functions without the need of providing an initial guess. Instead, a suitable starting point is constructed automatically from so-called local orbitals which are fundamental building blocks of the basis set within (linearized) augmented plane-wave methods. Our approach is applied to a vast variety of materials such as (semi)metals, bulk and low-dimensional semiconductors, and complex inorganic-organic hybrid interfaces. For the interpolation of electronic single-particle energies, an accuracy in the meV range can be easily achieved. We exemplify the capabilities of our method by the calculation of the joint density of states in aluminum, (generalized) Kohn-Sham and quasiparticle band structures in various semiconductors, and the electronic structure of beta-Ga2O3, including electron and hole effective masses.