returns Gibbs energy derivative (d2g/dd2)@T=const
function d2gd2T input EoS.HelmholtzDerivs f; output Der2EnergyByDensity2 d2gd2T; end d2gd2T;