Based on the best of my knowledge, fugacity coefficient for a component in a mixture is generally obtained from EOS, whereas the activity coefficient of a component in a solution is calculated using gibbs excess function. A united system of nonideality description for mixtures would be interesting and efficient from the thermodynamic point of view. 

Similar questions and discussions