Intramolecular bonding in a statistical associating fluid theory of ring aggregates

The first-order thermodynamic perturbation theory of Wertheim (TPT1) is extended to treat ring aggregates, formed by inter- and intramolecular association. The expression for the residual association contribution to the Helmholtz free energy for ring aggregates, incorporating the appropriate terms in Wertheim's fundamental graph sum of the TPT1 density expansion, is derived to calculate the distribution of the molecular bonding states. This requires the introduction of two new parameters to characterise each possible ring type: the ring size τ, which is equal to one in the case of intramolecular association, and a parameter W that captures the likelihood of two ring-forming sites bonding. The resulting framework can be incorporated in equations of state that account for the residual association contribution to the free energy, such as the statistical associating fluid theory (SAFT) family, or the cubic plus association (CPA) equation of state. This extends the applicability of these equations of state to mixtures with an arbitrary number of association sites capable of hydrogen bonding to form intramolecular and intermolecular rings. The formalism is implemented within SAFT-VR Mie to calculate the fluid-phase equilibria of model chain-like molecules containing two associating sites A and B, allowing for the formation of open-chain aggregates and intramolecular bonds. The effect of adding a second component that competes for the association sites that mediate intramolecular association in the chain is also examined. Accounting for intramolecular bonding is shown to have a significant impact on the phase equilibria of such systems.     

A New Association Term for SAFT Equation of State

This work is devoted to present a new expression for association part of SAFT equation of state. This expression is obtained by calculating the mole fraction of not bonded molecules through a general chemical theory of association. The new equation of state is named CSAFT and has one adjustable parameter less than original SAFT for pure associating fluids. The performance of CSAFT is investigated by PVT calculations for pure fluids and LLE calculations for binary and ternary mixtures. Results show that CSAFT correlates vapor pressure and liquid density data of pure associating components more accurate than SAFT. Also the accuracy of CSAFT model is much better than SAFT for LLE prediction of binary and ternary mixtures, considering SAFT has one pure component parameter more than CSAFT.     

Mutual solubilities of hydrocarbons and water with the CPA EoS

The knowledge of hydrocarbon/water phase equilibria is important in the design and operation of equipment for petroleum transport and refining and petrochemical plants. The presence of water in a hydrocarbon mixture can affect the product quality and damage the operation equipment due to corrosion and formation of gas hydrates. Tracing the concentration of hydrocarbons in aqueous media is also important for technical purposes like preventing oil spills and for ecological concerns such as predicting the fate of these organic pollutants in the environment.     

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

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], [7], [8] and [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₀) and underestimates the first order contribution (A₁) of Helmholtz energy relative to simulation. Nevertheless, the correlations are qualitatively reasonable. Significant inconsistencies arise when considering the second order contribution (A₂). For example, the PC-SAFT characterization of A₂ becomes larger than A₁ 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.     

pePC-SAFT: Modeling of polyelectrolyte systems: 1. Vapor–liquid equilibria

A molecular thermodynamic model for polyelectrolyte systems—called pePC-SAFT—is proposed. The effect of charged monomers within the polyelectrolyte chain is explicitly taken into account in the reference term by replacing the hard-chain contribution of the PC-SAFT model by a charged-hard-chain contribution. Moreover, counterion condensation is accounted for to determine the effective number of charges along the polyion as well as of free counterions. The electrostatic contribution of the free counterions is described by a Debye–Hückel term.     

Vapor liquid equilibrium modeling of alkane systems with Equations of State: “Simplicity versus complexity”

Two types of Equations of State (EoS), which are characterized here as “simple” and “complex” EoS, are evaluated in this study. The “simple” type involves two versions of the Peng–Robinson (PR) EoS: the traditional one that utilizes the experimental critical properties and the acentric factor and the other, referred to as PR-fitted (PR-f), where these parameters are determined by fitting pure compound vapor pressure and saturated liquid volume data. As “complex” EoS in this study are characterized the EoS derived from statistical mechanics considerations and involve the Sanchez–Lacombe (SL) EoS and two versions of the Statistical Associating Fluid Theory (SAFT) EoS, the original and the Perturbed-Chain SAFT (PC-SAFT).

