docs: def get_vib_energies(self):
          """Returns a list of vibration in energy units eV.
  
          Returns
          -------
          vibs : list of vibrations in eV
          """
          vibs = self.molecule_dict[self.name]['vibrations']
          vibs = np.array(vibs) * cm2ev
          return (vibs)