Computing the Helmholtz energy function of solids

The procedure is based in the techniques of thermodynamic integration. The usual implementations derive from the paper by Frenkel and Ladd (See Ref.1)