The evaluation of these two types of EoS is carried out with respect to their performance in the prediction and correlation of vapor liquid equilibria in binary and multicomponent mixtures of methane or ethane with alkanes of various degree of asymmetry. It is concluded that for this kind of systems complexity offers no significant advantages over simplicity. Furthermore, the results obtained with the PR-f EoS, especially those for multicomponent systems that are encountered in practice, even with the use of zero binary interaction parameters, indicate that this EoS may become a powerful tool for reservoir fluid phase equilibria modeling.     

Development of intermolecular potential models for electrolyte solutions using an electrolyte SAFT-VR Mie equation of state

ABSTRACT

We present a theoretical framework and parameterisation of intermolecular potentials for aqueous electrolyte solutions using the statistical associating fluid theory based on the Mie interaction potential (SAFT-VR Mie), coupled with the primitive, non-restricted mean-spherical approximation (MSA) for electrolytes. In common with other SAFT approaches, water is modelled as a spherical molecule with four off-centre association sites to represent the hydrogen-bonding interactions; the repulsive and dispersive interactions between the molecular cores are represented with a potential of the Mie (generalised Lennard-Jones) form. The ionic species are modelled as fully dissociated, and each ion is treated as spherical: Coulombic ion–ion interactions are included at the centre of a Mie core; the ion–water interactions are also modelled with a Mie potential without an explicit treatment of ion–dipole interaction. A Born contribution to the Helmholtz free energy of the system is included to account for the process of charging the ions in the aqueous dielectric medium. The parameterisation of the ion potential models is simplified by representing the ion–ion dispersive interaction energies with a modified version of the London theory for the unlike attractions. By combining the Shannon estimates of the size of the ionic species with the Born cavity size reported by Rashin and Honig, the parameterisation of the model is reduced to the determination of a single ion–solvent attractive interaction parameter. The resulting SAFT-VRE Mie parameter sets allow one to accurately reproduce the densities, vapour pressures, and osmotic coefficients for a broad variety of aqueous electrolyte solutions; the activity coefficients of the ions, which are not used in the parameterisation of the models, are also found to be in good agreement with the experimental data. The models are shown to be reliable beyond the molality range considered during parameter estimation. The inclusion of the Born free-energy contribution, together with appropriate estimates for the size of the ionic cavity, allows for accurate predictions of the Gibbs free energy of solvation of the ionic species considered. The solubility limits are also predicted for a number of salts; in cases where reliable reference data are available the predictions are in good agreement with experiment.     

Partitioning behaviour of organic compounds between ionic liquids and supercritical fluids

Applications and prospects of two-phase, tuneable solvent systems composed of ionic liquids (ILs) and supercritical fluids with an emphasis on supercritical carbon dioxide (scCO(2)) are reviewed. The IL-scCO(2) biphasic systems have increasingly been used in diverse fields of chemistry and technology, and some examples of these applications are mentioned here. Rational design of such applications can obviously benefit from pertinent data on phase equilibria including the partition coefficients of the prospective products and reactants between the two phases. Therefore, a reliable technique to measure the limiting partition coefficients would be of value. Here, the pros and cons of supercritical fluid chromatography in this respect are discussed. An overview of methods for predictive thermodynamic modelling of binary (IL-scCO(2)) and ternary (solute-IL-scCO(2)) systems is also included.     

Modelling the effect of methanol,glycol inhibitors and electrolytes on the equilibrium stability of hydrates with the SAFT-VR approach

In this work we integrate the statistical associating fluid theory for fluids interacting through potentials of variable range (SAFT-VR) into a traditional van der Waals and Platteeuw framework for modelling clathrate hydrates. We incorporate a new water–guest cell potential for the hydrate phase that can be related to the potential adopted in the familiar SAFT-VR equation of state for modelling fluids. We show how the ability of this equation of state to treat a wide range of complex fluids increases the scope of hydrate modelling to incorporate, in a single framework, the presence of various inhibitors (alcohols, glycols) or brines – or, indeed, any fluid for which a model is available (for use within SAFT-VR) or can be conveniently obtained. Agreement with experimental results is good throughout and, in many cases, excellent.     

Physical properties from association models

Equation of state models with a ‘chemical’ contribution that accounts for association and solvation effects, are computationally intensive as they have to solve an internal chemical equilibrium calculation.

Frequently, only the final results of this internal calculation are used in the subsequent evaluation of physical properties from the model. As a consequence, a substantial amount of unneeded work is performed. We show here how the state function minimization in the chemical equilibrium calculation can be utilized to simplify the calculation of physical properties like pressure and chemical potentials and the derivatives of these properties with respect to temperature, volume and composition.