A global investigation of phase equilibria using the perturbed-chain statistical-associating-fluid-theory approach

The recently developed perturbed-chain statistical-associating-fluid theory (PC-SAFT) is investigated for a wide range of model parameters including the parameter m representing the chain length and the thermodynamic temperature T and pressure p. This approach is based upon the first-order thermodynamic perturbation theory for chain molecules developed by Wertheim M. S. Wertheim, J. Stat. Phys. 35, 19 (1984); ibid. 42, 459 (1986)] and Chapman et al. G. Jackson, W. G. Chapman, and K. E. Gubbins, Mol. Phys. 65, 1 (1988); W. G. Chapman, G. Jackson, and K. E. Gubbins, ibid. 65, 1057 (1988)] and includes dispersion interactions via the second-order perturbation theory of Barker and Henderson J. A. Barker and D. Henderson, J. Chem. Phys. 47, 4714 (1967)]. We systematically study a hierarchy of models which are based on the PC-SAFT approach using analytical model calculations and Monte Carlo simulations. For one-component systems we find that the analytical model in contrast with the simulation results exhibits two phase-separation regions in addition to the common gas-liquid coexistence region: One phase separation occurs at high density and low temperature. The second demixing takes place at low density and high temperature where usually the ideal-gas phase is expected in the phase diagram. These phenomena, which are referred to as "liquid-liquid" and "gas-gas" equilibria, give rise to multiple critical points in one-component systems, as well as to critical end points and equilibria of three fluid phases, which can usually be found in multicomponent mixtures only. Furthermore, it is shown that the liquid-liquid demixing in this model is not a consequence of a "softened" repulsive interaction as assumed in the theoretical derivation of the model. Experimental data for the melt density of polybutadiene with molecular mass Mw=45,000 gmol are correlated here using the PC-SAFT equation. It is shown that the discrepancies in modeling the polymer density at ambient temperature and high pressure can be traced back to the liquid-liquid phase separation predicted by the equation of state at low temperatures. This investigation provides a basis for understanding possible inaccuracies or even unexpected phase behavior which can occur in engineering applications of the PC-SAFT model aiming at predicting properties of macromolecular substances.