Add thermodynamic potentials

Not ready to be merged. Opening for comments/suggestions.
I am worried about this function.  In the paper, it says that this is the thermal expansivity at ambient pressure.  Can you comment on whether this works for high pressure?