docs: def mw(mol, ndigits=2):
      """Return standard molecular weight
      :param ndigits: number of digits
      """
      mol.require("Valence")
      return round(sum(a.mw() for _, a in mol.atoms_iter()), ndigits)