Evaluating perturbation contributions in SAFT models by comparing to molecular simulation of n-alkanes

• Published in: Fluid phase equilibria 2011, Vol.306 (1), p.57-66
• Main Authors: Ghobadi, Ahmadreza F., Elliott, J. Richard
• Format: Article
• Language: English
• Subjects: Chemical thermodynamics; Chemistry; DMD simulation; Exact sciences and technology; General and physical chemistry; General. Theory; n-Alkanes; SAFT equation of state; Second order Thermodynamic Perturbation Theory
• ISSN: 0378-3812; 1879-0224
• DOI: 10.1016/j.fluid.2011.03.016

► Constraints are shown for trends of terms in thermodynamic perturbation theory (TPT). ► The critical region is strongly impacted by accurate characterization of TPT trends. ► Molecular simulation is extended to polymers by analyzing trends in TPT terms. ► Details of the PC-SAFT model can be refined by analyzing trends from molecular simulation. ► Significant improvement is provided by the SAFT-D formulation of Chapman et al. SAFT models are generally written as a perturbation series of the Helmholtz energy with reciprocal temperature as the argument. The perturbation coefficients are then functions of density and molecular size. The variation of the perturbation coefficients with molecular size is given primarily by Wertheim's theory [6–9], but there may be additional variations as in the PC-SAFT model. In the present work, we compare the characterization of perturbation coefficients inferred from PC-SAFT to those derived from molecular simulations. The molecular simulations are based on Discontinuous Molecular Dynamics (DMD) and second order Thermodynamic Perturbation Theory (TPT). DMD simulation is applied to the repulsive part of the potential model with molecular details like fused hard spheres for the interaction sites and 110° bond angles. The thermodynamic effects of disperse attractions are treated by rigorous application of TPT. The present work re-examines the related work of Elliott and Gray [35] in the low density and critical regions, focusing on n-alkanes with carbon numbers ranging from 3 to 80. We find that SAFT theory overestimates the repulsive contribution ( A 0) and underestimates the first order contribution ( A 1) of Helmholtz energy relative to simulation. Nevertheless, the correlations are qualitatively reasonable. Significant inconsistencies arise when considering the second order contribution ( A 2). For example, the PC-SAFT characterization of A 2 becomes larger than A 1 in the low density, long chain limit, raising concerns about the convergence of the series. Furthermore, fluctuations are underestimated in the critical region and overestimated in the liquid region. In each case, we can suggest improved characterizations. Altogether, these results suggest ways to modify the SAFT formalism to achieve greater consistency between atomistic and coarse-grained models.