Modeling the phase behavior of H2S+n-alkane binary mixtures using the SAFT-VR+D approach

A statistical associating fluid theory for potential of variable range has been recently developed to model dipolar fluids (SAFT-VR+D) [Zhao and McCabe, J. Chem. Phys. 2006, 125, 104504]. The SAFT-VR+D equation explicitly accounts for dipolar interactions and their effect on the thermodynamics and structure of a fluid by using the generalized mean spherical approximation (GMSA) to describe a reference fluid of dipolar square-well segments. In this work, we apply the SAFT-VR+D approach to real mixtures of dipolar fluids. In particular, we examine the high-pressure phase diagram of hydrogen sulfide+n-alkane binary mixtures. Hydrogen sulfide is modeled as an associating spherical molecule with four off-center sites to mimic hydrogen bonding and an embedded dipole moment (micro) to describe the polarity of H2S. The n-alkane molecules are modeled as spherical segments tangentially bonded together to form chains of length m, as in the original SAFT-VR approach. By using simple Lorentz-Berthelot combining rules, the theoretical predictions from the SAFT-VR+D equation are found to be in excellent overall agreement with experimental data. In particular, the theory is able to accurately describe the different types of phase behavior observed for these mixtures as the molecular weight of the alkane is varied: type III phase behavior, according to the scheme of classification by Scott and Konynenburg, for the H2S+methane system, type IIA (with the presence of azeotropy) for the H2S+ethane and+propane mixtures; and type I phase behavior for mixtures of H2S and longer n-alkanes up to n-decane. The theory is also able to predict in a qualitative manner the solubility of hydrogen sulfide in heavy n-alkanes.     

Application of a renormalization-group treatment to the statistical associating fluid theory for potentials of variable range (SAFT-VR)

An accurate prediction of phase behavior at conditions far and close to criticality cannot be accomplished by mean-field based theories that do not incorporate long-range density fluctuations. A treatment based on renormalization-group (RG) theory as developed by White and co-workers has proven to be very successful in improving the predictions of the critical region with different equations of state. The basis of the method is an iterative procedure to account for contributions to the free energy of density fluctuations of increasing wavelengths. The RG method has been combined with a number of versions of the statistical associating fluid theory (SAFT), by implementing White's earliest ideas with the improvements of Prausnitz and co-workers. Typically, this treatment involves two adjustable parameters: a cutoff wavelength L for density fluctuations and an average gradient of the wavelet function Φ. In this work, the SAFT-VR (variable range) equation of state is extended with a similar crossover treatment which, however, follows closely the most recent improvements introduced by White. The interpretation of White's latter developments allows us to establish a straightforward method which enables Φ to be evaluated; only the cutoff wavelength L then needs to be adjusted. The approach used here begins with an initial free energy incorporating only contributions from short-wavelength fluctuations, which are treated locally. The contribution from long-wavelength fluctuations is incorporated through an iterative procedure based on attractive interactions which incorporate the structure of the fluid following the ideas of perturbation theories and using a mapping that allows integration of the radial distribution function. Good agreement close and far from the critical region is obtained using a unique fitted parameter L that can be easily related to the range of the potential. In this way the thermodynamic properties of a square-well (SW) fluid are given by the same number of independent intermolecular model parameters as in the classical equation. Far from the critical region the approach provides the correct limiting behavior reducing to the classical equation (SAFT-VR). In the critical region the β critical exponent is calculated and is found to take values close to the universal value. In SAFT-VR the free energy of an associating chain fluid is obtained following the thermodynamic perturbation theory of Wertheim from the knowledge of the free energy and radial distribution function of a reference monomer fluid. By determining L for SW fluids of varying well width a unique equation of state is obtained for chain and associating systems without further adjustment of critical parameters. We use computer simulation data of the phase behavior of chain and associating SW fluids to test the accuracy of the new equation.     

Predicting mixture phase equilibria and critical behavior using the SAFT-VRX approach

