c_f = \frac{1}{\rho_f} \frac{\partial \rho_f}{\partial p} = -\frac{1}{V_{mf}} \frac{\partial V_{mf}}{\partial p}