code: def natoms(self):
          
          syms = [site.specie.symbol for site in self.structures[0]]
          return [len(tuple(a[1])) for a in itertools.groupby(syms)]
docs: Sequence of number of sites of each type associated with the Poscar.
          Similar to 7th line in vasp 5+ Xdatcar.