The SAFT-VRX equation of state combines the SAFT-VR equation with a crossover function that smoothly transforms the classical equation into a nonanalytical form close to the critical point. By a combinination of the accuracy of the SAFT-VR approach away from the critical region with the asymptotic scaling behavior seen at the critical point of real fluids, the SAFT-VRX equation can accurately describe the global fluid phase diagram. In previous work, we demonstrated that the SAFT-VRX equation very accurately describes the pvT and phase behavior of both nonassociating and associating pure fluids, with a minimum of fitting to experimental data. Here, we present a generalized SAFT-VRX equation of state for binary mixtures that is found to accurately predict the vapor-liquid equilibrium and pvT behavior of the systems studied. In particular, we examine binary mixtures of n-alkanes and carbon dioxide + n-alkanes. The SAFT-VRX equation accurately describes not only the gas-liquid critical locus for these systems but also the vapor-liquid equilibrium phase diagrams and thermal properties in single-phase regions.     

Equation of state and liquid-vapor equilibria of one- and two-Yukawa hard-sphere chain fluids: theory and simulation

The accuracy of several theories for the thermodynamic properties of the Yukawa hard-sphere chain fluid are studied. In particular, we consider the polymer mean spherical approximation (PMSA), the dimer version of thermodynamic perturbation theory (TPTD), and the statistical associating fluid theory for potentials of variable attractive range (SAFT-VR). Since the original version of SAFT-VR for Yukawa fluids is restricted to the case of one-Yukawa tail, we have extended SAFT-VR to treat chain fluids with two-Yukawa tails. The predictions of these theories are compared with Monte Carlo (MC) simulation data for the pressure and phase behavior of the chain fluid of different length with one- and two-Yukawa tails. We find that overall the PMSA and TPTD give more accurate predictions than SAFT-VR, and that the PMSA is slightly more accurate than TPTD.     

Solid-liquid equilibrium using the SAFT-VR equation of state: Solubility of naphthalene and acetic acid in binary mixtures and calculation of phase diagrams

A solid-liquid equilibrium (SLE) thermodynamic model based on the SAFT-VR equation of state (EOS) is presented. The model allows for the calculation of solid-liquid phase equilibria in binary mixtures at atmospheric pressure. The fluid (liquid) phase is treated with the SAFT-VR approach, where molecules are modelled as associating chains of tangentially bonded spherical segments interacting via square-well potentials of variable range. The equilibrium between the liquid and solid phase is treated following a standard thermodynamic method that requires the experimental temperature and enthalpy of fusion of the solute. The model is used to calculate the solubilities of naphthalene and acetic acid in common associating and non-associating organic solvents and to determine the solid-liquid phase behaviour of binary mixtures with simple eutectics. The SAFT-VR pure component model parameters are determined by comparison to experimental vapour pressure and saturated liquid density data with the choice of association models according to the nature of the molecule; in addition, an unlike adjustable parameter (k_ij) is used to model the solutions. The solubility data of naphthalene and acetic acid in both associating and non-associating solvents are reproduced essentially within the accuracy of the experimental measurements. The phase boundaries and the position of the eutectic points in the binary mixtures considered are, in most cases, reproduced with the accuracy commensurate with the industrial applications. Overall, the results presented show that the SAFT-VR EOS can be used with confidence for the prediction of the SLE of binary systems at atmospheric pressure.     

Semiclassical approach to model quantum fluids using the statistical associating fluid theory for systems with potentials of variable range

