AFreeenergydensityf computation