docs: def normalization_factor(self):
          """
          Sum of number of atoms minus the number of H and O in composition
          """
          return 1.0 / (self.num_atoms - self.composition.get('H', 0)
                        - self.composition.get('O', 0))