For the vapor–liquid equilibrium a heterogeneous approach was adopted: the physical interactions between all species are taken into account through activity coefficients calculation in the liquid phase and through an equation of state for the gas phase. In this study, chemical reaction equilibrium constants are not considered as variables.