if (itmp < 0 || itmp >= nreax_types)

      error->all("Non-existent ReaxFF type");

  }

  int n = atom->ntypes;