Thermodynamic properties of quantum fluids are described using an extended version of the statistical associating fluid theory for potentials of variable range (SAFT-VR) that takes into account quantum corrections to the Helmholtz free energy A, based on the Wentzel-Kramers-Brillouin approximation. We present the theoretical background of this approach (SAFT-VRQ), considering two different cases depending on the continuous or discontinuous nature of the particles pair interaction. For the case of continuous potentials, we demonstrate that the standard Wigner-Kirkwood theory for quantum fluids can be derived from the de Broglie-Bohm formalism for quantum mechanics that can be incorporated within the Barker and Henderson perturbation theory for liquids in a straightforward way. When the particles interact via a discontinuous pair potential, the SAFT-VR method can be combined with the perturbation theory developed by Singh and Sinha [J. Chem. Phys. 67, 3645 (1977); and ibid. 68, 562 (1978)]. We present an analytical expression for the first-order quantum perturbation term for a square-well potential, and the theory is applied to model thermodynamic properties of hydrogen, deuterium, neon, and helium-4. Vapor-liquid equilibrium, liquid and vapor densities, isochoric and isobaric heat capacities, Joule-Thomson coefficients and inversion curves are predicted accurately with respect to experimental data. We find that quantum corrections are important for the global behavior of properties of these fluids and not only for the low-temperature regime. Predictions obtained for hydrogen compare very favorably with respect to cubic equations of state.     

Thermodynamic modeling of hydrogen sulfide solubility in ionic liquids using modified SAFT-VR and PC-SAFT equations of state

Equations of state based on the statistical associating fluid theory for potentials of variable range (SAFT-VR) and the perturbed chain statistical associating fluid theory (PC-SAFT) have been used to model the PVT behavior of ionic liquids and the solubility of H₂S in six imidazolium-based ionic liquids. The studied systems included [bmim][PF6], [hmim][PF6], [bmim][BF4], [hmim][BF4], [bmim][NTF2] and [hmim][NTF2] at various temperatures and pressures.For pure components, parameters of the models have been obtained by fitting the models to experimental data on liquid densities; the average relative deviation between the calculated and experimental densities for ionic liquids is less than 2.42% in the PC-SAFT model and 5.44% in the SAFT-VR approach, the latter which incorporates the square-well potential for short-range interactions. In both models an additional term has been added to account for dipole-dipole interactions between solute molecules resulting from the permanent charges on the chain molecules of the solvents. The model parameters have also been correlated as functions of the molecular weight of the solvents. For binary mixtures of ionic liquids and H₂S, the association interactions between H₂S molecules and between the ionic liquids and H₂S molecules have also been taken into account in both approaches, using binary interaction coefficients. The results show an average deviation of less than 5% in the calculation of the mole fraction of H₂S in the ionic liquids. The effect of inclusion of the polar term has been studied for binary systems in both models.     

Phase behavior of self-associating fluids with weaker dispersion interactions between bonded particles

In this study, we explore the global phase behavior of a simple model for self-associating fluids where association reduces the strength of the dispersion interactions between bonded particles. Recent research shows that this type of behavior likely explains the thermodynamic properties of strongly polar fluids and certain micellar solutions. Based on Wertheim's theory of associating liquids [M. S. Wertheim, J. Stat. Phys. 42, 459 (1986); 42, 477 (1986)], our model takes into account the effect that dissimilar particle interactions have on the equilibrium constant for self-association in the system. We find that weaker interactions between bonded molecules tend to favor the dissociation of chains at any temperature and density. This effect stabilizes a monomeric liquid phase at high densities, enriching the global phase behavior of the system. In particular, for systems in which the energy of mixing between bonded and unbonded species is positive, we find a triple point involving a vapor, a dense phase of chain aggregates, and a monomeric liquid. Phase coexistence between the vapor and the monomeric fluid is always more stable at temperatures above the triple point, but a highly associated fluid may exist as a metastable phase under these conditions. The presence of this metastable phase may explain the characteristic nucleation behavior of the liquid phase in strongly dipolar fluids.     

Liquid mixtures of xenon with fluorinated species: xenon + sulfur hexafluoride

The vapor-liquid equilibrium of binary mixtures of xenon + SF6 has been measured at nine temperatures from 235.34 to 295.79 K and pressures up to 6.5 MPa. The mixture critical line is found to be continuous between the critical points of the pure components, and hence, the system can be classified as type I phase behavior in the scheme of van Konynenburg and Scott. The excess Gibbs free energies have been calculated, and the experimental results have been interpreted using the statistical associating fluid theory for potentials of variable range (SAFT-VR). Additionally, the SAFT-VR equation has been used to model other systems involving SF6 and alkanes, illustrating the predictability of the approach and further demonstrating the transferability of parameters between binary mixtures involving alkanes and xenon.     

