A more specific description of the orbital localization can be achieved through decomposition of IAO charges in terms of their shells, where a shell is defined as a full set of basis functions centered on the same atom for a given angular momentum (see Refs. 23,42). The penalty function S i then becomes