code: def get_helmholtz_energy(self, temperature, electronic_energy=0, verbose=False):
          
          thermo_object = HarmonicThermo(vib_energies=self.vib_energies,
                                         potentialenergy=electronic_energy)
          self.helmholtz_energy = thermo_object.get_helmholtz_energy(temperature=temperature, verbose=verbose)
          return (self.helmholtz_energy)
docs: Returns the Helmholtz energy of an adsorbed molecule.
  
          Parameters
          ----------
          temperature : numeric
              temperature in K
          electronic_energy : numeric
              energy in eV
          verbose : boolean
              whether to print ASE thermochemistry output
  
          Returns
          -------
          helmholtz_energy : numeric
              Helmholtz energy in eV