Augmented van der Waals equations of state: SAFT-VR versus Yukawa based van der Waals equation

Performance of the SAFT-VR equation of state developed for the hard sphere based simple fluids, namely the square-well, Sutherland and Yukawa fluids, is examined by comparing its results with simulation data and an augmented van der Waals (vdW) equation based on a Yukawa (Y) reference. Its shown that both for the equilibrium vapor-liquid data and data along selected isotherms in the liquid and supercritical fluid phases the vdW(Y) equation provides better results, particularly when going to lower temperatures.     

Phase behavior of model confined fluids. Influence of substrate-fluid interaction strength

We examine the relationship between the macroscopic phase behavior of nanoconfined fluids and the nature of microscopic interactions between a confining substrate and fluid. Two model slit-pore systems are explored using grand canonical transition-matrix Monte Carlo simulation. One system consists of a square-well fluid interacting with a square-well substrate, and the other contains an embedded point charge model of lysozyme interacting with a mica surface. Fluid phase diagrams are constructed for a broad range of substrate conditions. Our results indicate that one observes a maximum in the critical temperature of the fluid phase envelope upon variation of substrate strength for a given slit width. Both systems studied exhibit such maxima at intermediate wall strength. The physical rationale for this observation suggests that this behavior should be generally expected. We introduce two metrics that enable one to predict conditions that produce maxima in critical temperature. The first is related to the contact angle a fluid develops at a single confining substrate. The second is based upon virial coefficient information and requires knowledge of the substrate-fluid and fluid-fluid interaction potentials only.     

Simulation of chemical potentials and phase equilibria in two- and three-dimensional square-well fluids: finite size effects

We study the simulation cell size dependence of chemical potential isotherms in subcritical square-well fluids by means of series of canonical ensemble Monte Carlo simulations with increasing numbers of particles, for both three-dimensional bulk systems and two-dimensional planar layers, using Widom-like particle insertion methods. By estimating the corresponding vapor/liquid coexistence densities using a Maxwell-like equal area rule for the subcritical chemical potential isotherms, we are able to study the influence of system size not only on chemical potentials but also on the coexistence properties. The chemical potential versus density isotherms show van der Waals-like loops in the subcritical vapor/liquid coexistence range that exhibit distinct finite size effects for both two- and three-dimensional fluids. Generally, in agreement with recent findings for related studies of Lennard-Jones fluids, the loops shrink with increasing number of particles. In contrast to the subcritical isotherms themselves, the equilibrium vapor/liquid densities show only a weak system size dependence and agree quantitatively with the best-known literature values for three-dimensional fluids. This allows our approach to be used to accurately predict the phase coexistence properties. Our resulting phase equilibrium results for two-dimensional square-well fluids are new. Knowledge concerning finite size effects of square-well systems is important not only for the simulation of thermodynamic properties of simple fluids, but also for the simulation of models of more complex fluids (such as aqueous or polymer fluids) involving square-well interactions.     

A crossover SAFT-VR equation of state for pure fluids: preliminary results for light hydrocarbons

SAFT is perhaps the most versatile, fundamentally, based engineering equation of state in use today. However, in common with all analytic equations of state, SAFT exhibits classical behavior in the critical region rather than the non-analytical, singular behavior seen in real fluids. Recently, so-called crossover equations of state have been developed which solve this shortcoming by incorporating the scaling laws valid asymptotically close to the critical point while reducing to the original classical equation of state far from the critical point. We have combined the SAFT-VR equation of state with an analytical crossover technique to obtain the SAFT-VRX equation of state. The SAFT-VRX approach combines the accurate low temperature behavior of SAFT-VR with a precise representation of the critical region. Preliminary results are presented for hydrocarbon systems which illustrate the accuracy of the SAFT-VRX approach over the entire fluid phase region.