Construction of Maximally Localized Wannier Functions
  We present a general method of constructing maximally localized Wannier functions. It consists of three steps: (1) picking a localized trial wave function, (2) performing a full band projection, and (3) orthonormalizing with the Lowdin method. Our method is capable of producing maximally localized Wannier functions without further minimization, and it can be applied straightforwardly to random potentials without using supercells. The effectiveness of our method is demonstrated for both simple bands and composite bands. 
Junbo Zhu, Zhu Chen, Biao Wu
cond-mat.other