This method is based upon work on the projection of plane-wave calculations
onto atomic orbitals [166], which has been used to
analyse atomic basis-sets [167] and obtain local atomic
properties from the extended Kohn-Sham orbitals [168].
Here we review the method, for the special case of a -point
Brillouin zone sampling.