A Pipek-Mezey style localization scheme has also been recently proposed for obtaining Wannier functions. Contrary to the maximally localized Wannier functions (which are an application of the Foster-Boys scheme to crystalline systems), the Pipek-Mezey Wannier functions do not mix σ and π orbitals.