Add thermodynamic potentials

Not ready to be merged. Opening for comments/suggestions.
are these copied over from SLB_2011.py?