  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.
